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1.  INTRODUCTION 


This  report  summarizes  a  multi-year  program  to  analyze  upper  atmospheric  radiance  data 
to  deduce  product  distributions,  extract  the  operative  kinetic  mechanisms,  and  quantitatively 
determine  local  conditions  in  the  emission  volume.  Spectral  and  kinetic  models  were  used  to 
analyze  emissions  from  several  molecules  including  NO,  NO+,  ozone,  and  excited  nitrogen. 
Radiances  from  the  mesosphere  both  under  equilibrium  and  non-equilibrium  conditions  were 
analyzed.  The  radiance  changes  produced  by  non-equilibrium  excitation  mechanisms  such  as 
electrons,  kinetic  energy  from  fast  atoms,  and  solar-induced  transients  at  terminator  were  all 
addressed  by  multiple  analyses  during  this  program.  Effects  of  electron  excitation  were  isolated 
via  the  analysis  of  data  from  an  artificial  injection  experiment  and  in  the  aurorally  disturbed  and 
sprite-disturbed  mesosphere.  The  EXCEDE  in  electron  beam  injection  experiment  produced 
enhanced  emissions  from  NO  and  Ozone.  Our  findings  are  presented  in  Section  2.  Mechanisms 
responsible  for  auroral  production  of  NO+  are  presented  in  Section  3.  Because  fast  atoms 
produced  in  photolytic  decomposition  and  electron  impact/recombination  apparently  contribute 
to  observed  radiance  enhancements,  we  performed  laboratory  simulation  experiments  to  quantify 
fast  atom  product  distributions  and  confirm  radiometric  calibrations.  These  efforts  are  described 
in  Section  4.  Ozone  radiance  variability  on  global  spatial  scales  and  at  terminators  was  derived 
from  analysis  of  the  data  from  the  CIRRIS  1A  shuttle-based  experiment  and  reported  in 
Section  5.  We  participated  in  the  spectral  analysis  of  field  data  observing  the  transient  red  sprites 
phenomenon  placing  constraints  on  the  excitation  electron  distribution  based  on  prior  laboratory 
investigations  at  the  Research  Laboratory  at  Hanscom.  Our  efforts  are  summarized  in  Section  6. 
The  CIRRIS  1 A  mission  acquired  the  first  high  quality  infrared  atmospheric  emissions  global 
data  base  from  space  -  the  precursor  for  many  other  orbital  missions.  In  Section  7  we  report  a 
brief  investigation  of  IR  radiances  acquired  during  the  Delta  181  mission  providing  simultaneous 
observation  of  several  oxygen  emission  features  -  providing  the  ability  to  constrain  and  improve 
global  radiance  models.  In  Section  8  we  present  a  summary  of  our  participation  in  the  analysis  of 
the  MSX  satellite  data  base  observing  a  strong  auroral  event. 
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Although  the  results  of  several  of  these  analyses  should  be  viewed  as  preliminary,  taken 
together  they  represent  several  significant  steps  forward  in  our  understanding  of  the  non¬ 
equilibrium  upper  atmosphere. 
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2.  ELECTRON  INDUCED  FLUORESCENCE  OBSERVED  DURING 

THE  EXCEDE  3  MISSION 


The  EXCEDE  in  rocket  experiment  bombarded  the  nighttime  lower  thermosphere  with 
energetic  electrons  at  high  dosing  rates  and  observed  the  resulting  infrared  fluorescence  with  a 
suite  of  instruments  including  a  cryogenic  interferometer.1  Nitric  oxide  fluorescence  from  the 
electron  dosed  upper  atmosphere  was  analyzed  to  extract  excitation  mechanisms.  A  large  beam- 
induced  NO  production  was  observed.  Vibrational  /  rotational  distributions  from  previous 
COCHISE  laboratory  measurements  helped  in  mechanism  investigation.  Across  the  altitude 
region  covered  by  EXCEDE  HI,  different  classes  of  processes  will  contribute  to  the  fluorescence. 
The  rotationally  hot  component  appears  likely  to  arise  from  the  reaction  of  metastable  nitrogen 
atoms  with  02.  Production  of  this  component  decreases  as  dose  rate  increases  due  to  precursor 
quenching.  The  interim  findings  of  this  study  were  reported  at  the  Spring  AGU  Meeting  in  May 
1993.  The  viewgraphs  of  that  presentation  are  included  as  Appendix  1  to  this  report. 

Carrying  this  effort  further,  we  quantified  the  relative  contributions  of  the  vibrational 
state  populations  with  both  excited  and  extremely  rotationally  excited  states.  The  role  of 
metastable  and  translationally  excited  atoms  were  discussed.  A  poster  presentation  entitled 
"Kinetics  and  Mechanisms  for  the  Fragmentation  of  Ro  vibrationally  Excited  Nitric  Oxide  in  the 
EXCEDE  III  Rocket  Experiment"  by  Fraser,  Rawlins,  Green,  Paulsen,  Armstrong,  Upson, 
Blumberg,  and  Murphy  was'  given  at  the  1993  Fall  AGU  Meeting  and  is  included  as  Appendix  2 
to  this  report.  In  addition  we  supported  the  companion  paper  on  this  topic  led  by  Air  Force 
Research  Laboratory  researchers. 

Fourier  transform  spectra  in  the  9  to  1?  /im  spectral  region  above  the  mesopause  exhibit 
the  presence  of  vibrationally  excited  03(v3)  excited  both  by  quiescent  and  electron  initiated 
processes.  We  analyzed  the  fluorescence  near  10  fim  to  extract  electron-initiated  processes  from 
the  natural  background  chemiluminescence.  Populations  of  ozone  v3,v=l-4  were  observed  to  be 
enhanced  by  electron  bombardment.  Secondary  electron  excitation  of  ozone  seems  to  be 
occurring  in  agreement  with  laboratory  measurements.  The  poster  presentation  at  the  1993  Fall 
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AGU  Meeting  "Excitation  of  03(v3)  by  Electron  Bombardment  in  the  EXCEDE  HI  Rocket 
Experiment"  by  Rawlins,  Green,  Paulsen,  Blumberg,  and  Murphy  is  included  as  Appendix  3. 

As  other  analyses  of  the  EXCEDE  HI  rocket  profile,  dosing  rates,  and  atmospheric 
concentrations  arose  from  other  analyses  of  the  EXCEDE  in  mission,  a  coherent  picture  of  the 
ozone  production  emerged.  When  the  electron  beam  is  on,  ozone  emission  from  the  irradiated 
volume  is  more  intense  than  that  from  the  entire  ambient  atmospheric  column,  with  excitation 
into  at  least  four  vibrational  quanta  identified.  Direct  excitation  of  ozone  vibration  by  secondary 
electron  impact  plays  a  significant  role.  In  addition,  during  the  period  with  the  most  intense 
dosing,  ozone  is  formed  via  a  bi -molecular  reaction  of  beam-excited  metastable  precursors.  The 
kinetic  behavior  favors  the  OjfA^.v)  +  02  reaction  to  form  03  either  in  excited  states  or  in 
ground  state  subsequently  efficiently  excited  by  electron  impact.  A  journal  article  on  the  beam- 
produced  excitation  of  ozone  as  observed  in  the  EXCEDE  HI  rocket  experiment  is  included  as 
Appendix  4. 
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3.  VIBRATIONALLY  EXCITED  NO  ION  MISSION  FROM 
THE  AURORALLY  DISTURBED  THERMOSPHERE  AS 
OBSERVED  DURING  THE  ELIAS  MISSION 


The  Earth  Limb  Infrared  Atmospheric  Structure  (ELIAS)  experiment  successfully 
acquired  data  on  the  infrared  spatial  and  temporal  structure  created  by  an  aurora  in  the  polar 
upper  atmosphere.  The  sub-orbital  mission  conducted  in  Alaska  on  19  March  1983  was 
fortunate  to  observe  the  emissions  from  an  IBC  in+  aurora  for  several  hundred  seconds. 
Measurements  of  the  atmospheric  radiances  from  well  defined  spatial  volumes  were  acquired 
with  bandpass  radiometers  alternately  staring  at  a  spatial  position  and  with  the  radiometers 
scanned  over  the  spatial  scene.  Ultraviolet,  and  short-  and  mid-wave  infrared  emissions  were 
analyzed  and  correlated  to  extract  auroral  excitation  of  vibrationally  excited  NO  and  NO+. 
Through  modeling  we  determined  that  NO+(v)is  produced  in  aurora  from  three  potential  sources: 
the  reactions  of  N+  with  02,  of  N2+  with  O,  and  of  N  2P  with  O.  Modeling  indicates  that  the  N2+ 
channel  dominates  NO+  production.  The  auroral  production  efficiency  for  the  NO+  fundamental 
band  was  predicted  to  be  1.35  photons  per  ion  pair  with  the  N2+  +  O  reaction  being  the  dominant 
contributor.  The  modeling  indicated  that  about  1%  of  the  deposited  auroral  energy  appears  as 
NO+  as  opposed  to  4%  for  NO.  A  nonlinear  scaling  with  dose  rate  is  predicted  by  the  model. 
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4.  LABORATORY  STUDIES  IN  SUPPORT  OF  FIELD  OBSERVATIONS 


While  field  observations  provide  the  fundamental  observations  of  key  atmospheric 
processes,  controlled  laboratory  experiments  permit  hypotheses  to  be  tested  and  parameters  to  be 
systematically  varied.  As  part  of  this  program  we  performed  several  investigations  to  aid  in  the 
understanding  of  field  data  observations. 

Atoms  moving  with  significant  energies  (of  several  eV)  can  be  produced  in  the  natural 
atmosphere  by  photolysis,  in  the  aurorally  disturbed  atmosphere  by  electron-impact  followed  by 
dissociative  recombination,  and  as  an  artifact  of  space  based  observatories  arising  from  the 
interaction  of  outgassed  or  reflected  species  with  the  residual  natural  atmosphere.  The  role  of 
high  velocity  (hyperthermal)  atoms  in  enhancing  rate  coefficients,  opening  new  reaction 
pathways,  and  altering  vibrational  distributions  produced  by  chemiluminescent  reactions  was 
investigated  using  the  Physical  Sciences  Inc.  (PSI)  fast  atom  source.  In  particular,  the  role  of  fast 
atomic  oxygen  colliding  with  NO  to  produce  emission  from  highly  excited  NO(vJ)  states  in  the 
thermosphere  is  explored  via  controlled  laboratory  experiments. 

SWIR  measurements  were  performed  under  single  collision  conditions  representative  of 
those  encountered  by  molecules  outgassed  or  produced  on  or  above  spacecraft  surfaces  in  low 
earth  orbit  interacting  with  the  residual,  atomic  oxygen  dominated  atmosphere.  The  reaction  of 
O  with  N2  to  produce  vibrationally  and  rotationally  excited  NO  and  O  was  of  particular  interest 
for  both  spacecraft  interactions  as  a  precursor  for  “shuttle  glow”,  and  as  a  potential  contributor  to 
NO(v)  production  in  the  disturbed  atmosphere.  The  O  +  N2  reaction  is  endothermic  by  3.26  eV 
at  room  temperature.  The  center  of  mass  energy  with  O  moving  at  8  km/s  (orbital  velocity)  is 
just  comparable  to  this  endothermicity.  Thus  surrounding  spacecraft,  the  kinetic  energy  has  the 
potential  for  overcoming  the  reaction  endothermicity.  No  theoretical  predictions  were  available 
for  this  reaction  cross-section.  The  experiments  revealed  that  the  cross-section  for  NO(v) 
production  begins  to  rise  significantly  at  10  km/s  with  the  overtone  to  fundamental  intensity  ratio 
indicating  significant  vibrational  excitation.  Spectrally-resolved  emission  that  was  obtained 
using  a  circularly  variable  filter  substantiates  the  presence  of  highly  vibrationally  excited  NO. 
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Comparisons  with  field  data  were  performed.  The  cross-section  for  this  reaction 
produces  sufficient  NO  in  this  gas-phase  reaction  to  explain  the  shuttle  glow  intensities  via  gas 
phase  NO  production,  followed  by  surface  adsorption,  and  catalytic  recombination  with  O.  This 
cross-section  also  produces  NO(v)  in  good  agreement  with  the  SKIRT  mission  CVF  spectrally 
resolved  NO(v)  observations  (also  on  STS-39).2  Near  field  NO(v)  radiances  were  also  observed 
within  the  CIRRIS  1 A  data  set  at  radiance  levels  comparable  to  the  SKIRT  radiances  and 
predictions  based  on  the  laboratory  data.  The  spectral  distribution  observed  in  the  CIRRIS  1 A 
data  matches  the  laboratory  distribution  within  their  overlap  region  including  clear  indications  of 
rotationally  hot  NO  in  a  non-Boltzmann  distribution  are  produced  in  the  fast  O  +  N2  reaction. 

The  reaction  of  O  with  C02  was  also  investigated  to  determine  the  cross-section  as  a 
function  of  velocity  and  the  favored  reaction  pathway  at  8  km/s:  vibrational  excitation  of  C02(v3) 
or  reaction  to  form  CO(v).  The  cross-section  favored  the  theoretical  prediction  of  Schatz  and 
Redmon. 

The  ability  to  measure  cross-sections  for  these  hypervelocity  collisions  required  accurate 
knowledge  of  the  atomic  oxygen  flux.  A  description  of  the  atomic  oxygen  measurement 
technique  and  the  fast  atom  apparatus  producing  a  high  flux  of  5  to  12  km/s  velocity  atomic 
oxygen  and  its  performance  characteristics  is  given  in  Appendix  5.  While  total  flux  can  be 
estimated  by  inlet  oxygen  mass  consumption  and  flow  modeling,  the  spatial  distribution  in  the 
test  (collision)  region  needs  to  be  accurately  measured  to  permit  absolute  cross-section 
determination.  Silver  oxidizes  very  efficiently  upon  atomic  oxygen  impact:  gaining  mass 
through  oxidation;  changing  from  a  reflector  to  optically  transmissive;  and  changing  from  a 
conductor  to  an  insulator.  Space-based  and  ground  experiments  have  used  resistance  and  optical 
changes  of  thin  silver  films  to  monitor  atomic  oxygen  flux.  We  describe  in  Appendices  5  and  6 
measurements  of  atomic  oxygen  fluence  using  a  thin  silver  film’s  resistance  change  (using  an 
actinometer)  and  mass  gain  (using  a  Quartz  Crystal  Microbalance).  The  temperature  dependence 
of  the  silver  reaction  rate  was  measured.  Diffusion  of  O  into  the  silver  film  was  determined  to 
produce  non-linear  responses  at  film  thicknesses  greater  than  25  nm.  Correlating  silver  with 
carbon  coated  QCMs  permitted  absolute  fluence  determination.  We  also  report  of  lower  velocity 
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(2  to  4  km/s)  beams  created  to  study  the  interactions  and  signatures  of  plume  exhaust  species. 
Beams  of  C02,  CO,  and  H20  and  mixtures  were  successfully  formed. 

Metastable  nitrogen  atoms  in  the  2P  and  ^  states  are  important  precursors  to  NO  chemi¬ 
luminescence  in  the  electron  dosed  upper  atmosphere,  as  observed  in  auroral  field  missions  and 
in  artificial  auroral  experiments  such  as  EXCEDE.  Remote  measurement  of  N  C2!*)  concen¬ 
trations  relies  upon  knowledge  of  the  Einstein  coefficient  for  radiative  decay  to  turn  radiances 
into  absolute  column  densities.  Metastable  N  2P  atoms  have  transitions  that  emit  brightly  in  the 
ultraviolet  at  347  nm  to  the  ground  state  and  in  the  near  infrared  at  1040  nm  to  the  ^  state.  The 
relative  brightness  of  the  *P  emission  relative  to  N2+  B-X  First  Negative  and  N2  B-A  First 
Positive  transitions  has  been  used  to  gauge  the  depth  of  auroral  penetration  and  electron  energy. 
The  N  *P-  4S  transition  probability  has  only  been  estimated  in  two  theoretical  calculations.  We 
performed  laboratory  measurements  in  a  discharge-flow  reactor.  Metastables  were  generated  in  a 
microwave  discharge,  and  N  2P  detected  via  emission  at  346.6  nm,  and  via  VUV  resonance 
absorption  at  174  nm  using  monochromators.  The  visible  monochromator  was  calibrated  for 
absolute  photon  emission  rate  measurements  by  observing  the  intensity  of  the  O/NO  air 
afterglow.  The  N  2P  number  densities  were  determined  by  analysis  of  the  fractional  absorption  at 
174  nm.  The  346  nm  absolute  emission  from  the  same  volume  was  determined  to  scale  linearly 
with  N  2P  concentration,  permitting  determination  of  an  Einstein  coefficient  of  0.0054  s'1  with  a 
statistical  uncertainty  from  all  sources  of  ±  26%.  This  agrees  with  both  theoretical  estimates 
within  our  experimental  determination  error  bounds.  These  measurements  are  described  more 
fully  in  Appendix  7. 

Laboratory  measurements  of  absolute  emission  rates  often  rely  upon  a  standard  well- 
characterized  reaction  to  account  for  optical  collection  efficiencies,  optical  element  trans¬ 
missions,  and  geometric  factors,  and  to  determine  the  spectral  responsivity  variation  with 
wavelength.  The  O/NO  air  afterglow  is  probably  the  most  widely  used  such  calibration  standard. 
We  have  used  it  to  determine  many  absolute  decay  rates  and  reaction  efficiencies  (such  as  for 
N  2P  presented  above).  Upon  critically  reviewing  the  literature  for  the  O/NO  spectral  distribution 
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(and  converting  them  to  comparable  units),  we  discovered  significant  variations  at  wavelengths 
longer  than  700  nm. 

We  performed  careful  measurements  of  the  O/NO  intensity  distribution  over  the  400  to 
1750  nm  spectral  region.  We  generated  the  precursors  using  multiple  independent  techniques, 
and  calibrated  the  spectral  response  of  the  spectrometer/photomultiplier  using  several  different 
radiation  standards  with  calibrations  traceable  to  NIST.  We  find  the  peak  of  the  distribution  to 
be  around  670  nm.  Detailed  comparison  with  the  distributions  determined  by  eight  other  groups 
is  presented  in  Appendix  8  along  with  our  tabulated  determination.  Although  the  agreement  is 
often  excellent,  significant  discrepancies  occur. 

We  also  provide  a  re-analysis  of  their  data  to  interpret  the  differences.  Again  comparison 
of  the  absolute  rate  coefficient  reported  by  the  eight  groups  is  confounded  by  inconsistent 
reporting.  Careful  reanalysis  and  explanation  of  their  data  produces  a  consistent  value  for  five  of 
the  groups.  We  recommend  a  value  of  (1.17  ±  0.1 1)  x  10'19  cm3  molecule'1  s'1  nm"1  at  580  nm, 
with  a  value  integrated  over  all  wavelengths  of  9  x  10"17  cm3  molecule'1  s'1. 


9 


5.  OZONE  VARIABILITY  NEAR  THE  MESOPAUSE  FROM  THE 

CIRRIS  1A  DATA  BASE 


The  CIRRIS  1 A  mission  acquired  the  first  global  infrared  atmospheric  emissions  data 
base  during  a  several  day  mission  in  1991.3  They  successfully  obtained  over  30  hours  of 
observation  time  of  emissions  from  the  daytime  and  nighttime  quiescent  and  aurorally  disturbed 
mesosphere  and  thermosphere  using  a  cryogenically-cooled,  telescoped  interferometer  and 
radiometers.  In  earlier  work  we  quantified  the  optical  contamination  of  CIRRIS  data  by  near 
field  particles  and  high-velocity  gas  phase  interactions.  As  part  of  this  program  we  investigated 
upper  atmospheric  ozone  radiant  signatures  obtained  with  the  CIRRIS  1 A  interferometer. 

First  data  from  a  few  high  quality  observations  were  thoroughly  analyzed  to  identify 
ozone  vibrational  states  contributing  to  the  observed  emission  signatures.  Limb  emission  spectra 
from  three  tangent  point  altitudes  of  71,86,  and  98  km  were  selected  to  span  the  altitude  regions 
of  most  interest.  The  spectroscopic  constants  and  radiative  lifetimes  were  acquired  under 
previous  laboratory  investigations  in  the  COCHISE  Facility  at  the  Research  Laboratory  at 
Han  scorn.4  Both  spectral  smoothing  and  baseline  corrections  were  applied  before  spectral  fitting. 
We  attempted  to  fit  the  LWIR  ozone  emission  between  820  and  1060  cm'1  via  three  approaches: 
including  seven  excited  states  with  only  v3  quanta  excited;  including  seven  v„v3  coupled  dyads 
(101,002;  102,003;...);  and  including  twelve  v,  v3  and  vzv3  coupled  dyad  states  (101,002; 

012,1 11; ...)  in  the  analysis.  We  determined  that  the  spectral  signature  was  well  described  by 
including  the  v„v3  and  v2,v3  dyads  in  the  spectral  fitting.  Because  addition  of  v2,v3  dyad  emission 
does  not  improve  the  quality  of  the  least-squares  agreement  with  the  data,  we  feel  its  inclusion  is 
not  required. 

Hot  band  emission  and  spectral  structure  from  ozone  excited  states  was  detected  out  to 
beyond  12  /im,  demonstrating  the  excellent  quality  of  the  data  base.  The  vibrational  population 
distribution  observed  agreed  well  with  modeling  predictions  and  the  distribution  derived  from  the 
earlier  SPIRE  mission  data.  The  ozone  spectral  distribution  changed  with  altitude  indicating 
vibrational  level  dependent  quenching  occurs.  The  preliminary  quenching  coefficients  for  v^  7 


10 


agree  with  both  LAB  CEDE  laboratory5  and  SPIRE6  field  observations.  A  preliminary 
examination  of  ozone  4.8  /tm  emission  supported  these  conclusions.  These  findings  were 
summarized  at  the  CIRRIS  1A  Program  Review  on  September  9, 1993,  and  appeared  in  the 
Proceedings  of  that  meeting.  The  viewgraph  presentation  material  in  given  in  Appendix  9. 

Spectral  fitting  of  additional  selected  high  resolution  data  from  the  wider  CIRRIS  1A  data 
base  permitted  internal  distributions  to  be  extracted.  Emission  from  (v„v3)  and  v3  states  were 
observed  in  the  data.  The  CIRRIS  1A  data  was  of  excellent  quality,  however,  emission  did  not 
extend  beyond  the  laboratory  data  base.  No  new  spectroscopic  information  on  ozone  highly 
excited  vibrational  states  could  be  derived.  The  observed  distributions  are  consistent  with  three- 
body  recombination  as  the  source  of  ozone  chemiluminescence.  While  the  CIRRIS  data  agreed 
well  with  SPIRE  and  EXCEDE  El  background  radiances,  it  highlighted  the  need  for  improved 
ozone  photochemical  models 

Armed  with  this  knowledge  of  the  ozone  excited  states  contributing  to  the  observed 
radiant  signatures,  we  next  undertook  a  screening  of  the  entire  CIRRIS  1 A  data  base  to  extract 
variabilities  due  to  spatial,  latitudinal,  and  solar  effects.  With  the  assistance  of  data  base 
screening  tools  developed  by  coworkers  at  the  Air  Force  Research  Laboratory,  we  retrieved  the 
LWIR  ozone  data  “scans”  at  three  spectral  resolutions  for  all  observations  with  tangent  point 
altitudes  from  70  to  1 10  km  and  transferred  the  data  to  PSI  for  analysis. 

Initial  batch  processing  of  the  data  sets  was  performed  to  produce  radiance  altitude 
profiles.  These  profiles  did  not  follow  a  smooth  altitude  scaling.  The  profiles  exhibited  wide 
scatter  between  the  data  with  radiances  from  many  scans  an  order-of-  magnitude  larger  than 
expected.  Upon  visual  inspection,  LWIR  spectra  in  the  9  to  12  pun  region  of  the  outlying  data 
points  were  observed  to  be  not  ozone,  nor  even  resemble  molecular  emission.  We  performed  the 
arduous  task  of  visually  reviewing  each  of  4000  spectra  in  the  entire  data  base  to  identify  and 
remove  spectral  scans  containing  artifacts.  These  are  likely  arise  from  either  interferometer 
scanning  anomalies  or  from  shuttle  environmental  perturbations  (such  as  thruster  firings  required 
to  maintain  pointing  accuracy)  affecting  the  ability  to  observe  the  far-field  atmospheric  signature. 
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Over  32%  of  the  spectra  were  eliminated  via  this  screening  process,  most  for  containing  excess 
noise.  These  removed  spectra  had  the  wrong  spectral  shape,  high  or  low  frequency  oscillations 
running  across  the  spectra,  baseline  offsets,  or  other  clear  indications  that  the  spectra  were 
corrupted  and  not  representative  of  remote  mesospheric  emissions.  Due  to  the  regular  periodicity 
of  the  “noisy”  spectra,  we  feel  these  are  contaminated  radiance  from  the  shuttle  thrusters  rather 
than  interferometer  scan  anomalies. 

The  remaining  2740  spectra  were  too  numerous  to  permit  detailed  spectral  fitting.  To 
characterize  the  large  data  base,  we  integrated  the  spectra  over  ten  bandpasses  within  the  ozone 
feature  to  permit  trends  in  v3,v=l,  hotband,  and  C02  radiances  to  be  extracted.  The  bandpass 
wavelengths  were  selected  based  on  our  spectral  fitting  analysis  to  isolate  the  ozone  v=l  R 
branch,  v=2,  and  other  ozone  hotbands  emissions,  while  carefully  accounting  for  baseline  shifts 
and  ramps  and  spectral  interferences  from  mesospheric  C02  emission  under  sunlit  conditions. 

The  integrals  were  used  to  correlate  ozone  relative  excited  state  densities.  The  integrals  from  the 
2740  spectral  scans  were  incorporated  into  a  single  spreadsheet  and  manipulated  to  display 
profiles  as  a  function  of  altitude,  latitude,  solar  zenith  angle,  and  longitude. 

Our  review  of  this  data  base  indicates  that  the  global  variability  is  large  compared  to  the 
local  variability.  Over  an  order-of-magnitude  variability  was  observed  in  the  ozone  v3,v=l 
radiances  (column  densities)  at  altitudes  above  100  km  even  after  background  correction.  Ozone 
fundamental  band  radiance  was  observable  above  the  noise  for  tangent  point  altitudes  up  to 
109  km.  Excellent  quality  terminator  series  were  obtained  after  correction  for  shuttle  pointing 
jitter.  Both  1 -quanta  states  (proportional  to  total  ozone  concentration)  and  3-quanta  states 
(proportional  to  atomic  oxygen  concentration)  showed  clear  terminator  effects  even  up  to  100  km 
tangent  point  altitudes.  The  atmospheric  C02  radiances  were  used  to  discriminate  local 
concentration  variances  from  larger  scale  atmospheric  transport  (heave).  Ozone  variability  was 
far  in  excess  of  C02  variability.  The  ozone  variability  does  not  seem  to  be  due  to  kinetic  rate 
uncertainties  but  rather  due  to  dynamic  processes.  A  very  significant  local  enhancement  in  ozone 
and  atomic  oxygen  was  observed  later  in  time  and  equatorward  of  an  intense  auroral  arc, 
suggesting  large  scale  atmospheric  air  parcel  transport.  Ozone  variability  was  observed  to  be 
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greatest  during  high  Kp  periods.  There  were  no  clear  trends  with  latitude  or  solar  zenith  angle 
within  the  day  and  nightside  data  sets,  but  substantial  scatter  existed  between  the  values  from 
different  data  collection  times.  The  transition  from  day  to  night  distributions  appears  to  occur 
rather  rapidly.  Crude  temperatures  were  obtained  from  bandpass  ratios.  These  results  were 
presented  at  the  21st  IUGG  Meeting  in  July  1995  using  the  material  in  Appendix  10.  A 
preliminary  draft  of  a  journal  article  summarizing  these  analysis  efforts  is  presented  in 
Appendix  11. 
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6.  SPECTRAL  ANALYSIS  OF  THE  EMISSION  FROM  RED  SPRITES 


During  the  course  of  this  program  a  new  class  of  atmospheric  phenomena  was  first 
systematically  observed:  sprites.  These  transient  spatially  structured  radiant  columns  connect  the 
tops  of  intense  thunderstorms  to  the  bottom  of  ionosphere  and  thus  are  felt  to  play  a  significant 
role  in  the  global  electrical  circuit.  The  first  spectrally  resolved  observations  of  sprites  exhibited 
a  visible  spectra  dominated  by  nitrogen  First  Positive  emission.  PSI  has  extensive  experience  in 
the  analysis  of  radiance  from  the  N2  B-A  system  to  extract  the  excited  state  populations  present 
in  the  excited  volume.  Using  PSI  internal  funding  we  supported  meetings  sponsored  by  AFRL 
and  performed  a  spectral  fitting  analysis  on  the  first  spectrally  resolved  sprite  spectrum  to  deter¬ 
mine  the  B  state  distribution  in  the  red  sprite  field  observations.  Emission  from  570  to  820  nm 
exhibited  N2  B-A  in  the  Av=2,3,4  sequences  from  vibrational  levels  up  to  v=l  1.  N2+  A-X  Meinel 
band  emission  was  not  observed,  only  an  upper  bound  on  its  column  density  could  be  set. 
Assumption  of  a  Boltzmann  distribution  for  the  electron  energy  permitted  the  B-state  vibrational 
distribution  to  be  predicted  as  a  function  of  electron  temperature  based  on  previous  laboratory 
measurements.  We  deduced  that  electrons  with  a  characteristic  temperature  of  1  eV  were 
producing  the  observed  sprite  distribution,  and  that  no  ion  state  emission  was  expected  (nor 
observed)  at  these  electron  energies.  This  work  was  completed  under  AFRL  sponsorship  and  the 
complete  text  of  a  Geophysics  Research  Letters  article  is  included  as  Appendix  12. 

As  part  of  this  program  we  undertook  a  more  systematic  fitting  of  sprite  data  acquired  by 
the  University  of  Alaska  group  during  the  1995  summer  campaign.  We  determined  the  internal 
vibrational  and  electronic  state  energy  distributions  of  the  large  spectrometric  observations  data 
base  to  explore  red  sprite  variability.  This  data  set  exhibited  a  very  large  signal  to  noise  dynamic 
range  of  over  1000.  Only  small  variations  in  the  spectral  distributions  were  obtained  for  all 
sprites  analyzed.  The  derived  B-state  distribution  agreed  with  that  extracted  from  the  earlier 
analysis  of  the  first  sprite  spectrum.  These  findings  were  summarized  in  a  presentation  at  the 
Fall  AGU  Meeting  which  is  included  as  Appendix  13. 
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This  analysis  was  extended  to  predict  the  likely  infrared  radiances  within  the  red  sprite 
column.  We  considered  several  classes  of  processes  that  could  produce  IR  emission.  Electron 
impact  excitation  of  Nitrogen  electronic  states  can  directly  produce  W-B  emission  as  observed  in 
EXCEDE5  and  COCHISE7  experiments.  Based  on  observed  600  kR  B-state  radiances,  emission 
from  the  nested  W  state  can  be  accurately  predicted  as  only  1  kR  around  3  /xm.  Electron  impact 
dissociative  excitation  of  N  atom  metastable  states  to  produce  NO  is  predicted  to  yield  small 
MWIR  radiances  except  for  long-lived  sprite  columns.  Electron  impact  efficiently  produces 
vibrationally  excited  N2.  In  the  mesosphere  vibra-luminescent  near-resonant  energy  transfer 
from  N2(v)  to  C02(v3)  is  the  only  significant  relaxation  mechanism.  We  presented  these  findings 
at  a  Sprites  Workshop  at  Research  Laboratory,  Han  scorn  in  April  1997.  The  viewgraphs  from 
this  presentation  are  included  as  Appendix  14  of  this  report. 

We  next  undertook  more  detailed  modeling  of  the  electron  distribution  including  all 
excitation  pathways  to  predict  absolute  excitation  rate  as  a  function  of  atmospheric  field  strength 
E/N.  Energy-dependent  cross-sections  for  vibrational  levels  in  all  related  nitrogen  electronic 
states  were  assembled  as  input  into  the  ELENDIF  Boltzmann  transport  model.8  Excitation  into 
each  excited  state  of  nitrogen  and  oxygen  was  computed  as  a  function  of  E/N.  We  learned  that 
the  observed  sprites  N2B  state  distribution  favors  E/N  in  the  3  to  5  x  10‘16  V/cm  range.  Moreover 
the  ion  states  are  predicted  to  be  formed  at  completely  negligible  levels.  UV  emissions  from  the 
long-lived  sprite  component  are  expected  small  and  due  to  N2  C  fluorescence.  Any  UV  ion 
emissions  must  arise  from  the  short  lived  (<  1  ms)  precursor  pulse. 

The  calculated  rates  from  the  code  permit  the  required  electron  density  to  be  estimated 
from  the  observed  sprite  intensities  as  a  function  of  E/N.  E/N  of  3  to  5  x  10'16  require  electron 
densities  in  the  range  of  104/cc,  a  value  in  agreement  with  measurements.  Thi  s,  we  feel  the 
physical  regime  where  sprites  can  exist  are  tightly  bounded  by  the  current  visible  observational 
data.  This  material  was  presented  at  a  Sprites  Meeting  at  Air  Force  Research  Laboratory  at 
Hanscom  in  June  1998  and  is  included  as  Appendix  15. 
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We  feel  these  observations  should  help  constrain  sprite  mechanisms  and  provide 
guidance  into  formation  processes.  Based  on  these  conclusions,  we  feel  the  emphasis  for  future 
measurements  should  be  on  performing  1R  measurements.  A  draft  journal  article  describing 
these  results  is  has  been  prepared  for  journal  submission  and  is  included  as  Appendix  16  of  this 
report. 
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7.  DIURNAL  VARIABILITY  OF  INFRARED  ATMOSPHERIC  OXYGEN 
EMISSIONS  AS  OBSERVED  FROM  THE  DELTA  181  MISSION 


The  Delta  181  satellite  was  a  quick  response  multipurpose  mission  that  operated  for 
several  days  in  a  low  inclination  orbit  at  about  300  km  altitude.  Two  cryogenically  cooled  CVF 
spectrometers  providing  continuous  spectral  coverage  from  700  to  22000  nm  (visible  to  LWIR) 
were  part  of  the  instrumentation  package.  Several  observations  of  earthlimb  emissions  were 
successfully  completed  during  the  brief  mission.  Daytime,  nighttime  and  dawn  and  dusk  term¬ 
inator  emissions  were  observed  with  full  spectral  coverage.  Although  the  CVF  spectrometers 
had  poor  off-axis  rejection  and  a  large  footprint  at  the  tangent  point,  this  data  base  offers  the 
unique  opportunity  to  understand  oxygen  chemistry  near  the  mesopause.  The  Delta  181  data 
base  provides  a  simultaneous  observation  of  ozone  v3  LWIR  emission  (proportional  to  ozone 
density),  of  02  ‘A  emission  at  1.58  /im  (proportional  to  ozone  during  the  day),  and  OH  vibrational 
bands  (proportional  to  H  and  ozone).  Surprisingly,  neither  CIRRIS  1 A  nor  MSX  have 
instrumentation  with  the  ability  to  perform  these  observations.  Although  we  had  reduced  and 
analyzed  the  data  as  part  of  a  previous  program,  we  undertook  comparison  to  photochemical 
models  and  previous  field  radiance  observations  as  part  of  this  program.  The  Delta  181  data 
again  emphasizes  that  atmospheric  ozone  photochemical  models  appear  unable  to  accurately 
predict  nighttime  ozone  concentrations.  This  brief  effort  is  summarized  in  a  draft  article  for 
journal  submission  and  included  as  Appendix  17. 
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8.  MSX  EARTHLIMB  ANALYSES  AND  SUPPORT 


The  Midcourse  Space  Experiment  provided  an  excellent  multispectral  earthlimb  data  base 
during  the  one  year  long  cryogenic  phase  of  its  mission.  Calibrated  spectrometers  and  imagers/ 
radiometers  provided  nearly  continuous  spectral  coverage  over  the  110  to  26000  nm  region.9  PSI 
provided  technical  assistance  to  the  MSX  Earthlimb  team  during  the  course  of  this  program.  The 
work  involved  coordinating  and  analyzing  data  from  several  Earthlimb  experiments.  Particular 
attention  was  paid  to  interpreting  the  results  of  the  aurora  of  10  November  1996  that  was 
observed  over  Scandinavia.  Aurora  are  important  upper  atmospheric  excitation  phenomena  with 
critical  systems  implications.  The  Midcourse  Space  Experiment  provided  the  first  space-based 
observations  of  this  phenomena  with  co-aligned  and  time  synchronized  instruments  covering 
wavelengths  from  the  hard  ultraviolet  into  the  long-wave  infrared. 

The  aurora  was  bright,  measuring  in  excess  of  200  kR  in  the  3914  A  band,  and  was 
observed  for  a  total  of  19  minutes.  The  field-of-view  of  the  sensors  was  maintained  at  a  constant 
tangent  height  of  105  km.  The  UVISI  ultraviolet  and  visible  imagers  acquired  both  narrow  and 
wide-field  data  at  a  rate  of  two  frames  per  second.  SPIRIT  El  acquired  data  in  all  five  IR  bands 
with  the  minor  scan  on  creating  3  degree  images.  All  five  of  the  UVISI  spectrographic  imagers 
were  employed  as  was  the  space-based  visible  camera. 

Substantial  auroral  enhancements  were  observed  in  the  ultraviolet  and  visible  wave¬ 
lengths.  Auroral  effects  were  also  observed  in  three  of  the  SPIRIT  HI  infrared  radiometer  bands. 
SPIRIT  IH  observed  a  strong  4.3  fim  enhancement  of  10-20  times  the  ambient  background  above 
100  km  with  smaller  enhancements  at  lower  altitudes  (a  fourfold  enhancement  at  95  km  and  a 
two-fold  enhancement  at  90  km).  The  strongest  dosing  regions  in  the  UV  and  IR  appear  to 
coincide.  A  video  presentation  from  all  primary  sensors  [SBV  (visible),  SPIRIT  III  (IR 
radiometer),  and  UVISI  (images  and  spectral  resolved  images)]  clearly  displayed  this  correlation. 
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A  preliminary  modeling  evaluation  determined  that  the  strong  C02  enhancements  at 
lower  altitudes  requires  300  to  500  seconds  pre-dosing  (confirmed  in  late  time  observations). 
Preliminary  AARC  model  comparisons  produce  a  consistent  altitude  profile  with  the  enhanced 
C02  data.  The  4.3  fim  enhancement  is  consistent  with  known  sources  and  excitation  mechan¬ 
isms.  The  principal  contributor  is  believed  to  be  v-v  transfer  from  vibrationally  excited  N2. 
Another  contributor  to  this  wavelength  region  is  NO+  emissions  arising  from  ion-molecule 
reactions. 

The  enhancement  in  Band  A  (6  to  10.9  /xm)  is  attributed  to  NO  radiance  above  ambient 
levels  produced  from  Nf^D)  +  02.  The  IR  enhancements  are  dynamic  with  rapid  temporal 
variations  and  strong  gradients. 

This  work  was  presented  at  the  Spring  AGU  Conference  in  May  1997  and  at  the  1997 
Meeting  of  the  MSX  Interim  Results  Review  (IRR)  in  June  1997.  The  viewgraphs  from  these 
presentations  are  included  as  Appendices  18  and  19  respectively.  Future  analyses  will  rely  upon 
insights  gathered  from  sub-orbital  and  early  orbital  data  bases,  but  the  depth  and  spatial/temporal 
coverage  of  MSX  data  will  supply  the  most  useful  data  sets  for  model  refinement  and  validation 
for  many  years  to  come. 
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APPENDIX  1 


Dosing.  Altitude  Dependence,  and  Excitation  Mechanisms  of  NO  Vibraluminescence 

Observed  in  the  EXCEDE  III  Mission 
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Sponsored  by  the  Defense  Nuclear  Agency 


EXCEDE  III  NITRIC  OXIDE  ANALYSIS  PLAN 
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Scalings  provide  insight  into  excitation  mechanisms  and  precursor  fates 


NO  EXCITATION  MECHANISMS 
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•  Possible  but  never  observed 
O(fast)  +  NO,  N2 

•  Experiments  in  progress 


VIBRATIONAL  DISTRIBUTIONS:  UPLEG  LOW-J 


Vibrational  Level 
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DOSING  EFFICIENCY:  THERMAL  COMPONENT 


KEY  PROCESSES  FOR  EXCEDE 
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N0(v,J)+M  J  cooling,  J  — »  V? 
Radiative  cascade 

Transport  controls  yield,  distributions  for 
apogee  and  downleg 


EXCEDE  III  (NO)v  LOSS  RATES 
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Rate  (s*1) 


PRECURSOR  REACTION  RATES 
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CONCLUSIONS  AND  REMAINING  ISSUES 
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Hot  N*  may  affect  distributions,  production  rates 


APPENDIX  2 


Kinetics  and  Mechanisms  for  the  Formation  of  Rovibrationallv-Excited 
Nitric  Oxide  in  the  EXCEDE  in  Rocket  Experiment 
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Sponsored  by  the  Defense  Nuclear  Agency 

December  1993 
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modeling  calculations  suggest  that  rotationally/vibrationally  excited  NO(v,J)  ii 
rapidly  formed  at  high  altitudes  by  reactions(s)  of  translationally  hot  N(2D 
and/or  4S)  with  O2.  We  discuss  the  possible  excitation  reactions  in  terms  of 
the  kinetic  constraints  imposed  by  the  spectral  measurements 
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2  cnr1  cryogenic  interferometer  views  excitation  volume 
Copious  NO(v, J)  excitation  95  to  1 1 5  km 


EXCEDE  III  NITRIC  OXIDE  ANALYSIS 
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Scalings  provide  insight  into  excitation  mechanisms  and  precursor  kinetics 
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Intensity  (W/cm2  sr 
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Short  residence  time  controls  yield,  distributions  for  apogee  and  descent 
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source  of  N0(v,high-J) 
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Collisions  reduce  role  for  lower  altitudes  and  longer 
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NO(v)  PRODUCTION  EFFICIENCIES:  ASCENT 

Comparison  To  Model  Predictions 
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NrD,  ^S)  dominates  N0(v,J)  production  at  high  altitudes 
— >  relevant  to  dayglow,  auroras  above  120  km 


APPENDIX  3 


Excitation  of  P3(v3)  by  Electron  Bombardment  in  the  EXCEDE  HI  Rocket  Experiment 

Presentation 


73/74 


EXCITATION  OF  Oa(vr»)  BY  ELECTR 


T- 18795 


■£  © 
i  c 

CD  g 
C/D' 


CD 


0 

i— 

O 

O  ©  O  § 


0 


0  0 


=  5  0  0  S 


0  _  u>  vj  ^  7= 

a.fc  E  Q  ®  E  2 

0  3° 

Oqzo 

E  -o  ° 

©fi^aS 

0  © 

•4— »  l. 


0 


0 

0 

0 


O 

Q.  0 


>  C 


f  i 

C 

©  0 

O  0 

O-jc 
E 
o 
o 


c 

o 


E$Ep| 


0  3-^0  £  ^ 
(1)  •  *  TT  O) aM 

®  ©  $  2  0  ©  © 


0  .i=  O  _  ~  _  h-  « 

I  5>g  It®  2--S 


o 

©  TD  0 
£0 

o  J?  ©  — 

2.0  8.-S 
.i®©g 

c  E  E  © 

O  Q-o ^ 

i=M-  0  0 
°  •—  -^ 
0  o  o  = 


JC  “  co*  0  CD  x  D)C  O  ©  g  0-0)# 

•2>®  J=g-oJ=  ®  =  «  o  ®  ®  ®.g  g 

C£  g?c“®s?E "£ >  i 


®"5  o.-a  ® 

£®“cc2®®S® 

C 


^  o  .  -  O 


_  £ 

m  © 
©  0 


§  § 


0 

c  TO  0  ^  -i=  CD  *>  —  ■*-  uj 

im©o0O-5CO-0  co©  ©  ©  co 


0 


>>ja  > 


>  © 
0  cr 
>•—•*- 


0 


*  e  a-S 


d)5“  So  =_ 


(B  °  »  C3  g  c  »  <D  g’S.-K^  C  0  = 

—  ©©crEo0"HJ950c0©:0LLi 

~  ©  O)  0  £  TJ  m-—  O  V_  Q)1-1 
O  ^ ^  v  o  C/) 

O  n- o  0  0  0 

<-r-LJ  t—  cn  -*— 


E  c  2  ©  ~  E 
o  © ~ ££  © 
-Q  ©  0  0  S 

0  -Dif  _Q  -Q 


©.^0-0  £,p  ©  f0  o 

q0  0  0  coE  •—  0  £  0 

0S-t.£0  s  £2  £  ®  D 

®oo£?s- 

r® 8|g8 

<B-°  ®  c  ^  ® 


f-  w 

0-0  oo 

= soils' 

y£.9  -§±:'D  —  Oc  — 

Q^CCXo-CO^mO0 
LL|00O©0-^0©  q.*£  - 
Om  O^OD>>a3  E  0  0  0- 

S^a'O0=^0_  —  7-5 

0  0£ 

C  —  Cl*' 

o 


X  o 

LU  .9 


0  0 


©75 


'W  w 

E  0?  ^ 


0  C  •—  © 


©°  O©  00^.2 


U/  w  w  —  w  ^  - 

0D0.J-+-2  C  0  (0t£  0 


0 


77  ts  s  u/  "s;  w-'  ■*— •  CZ  0  _  w 

ajr  ©  2  2  §0  5-o2.SC  o  o 

c  ~  0.0  f  S>jDr  (D.t:OO0p 
©  ^  0  ’>  ^  .E  .©  ’>  0  '>*£  DD  0  Q.  Q. 


76 


ELECTRON  EXCITATION  OF  03(v) 


CO 

h- 

CD 


CO 

CL 


0 

c 

0 

E 

0 

i— 

3 

0 

0 

0 

E 

0 

g  i$ 

.2 

IS  ®  7= 

0  v*  Z 
>  O  w— 

»-  T5 


0 
( f. ) 
JD 
O 


Ol 


CD 

CD 


0 


C 

0 


■55  o 
>  c 


W  C  I 

5  °  — 

O  -^=;  CD 

>  -225 
£  *5® 

CL  X  3- 
0  _ 

°  ib 
g*  oCC 

£  DO. 

|  <c o 

3 
CO 


Q 

CD 

0 

C 

0 

£ 

U- 

0 

CL 

X 

0 


O 

■*—> 

0 

O 

.Q 

0 


E 

.x 

o 

o 


E 

O  cc 

CO 

T" 

I 

o 

Al 

D 


0 

c 

0 

0 


0 

w  <*— • 


3 

0 


C  = 
0 


0  o  o 

6  m  £9 


®  c 
0  — 
—  •*-> 

.  O 

0  ® 
+  0 

0 

CD 
1— 

CO 


CO 


o 

t  ft 

CO  0> 

£>  O  c 


0 
+  "O 

•  g 

0  E 


O 

+ 

^co 

o 

1 

'Kj 

O 

+ 


C\J 

O 

0 

0 

0 

O 

O 


0 

O 

£ 

0 

o 

>* 

T3 

O 

-Q 

I 

O 

$ 

+-» 

0 

jD 

‘0 

0 

O 

Q_ 


c 
0 

E 

i_ 

0 

CL 
X 
.  0 

0 
X 
O 
O 

v_ 

E 
0 
0 
-Q 

C 

2| 

8  §. 
Q  O 
_  CD 
=  CD 

LU  Z 

Q  Q. 
< 

Is- 
C\J 

U 
0 


0 

~U 

c 

0 

CO 

0 


E 

0 

0 

X5 

C 

o 

1— 

o 

0 


0 

i_ 

0 

0 
0  2 

0  o 

E  ® 
o  Q- 

0-r 

1* 

oO 

C  vC 
0  0 

g>0 

frg 

O  O 
«•— * 

I® 

18- 

C\I  > 

2=) 

> 

0 
•«— > 

c 

0 

E 
0 


3 
0 
0  0 


0 

i— 

0 

0 

E 

o 

■4— » 

O 


CO 

I 

E 

o 

CD 

o 

T 

0 


0 
CO 

E  ° 

5  8 

2  o 

§.s 

c  w 
o 


0 

E 


LU 

o 

X 

LU 

0 


CL 

E 

< 

i 

co 


Si  ^  > 


0 

E 

0 

O 

c 

0 

O 

0 


ZS  g 
0  ^ 

0  0 


0  0 
O 

LD  3 


CL 

0~ 

»_ 

0 

0 

E 

o 

T3 

0 

i_ 

0 

3 


E 

x 

ID 


o 
■*— » 

LO 
CO 

0 
0 
"D 
3 


0 
0 

£  E 
o  ~ 

T-  0 

x  o 
x  c 
CO  0 
H  -D 

11 
I  s 

^  S 


£ 

0 

2 


—I  C\J  LL  CL  < 


E 

0 

0 

JD 


C 

0 

O 

0 

0 

U 

~u 

c 

0 

C 

0 

O 

0 

0 

E" 

-X 

LO 

o 


LO 

00 

0 

v. 

4— > 

o 

0 

CL 

0 

0 

O 

C 

0 

O 

0 

0 


C 

O 


E 
o  0 

30 


0 

3 

0 

L. 

0 


CO 


02 


E 

x 

o 

o 


0 
0 
C 
■#— < 
c 
0 

E 

0 

O 

c 

0 

c 

0 

"O 

0 

O 

3 

"O 

C 

*  i 

E 

0 

0 

-Q 

* 

C 

0 

0 

CL 

CL 

< 


77 


T- 18797 


0) 


CD 

~  i 

E  ~ 

§  2 

CD  2 

■°  © 
C  -Q 

2  CO 

■5  2 

®  5 

a>  o 

CD  5= 
P  CO 
I—  CO 
O) 

CD  CD 
O  JZ 
C  ■ 


o> 

12 
13  ■*“* 
CO  O 
CD  ~ 
*-  CD 

CO  2 
X  © 

10  w 

l§ 

CD 

X)  CO 
0  "u. 

*=  .c  o 


CD  co 
O  CD 
CO  C 
0  •= 
V-  F 

o  E 

3  0 


CO 

-S 

0 

<: 

o 


0 

.2-0 

0 


o 

0 

> 


o 

o 


co  <d 

_  z  > 

O  o  0 
0.0  o£ 

CO  0  ©  - 

O  >  > 

E  >'^'o> 

0  0  0  5= 

"O  —  —  © 
q)  (D  0  r 
&  >  >  ~ 
2  0  0  0 

?h^.o 

"  S2 


0  CO  m  .E 
0  0  o  0 

-Q  c  (O  E 

-  =  0  O 


X3 

£  12  c 


CO 

0 


(1)  —  0  0 
h=  >r  E 

t  n  »-  **'  •_= 


_ 

0  o  ^  Z  *■*= 

£  "E  llj  °  © 
g  ©mu_  o 

O)  SP©  to  © 

3  r  (I) 

■>  w  o  =  ® 

.2#  c3? 

>  0  o  c  2 

0  LU  w  w 

o  r.2T? 

c  C  •-  1  -I 

•2  o  E  §  3 

s  *  SI  s 

<0  0  -O  0 

—  o  ©  *“  £ 

=  0£  o>0 

=5  .c  _o  _2 


78 
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Abstract 


The  EXCEDE  IH  rocket  experiment  bombarded  the  nighttime  lower  thermosphere  with 
energetic  electrons  at  high  dosing  rates  and  observed  the  resulting  infrared  fluorescence  with  a 
suite  of  instruments  including  a  cryogenic  interferometer.  Fourier  transform  spectra  in  the 
900  to  1 100  cm'1  (9  to  1 1  pm)  spectral  region,  obtained  at  altitudes  between  85  and  105  km, 
show  the  presence  of  vibrationally  excited  03(v3)  excited  by  both  quiescent  (ambient)  and 
electron-initiated  processes.  When  the  electron  beam  is  on,  03(v)  emission  from  the  irradiated 
volume  is  more  intense  than  that  from  the  ambient  background,  and  the  excitation  of  at  least  four 
vibrational  quanta  can  be  identified.  Least  squares  spectral  fitting  analysis  determines  the 
vibrational  state  column  densities  for  both  the  ambient  and  beam-excited  components,  and  their 
dependencies  on  altitude  and  dosing.  Systematic  kinetic  analysis  of  these  results  yields  new 
information  on  the  processes  responsible  for  electron-initiated  ozone  excitation.  Direct 
excitation  of  03(v)  by  secondary  electron  impact  on  ground  state  ozone  clearly  plays  a  key 
role.  In  addition,  the  observations  at  the  highest  dosing  levels  (- 102  km)  indicate  the  formation 
of  additional  ozone  in  the  beam  through  bimolecular  reactions  involving  beam-excited 
metastable  precursors.  The  observations  match  well  with  the  expected  kinetics  of  the  reaction 
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of  beam-created  02(A3\*,v)  with  02  to  produce  03,  which  is  either  formed  in  vibrationally 
excited  states  or  subsequently  excited  by  electron  impact. 

1.  Introduction 

The  EXCEDE  HI  upper  atmospheric  rocket  experiment  observed  spectrally  resolved 
atmospheric  infrared  fluorescence  excited  in  the  lower  thermosphere  (85  to  1 15  km)  by  a  pulsed, 
2.5  kV,  18  A  electron  beam  generated  by  a  rocketbome  accelerator  module.  A  separate, 
extensively  instrumented  sensor  module  included  a  cryogenic  scanning  Michelson  interferometer 
which  produced  Fourier  transform  spectra  of  the  ambient  and  beam-created  fluorescence  over  the 
2.4  to  21  fi m  (4200  to  475  cm'1)  spectral  region  with  ~2  cm'1  spectral  resolution.  The  spectra 
reveal  a  rich  assortment  of  beam-created  excited  state  species,  including  NO,  NO+,  COz,  and  03. 
This  paper  describes  the  observations  and  analysis  of  beam-created  fluorescence  from  the 
asymmetric  stretching  fundamental  (v3)  of  vibrationally  excited  03(v)  near  10  fim,  observed  at 
altitudes  near  100  km. 

The  kinetics  and  mechanisms  of  non-equilibrium  03(v)  excitation  in  the  quiescent  upper 
mesosphere  and  lower  thermosphere  are  fairly  well  understood  [Rawlins,  1985;  Solomon  et  al, 
1986].  However,  the  question  of  auroral  excitation  or  enhancement  of  ozone  is  still  unresolved. 
This  issue  is  also  relevant  to  mechanisms  for  electrical  and  photochemical  production  of  ozone  in 
the  middle  and  lower  atmosphere.  In  the  quiescent  upper  atmosphere,  O3(001)  is  populated  by 
thermalizing  collisions  and  by  resonant  absorption  of  upwelling  10  fim  Earth  and  atmospheric 
radiation,  and  03(v3s  1)  is  generated  by  three-body  recombination  of  O  and  02.  Additional 

2 

112 


processes  for  excitation  of  03(v)  in  electron-bombarded  air  have  been  suggested,  but  the  avail¬ 
able  auroral  data  base  has  not  clearly  demonstrated  their  occurrence.  Potential  03(v)  formation/ 
excitation  mechanisms  include  direct  electron  impact  excitation  of  03,  energy  transfer  to  03  from 
metastable  electronically  or  vibrationally  excited  species,  and  chemiexcitation  reactions  of  the 
type  02*  +  02  -  03(v)  +  O;  these  mechanisms  were  discussed  by  Rawlins  [1985]. 

Because  of  the  low  precursor  concentrations  for  these  processes,  their  effects  are  likely  to 
be  observable  above  ambient  only  for  relatively  strong  dosings  at  low  altitudes,  e.g.,  an  IBC  HI  or 
nr  aurora  near  90  to  100  km.  These  conditions  are  rarely  encountered  by  auroral  rocketbome 
flight  measurements,  and  the  limited  existing  data  are  qualitative  and  inconclusive.  The  HIRIS 
rocketbome  interferometer  spectrometer  observed  a  very  intense  IBC  III  to  IV  arc  [Stair  et  al, 
1983;  Rawlins  et  al,  1985],  with  an  apparently  large  03(v)  excitation  efficiency  of  2  to  5  vibra¬ 
tional  quanta  per  ion  pair  near  100  km.  Unfortunately,  this  cannot  be  taken  to  be  a  definitive 
result,  since  tumbling  motion  of  the  payload  introduced  a  large  uncertainty  in  the  altitude  of  the 
amorally  excited  03(v)  emission.  A  more  recent  measurement  by  the  SPIRIT  1  rocketbome 
interferometer  [Rawlins  et  al,  1993]  observed  03(v)  emission  near  100  km  in  the  vicinity  of  an 
IBC  EQ  arc,  however  the  instrument  did  not  get  a  sufficiently  direct  view  of  the  arc  to  con¬ 
clusively  distinguish  the  observed  03(v)  emission  as  auroral  rather  than  intervening  ambient 
fluorescence.  These  observations  resulted  in  an  upper  bound  of  0.8  to  2.3  quanta  per  ion  pair, 
which  is  consistent  with  HIRIS  but  nevertheless  inconclusive.  In  order  to  clearly  assess  the 
possible  excitation  kinetics  and  scalings  with  dosing  conditions,  it  is  essential  to  obtain 
quantitative  observations  of  03(v)  emission  from  a  dosed  volume  at  a  well  determined  altitude. 
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A  recent  series  of  laboratory  electron-beam  dosing  experiments  provides  considerable 
insight  into  the  excitation  of  03(v)  fluorescence  in  the  bombardment  of  atmospheric  species  by 
energetic  electrons  [Upschulte  et  al.,  1989;  Holtzclaw  et  al,  1992].  These  experiments, 
performed  in  the  cryogenic  LABCEDE  facility  [Upschulte  et  al,  1994a]  at  the  Phillips 
Laboratory/Geophysics  Directorate,  observed  infrared  03(v)  fluorescence  from  electron- 
irradiated  Oj/Ar,  Oj/Nj,  Oj/Ar,  and  0^/N2  mixtures  at  pressures  of  0.1  to  80  mtorr  and 
temperatures  near  90  K.  The  measurements  demonstrated  the  following: 

(1)  Excitation  of  03(v)  by  secondary  electron  impact  on  03  is  highly  efficient,  with  a 
cross  section  approximately  10  times  larger  than  that  of  the  analogous  excitation 
of  C02(v3); 

(2)  03  is  formed  from  electron  irradiation  of  02  at  low  pressures,  apparently  by  a 
two-body  reaction  involving  beam-created  metastable  oxygen  precursors. 

The  results  from  the  e'  +  03  excitation  measurements  indicate  an  auroral  excitation  rate  which 
scales  with  secondary  electron  number  density  rather  than  with  dosing  strength  (ion  pair 
production  rate),  and  which  should  be  readily  observable  in  zenith  or  limb  views  of  BBC  HI  or 
greater  auroras  at  » 100  km  [Blumberg  et  al,  1990]. 

Recent  laboratory  kinetics  investigations  reported  by  Copeland  and  coworkers  [Copeland, 
1994;  Copeland  et  al,  1994;  Knutsen  et  al,  1994]  have  determined  room-temperature  rate 
coefficients  for  the  quenching  of  the  metastable  02(A%*,  6zvz9)  by  02.  These  investigators 
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find  the  quenching  rate  coefficients  to  be  orders  of  magnitude  faster  than  previously  believed 
[Kenner  and  Ogryzlo,  1980],  and  observe  substantial  production  of  atomic  oxygen,  suggesting  a 
high  probability  for  reactive  formation  of  ozone.  Thus  both  the  electron  impact  excitation  and 
the  two-body  reactive  production  of  ozone  now  appear  to  be  considerably  more  efficient  than 
previously  suggested  by  Rawlins  [1985],  and  one  can  expect  that  these  reactions  should  con¬ 
tribute  to  03(v)  excitation  in  the  electron-irradiated  upper  atmosphere. 

In  the  present  paper,  we  report  observations  from  a  high-altitude,  rocketbome  electron 
irradiation  experiment  in  which  spectrally  resolved  infrared  fluorescence  was  observed  from  the 
bombardment  of  the  ambient  atmosphere  by  a  pulsed  electron  beam  at  altitudes  of  85  to  1 15  km. 
This  bombardment  produced  a  maximum  dosing  level  of  6  x  10n  ion  pairs/(cm3s)  near  103  km 
and  electron  number  densities  as  high  as  1.5  x  109  cm'3,  with  gas  residence  times  of  20  to  600  ms 
in  the  irradiated  volume.  These  measurements  provide  the  first  conclusive  and  quantitative 
atmospheric  observations  of  03(v)  excitation  by  electron-initiated  processes.  The  results  are 
consistent  with  the  recent  laboratory  observations  of  efficient  e'  +  03  excitation  and  rapid 
reaction  of  beam-created  02*  metastables  with  02  to  form  03. 

2.  Summary  of  Flight  Measurements 

The  instrumentation  and  flight  of  EXCEDE  HI  are  described  elsewhere.  The  payload  was 
launched  aboard  an  Aries  booster  on  27  April,  1990,  at  White  Sands,  NM,  01:00:59  MDT 
(07:00:59  GMT),  and  achieved  an  apogee  of  1 15  km.  The  booster  was  separated  from  the 
payload  at  68  km.  At  82.7  km  on  ascent,  the  payload  separated  into  two  modules,  one  containing 
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the  electron  beam  accelerator  system  and  supporting  instrumentation,  and  one  containing  a  suite 
of  optical  spectrometers  and  photometers  for  viewing  the  beam-excited  volume.  The  accelerator 
system  generated  pulses  of  nominally  2.6  kV  and  18  A  in  a  series  of  7.1-second  cycles  (4.73  s  on, 
2.37  s  off).  Eighteen  usable  electron  beam  cycles  occurred  between  -95  km  on  ascent,  115  km 
apogee,  and  87  km  on  descent. 

For  the  present  analysis,  the  primary  instrument  aboard  the  sensor  module  was  a  cryo¬ 
genic  infrared  interferometric  spectrometer  (2.4  to  21  pm),  which  was  used  to  obtain  high- 
resolution  (=2  cm1)  fluorescence  spectra  of  the  excited  03(v3)  band  between  9  and  12  pn l. 
Additional  sensors  included  three  circular- variable-filter  spectrometers  (2.4  to  21  pm),  a  visible 
grating  monochromator  (300  to  800  nm),  an  ultraviolet  grating  spectrometer  (180  to  300  nm), 
two  X-ray  proportional  counters  [Rappaport  et  al,  1993],  four  spatially  scanning  filter 
photometers  centered  at  391.4, 380.4, 557.7,  and  276.1  nm,  additional  fixed  filter  photometers  at 
391.4, 520.0,  and  522.8  nm,  an  atomic  oxygen  resonance  fluorescence  light  source  and  detection 
system,  and  video  and  film  cameras.  For  the  ozone  excitation  analysis  described  here,  the  key 
supporting  sensors  were  the  scanning  391.4  nm  photometer,  which  determined  the  local  ion  pair 
production  rate  in  the  interferometer's  field  of  view  for  each  beam  pulse,  and  the  vacuum- 
ultraviolet  (130  nm)  resonance  fluorescence  diagnostic,  which  provided  estimates  of  the  ambient 
atomic  oxygen  number  densities  [Nelson  et  al,  1992]. 

The  interferometer  was  a  cryogenic,  telescoped,  scanning  Michelson  interferometer  with 
dual  optical  channels  [Thurgood  and  Huppi,  1989].  The  instrument  was  configured  to  operate 
simultaneously  in  four  spectral  bandpasses:  2.4  to  4.0  pm,  4.0  to  7.5  pm,  7.0  to  12.5  pm,  and 
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12.5  to  21  fim.  Careful  selection  of  the  cut-on  and  cut-off  wavelengths  of  the  bandpass  filters 
allowed  sensitive  observations  of  the  principal  high-altitude  radiators  (NO,  03,  and  COj)  without 
interferences  from  detector  non-linearity  or  dynamic  range  limitations.  The  interferometer's 
scanning  cycle  and  the  electron  beam  cycle  were  synchronized  such  that,  in  each  beam  cycle, 
there  were  two  complete  interferometer  scans  with  the  beam  on  and  a  third  complete  scan  and 
rapid  retrace  with  the  beam  off.  The  instrument  produced  single-sided  interferograms,  which 
were  detrended,  phase-corrected,  and  transformed  using  triangular  apodization.  Typical  spectral 
resolutions  of  the  Fourier  transform  spectra  were  1.7  to  1.9  cm'1,  depending  on  the  number  of 
data  points  in  the  interferogram. 

The  field  of  view  of  the  interferometer  was  3.6  deg  x  3.6  deg  along  an  upward-looking 
line  of  sight  as  illustrated  in  Figure  1.  The  beam  elevation  was  61.1  deg,  constrained  by  the  local 
magnetic  field  lines.  Throughout  most  of  the  measurements,  the  included  angle  between  the 
beam  and  the  line  of  sight  was  =47  deg;  however  this  angle  increased  late  in  the  descent  as  the 
sensor  scanned  back  toward  the  accelerator  module.  The  spatially  scanning  391.4  nm 
photometer  provided  determinations  of  the  energy  deposition  within  the  observed  beam  volume 
throughout  the  flight.  The  volumetric  ion  pair  production  rates  were  determined  from  the 
photometric  data  assuming  14.1  N2+  ions  per  391.4  nm  photon  [Borst  and  Zipf,  1970]  and  0.62 
N2+  ion  per  total  ion  [Rees  and  Jones,  1973],  and  using  the  observed  beam  widths  from 
generalized  Gaussian  fits  to  the  spatial  distributions.  Typical  beam  diameters  were  =8  m.  The 
local  volumetric  ion  pair  production  rates  are  plotted  for  ascent  and  descent  in  Figure  2. 

Figure  1  and  Figure  2 
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During  the  early  stage  of  the  ascent,  the  velocity  vector  of  the  accelerator  module  was 
directed  along  the  direction  of  beam  propagation,  resulting  in  longer  gas  residence  times  in  the 
irradiated  volume  for  altitudes  up  to  about  110  km.  As  the  module  approached  apogee  and 
subsequently  descended,  the  flow  direction  became  more  nearly  orthogonal  to  the  beam,  and  the 
residence  times  were  correspondingly  shorter.  Flow  times  through  the  viewed  dosed  volume 
were  determined  from  the  velocities  and  beam  dimensions,  and  are  plotted  in  Figure  3.  The 
characteristic  flow  times  on  ascent  up  to  105  km  were  well  in  excess  of  100  ms,  in  contrast  to 
values  of  =30  ms  or  less  near  apogee  and  throughout  the  descent.  The  long  residence  times  on 
ascent  resulted  in  the  largest  energy  depositions  of  the  flight,  in  excess  of  1012  eV/cm3  below 
105  km  and  peaking  at  =  1013  eV/cm3  near  101  km.  The  difference  in  residence  times  between 
the  maximum  dosing  on  ascent  below  105  km  and  the  remainder  of  the  flight  is  a  key  element  in 
the  kinetic  interpretation  of  the  observed  excitation  phenomena.  As  discussed  further  below,  the 
residence  times  near  maximum  dose  are  long  enough  for  the  beam  excitation  processes  to  closely 
approach  steady  state  kinetics  conditions,  while  the  shorter  times  near  apogee  and  throughout 
descent  more  closely  sample  the  initial  excitation  rates. 

Figure  3 

3.  03  Spectral  Data  and  Analysis 

The  spectral  analysis  centered  on  the  asymmetric  stretch  (v3)  band  of  03  near  1000  cm'1 
(10  fim).  Example  spectra  of  03(v)  fluorescence,  observed  on  ascent  with  the  beam  off,  are 
shown  in  Figure  4a.  At  altitudes  below  =  105  km,  the  (001)  -  (000)  band,  centered  at  1042  cm'1, 
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is  clearly  identifiable  by  its  characteristic  R,  Q,  and  P  branch  rotational  structure.  The  vibrational 
"hot  bands”,  resulting  from  Av3=l  transitions  from  upper-state  vibrational  levels  at  higher  energy 
than  (001),  are  discernible  at  transition  frequencies  below  about  1000  cm'1.  Fluorescence  from 
these  states  is  detectable  at  altitudes  below  =97  km,  and  is  below  the  instrument  noise-equivalent 
signal  radiance  (NESR)  above  that  altitude.  The  fluorescence  from  the  (001)  and  higher  states 
follows  the  spectral  and  altitude  behavior  expected  for  ambient  nighttime  atmospheric  03 
fluorescence  [Green  et  al,  1986;  Rawlins  et  al,  1993].  The  spectra  also  show  a  weak,  broad 
feature  centered  near  1070  cm'1.  The  intensity  of  this  feature  is  essentially  independent  of 
altitude  up  to  115  km,  a  behavior  signifying  a  rocketbome  contaminant  species.  Based  on  prior 
experience  with  rocketbome  infrared  spectrometer  measurements,  this  feature  is  probably  due  to 
earthshine-pumped,  near-Held  fluorescence  from  outgassing  products  of  silicon  oxide  sealants 
used  on  the  payload  [Smith,  1988].  The  presence  of  this  feature  does  not  significantly  affect  the 
analysis  of  the  03(v)  fluorescence.  As  shown  by  the  107  km  spectrum  in  Figure  4a,  the  spectral 
scans  for  altitudes  above  105  km  show  a  constant  baseline  radiance  of  =  10'12  W/(cm2  sr  cm'1) 
between  900  and  1040  cm'1,  signifying  the  effective  detection  limit  of  the  instrument  in  this 
spectral  region. 

Figure  4 

Fluorescence  spectra  observed  during  ascent  with  the  electron  beam  on  are  shown  in 
Figure  4b,  illustrating  the  behavior  above  and  below  the  altitude  threshold  for  detectable  03(v) 
emission.  The  107  km  "baseline"  spectrum  exhibits  the  same  noise  floor  between  900  and 
1040  cm'1  as  with  the  beam  off,  with  the  addition  of  a  few  sharp  atomic  lines  due  to  0(1)  Rydberg 
transitions.  The  contamination  feature  centered  near  1070  cm'1  remains  unchanged,  however  the 
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intensity  above  1 100  cm'1  is  enhanced  by  the  election  beam.  This  fluorescence  has  been 
identified  as  high-J'  transitions  of  rovibrationally  excited  NO(vJ)  formed  by  beam-initiated 
chemical  processes  [Fraser  et  al.,  1993;  Upson  et  alM  1993;  Armstrong  et  al.,  1995].  The 
majority  of  the  fluorescence  from  this  band  system  lies  at  higher  frequencies  (roughly  1600  to 
2000  cm'1),  and  will  not  be  discussed  further  here.  The  NO(v,J)  fluorescence  does  not  appear  to 
be  significant  below  « 1 100  cm'1,  and  thus  does  not  impact  the  analysis  of  the  03(v)  system.  As 
shown  by  the  102.5  km  spectrum  in  Fig.  4b,  the  03(v)  fluorescence  spectra  with  the  beam  on 
appear  qualitatively  similar  to  those  observed  from  the  ambient  atmosphere  with  the  beam  off. 
However,  the  beam-on  spectra  show  significantly  more  fluorescence  from  the  higher  ozone 
vibrational  levels  (900  to  1000  cm'1)  than  the  ambient  spectra  at  similar  altitudes,  as  illustrated  by 
comparison  of  the  102.5  km  spectrum  in  Fig.  4b  with  the  101  km  spectrum  in  Figure  4a. 

The  spectral  data  were  subjected  to  more  detailed  analysis  in  order  to  determine 
systematic  trends  in  03(v)  formation  with  altitude  and  dosing.  After  subtracting  the  instrumental 
baseline,  the  spectra  were  analyzed  by  a  least  squares  spectral  fitting  method  developed  for  the 
analysis  of  03(v)  spectra  by  Rawlins  and  Armstrong  (1987)  and  improved  by  Rawlins  et  al. 
(1993).  The  03(v)  spectrum  was  simulated  using  two  models,  one  incorporating  nine  vibrational 
bands  and  one  incorporating  eighteen.  The  9-band  model  includes  transitions  from  the  (00v3) 
and  (10(v3-l))  vibrational  states,  1  £  v3  s  5,  with  the  relative  populations  of  each 
(00v3y(10(v3-l))  dyad  set  to  local  thermodynamic  equilibrium  values.  This  accounts  for  the 
likely  excitation  of  the  Vj  mode  through  near-resonant  vibrational  energy  transfer,  and  reduces  the 
complexity  due  to  spectral  overlap  between  neighboring  bands.  The  least  squares  fitting 
procedure  then  gives  the  total  column  density  for  each  dyad  required  to  minimize  the  variance 
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between  the  computed  and  observed  spectra.  The  18-band  model  is  similar  and  includes 
additional  transitions  from  (01v3y(ll(v3-l))  dyads,  thus  including  the  possibility  of  excitation  in 
the  v2  (bending)  mode.  Typical  spectral  fitting  results  are  illustrated  in  Figure  5  for  the  18-band 
model. 

Figure  5 

In  general,  the  two  spectral  models  give  veiy  similar  results.  The  beam-off  spectra  tend 
to  be  more  consistent  with  the  9-band  model,  and  the  beam-on  spectra  are  equally  consistent  with 
both  models.  However,  it  is  not  possible  to  distinguish  clearly  whether  any  v2  combination  states 
other  than  (01 1)  are  significantly  populated  in  the  ambient  atmosphere  or  in  the  beam-irradiated 
volume.  Nevertheless,  as  in  previous  analyses  of  laboratory  [Rawlins  and  Armstrong,  1987]  and 
atmospheric  [Rawlins  et  al,  1993]  03(v)  spectra,  it  is  clear  from  the  EXCEDE  El  spectra  that 
population  of  at  least  the  (00v3)  and  (10(v3-l))  states  is  required  to  account  for  the  nearly  con¬ 
tinuous,  unstructured  spectral  distributions  in  the  "hot  band"  region.  When  the  populations  are 
lumped  into  a  common  upper  state  v  =  v,  +  v2  +  v3,  the  values  for  the  beam-off  condition  agree 
well  with  the  relative  population  distributions  predicted  by  the  steady-state  kinetic  model  of 
Rawlins  (1985)  for  2  £  v  <;  5.  Similar  agreement  was  found  for  SPIRIT  1  03(v)  spectra 
[Rawlins  et  al,  1993].  In  contrast,  the  lumped  populations  for  the  beam-on  case  indicate  greater 
excitation  at  the  higher  levels,  as  discussed  further  below. 
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4. 


03  (v)  Altitude  Profiles 


In  order  to  quantify  the  03(v)  component  excited  in  the  electron  beam,  it  was  first 
necessary  to  determine  the  contribution  from  the  ambient,  quiescent  atmosphere.  The  beam-off 
03(v)  column  densities  for  ascent  and  descent  were  fit  to  parametric  expressions  for  altitudes 
from  84  to  105  km.  These  expressions  were  then  used  to  correct  the  beam-on  quantities  for  the 
ambient  component.  The  residuals  from  these  corrections  give  a  measure  of  the  amount  of  03(v) 
produced  by  beam-initiated  processes. 

The  observed  column  densities  for  ambient  03(v)  on  ascent  and  descent  are  plotted  in 
Figure  6.  The  03(v=l)  column  densities  are  fairly  well  (±10%  or  better)  determined,  especially 
above  92  km.  However,  only  four  beam-off  spectra  had  sufficient  signal/noise  ratio  to  permit 
determination  of  the  03(v^2)  column  densities,  and  values  above  97  km  were  below  the 
detection  limit.  The  solid  curves  in  Figure  6  represent  polynomial  fits  to  the  data  points.  The 
03(v=l)  data  were  fit  to  a  cubic  polynomial  in  log(column  density)  vs  altitude,  and  the  03(v^2) 
data  were  fit  to  quadratic  polynomials.  The  03(v=5)  column  densities  are  highly  uncertain, 
especially  above  95  km;  however  the  03(v=2-4)  values  are  well  determined  and  exhibit  an 
altitude  dependence  characteristic  of  quiescent  nighttime  03(v). 

Figure  6 

The  03(v=1-3)  column  densities  observed  on  ascent  and  descent  with  the  beam  on  are 
shown  in  Figure  7,  overlaid  with  the  curvefits  to  the  ambient  values  from  Figure  6.  This 
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comparison  illustrates  the  departure  from  the  ambient  vibrational  distributions  and  altitude 
dependencies  observed  with  the  electron  beam  on,  especially  during  ascent  between  95  and 
103  km  when  the  dosing  residence  times  are  longest.  In  particular,  the  beam-induced  enhance¬ 
ment  in  03(v^2)  above  98  km  is  readily  apparent.  Figure  8  shows  the  beam-on  column  densities 
for  v=l,3  plotted  as  ratios  to  the  corresponding  beam-off  values  interpolated  from  the  curvefits. 
This  figure  shows  that  the  apparent  enhancement  above  95  km  is  as  large  as  a  factor  of  2  for  v=l 
and  a  factor  of  =9  for  v=3.  Since  the  optical  path  through  the  electron  beam  is  much  shorter  than 
the  viewed  ambient  atmospheric  column,  these  observations  indicate  dramatic  enhancement  in 
the  03(v)  number  densities  within  the  beam-irradiated  volume. 

Figure  7  and  Figure  8 

The  path-averaged  03(v)  number  densities  excited  by  the  electron  beam  between  98  and 
105  km  are  shown  in  Figure  9.  These  values  were  determined  by  subtracting  the  curvefit 
ambient  column  density  at  each  altitude  from  the  observed  beam-on  column  density,  and  dividing 
the  result  by  the  viewed  path  length  across  the  beam.  Total  03(v)  number  densities  of  ~2  x 
107  cm'3  were  generated  by  the  electron  beam  in  the  maximum-dose  configuration  on  ascent,  and 
=3  x  106  cm'3  were  produced  near  104  km  for  short  dosing  time  on  descent.  As  can  be  seen  from 
Figs.  7  and  8,  there  is  no  clear  indication  for  any  beam-enhanced  03(v)  on  descent  for  altitudes 
between  85  and  100  km,  i.e.  the  column  densities  observed  with  the  beam  on  are  within  ±20%  of 
the  curvefit  ambient  values.  The  final  descent  spectrum,  at  84  km,  appears  to  indicate  observable 
03(v)  enhancement  above  the  ambient  level,  however  the  field  of  view  through  the  beam-dosed 
region  is  large  and  highly  uncertain,  precluding  quantitative  analysis  for  that  case. 
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Figure  9 


5.  Excitation  Kinetics  and  Mechanisms 

The  beam-produced  03(v)  number  densities  shown  in  Figure  9  are  extremely  large. 
Based  on  MSIS90  [Hedm,  1991]  atmospheric  profiles  and  the  03(v)  kinetics  discussed  by 
Rawlins  et  al.  (1993),  the  ambient  03(v)  number  densities  at  103  km  should  be  =2  x  104  cm'3  for 
v=l  and  ^80  cm'3  for  v^2.  Thus  03(v=l)  in  the  dosed  volume  is  enhanced  by  up  to  three  orders 
of  magnitude,  and  03(v^2)  is  enhanced  by  more  than  four  orders  of  magnitude.  Indeed,  the  total 
03(v)  number  density  produced  by  the  beam  in  the  maximum-dose  configuration  at  103  km  is 
comparable  to  the  total  03  number  density  expected  for  that  altitude  for  normal  quiescent 
conditions.  Clearly,  this  cannot  occur  unless  chemical  processes  initiated  by  the  electron  beam 
are  generating  new  ozone  at  a  substantial  rate.  It  is  also  obvious  that  three-body  recombination 
of  O  and  02  cannot  account  for  this  observation.  We  present  here  a  simple  kinetic  analysis  to 
identify  the  03(v)  formation  mechanisms,  focussed  on  the  data  for  103  km  on  ascent  and  104  km 
on  descent. 

We  summarize  our  analysis  first  and  then  present  a  more  detailed  discussion.  To  first 
order,  the  primary  deactivation  process  for  03(v=l)  at  103-104  km  is  radiative  decay  at  =11  s'1. 
As  discussed  further  below,  superelastic  deactivation  by  electrons  contributes  another  =5  s'1.  At 
104  km  on  descent,  the  gas  residence  time  in  the  beam  is  only  =20  ms  (cf.  Figure  3),  which  is 
substantially  less  than  the  kinetic  lifetime  of  03(v=l).  Under  these  conditions,  the  chemical 
losses  during  the  residence  time  are  relatively  small,  and  the  data  provide  a  measure  of  the  initial 
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production  rate  of  03(v=l).  This  rate  is  then  ~2  x  108  cm'V1,  and  is  likely  due  to  vibrational 
excitation  by  secondary  electron  impact  on  ambient  ozone.  In  contrast,  for  103  km  on  ascent,  the 
gas  residence  time  substantially  exceeds  the  chemical  lifetime,  so  the  system  is  approaching 
steady  state.  From  the  steady  state  approximation,  the  apparent  03(v)  production  rate  at  103  km 
on  ascent  is  then  =3  x  108  cm'V1,  comparable  to  the  value  estimated  from  the  initial  rate. 
However,  excitation  of  ozone  at  this  rate  will  deplete  O3(v=0)  on  the  100  ms  time  scale,  such  that 
the  actual  excitation  rate  would  have  to  be  significantly  larger  to  achieve  the  observed  03(v=l) 
number  densities.  Indeed,  most  or  all  of  the  ambient  03  in  the  beam  would  have  to  be 
vibrationally  excited,  a  condition  that  cannot  be  achieved  because  of  radiative  and  superelastic 
deactivation  from  v=l  back  to  v=0.  Thus  either  the  ambient  ozone  number  density  is  much 
larger  than  inferred  from  the  ambient  radiance  data,  or  a  beam-initiated  chemical  process  makes 
more  ozone  in  the  dosed  volume.  Based  on  obvious  kinetic  considerations,  such  a  beam-initiated 
process  would  likely  be  a  two-body  reaction  involving  at  least  one  beam-created  metastable 
precursor,  e.g.  02*  +  02. 

Electron  Impact  Excitation 

We  have  modeled  the  kinetics  of  the  electron  beam  excitation  process  using  model- 
derived  species  profiles,  secondary  electron  distributions,  and  the  reaction  mechanisms  listed  in 
Table  1.  Concentration  profiles  for  03  and  O  were  derived  from  the  ambient  03(v)  data  using 
the  density  and  temperature  profiles  given  by  MSIS90  together  with  the  kinetic  analysis 
described  by  Rawlins  et  al.  (1993).  For  this  purpose,  the  ambient  03(v)  curvefit  column  densities 
were  corrected  to  zenith  values  and  were  differentiated  with  respect  to  altitude  to  give  local 
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number  density  profiles.  03(v=l)  number  densities  were  used  to  compute  O3(v=0)  number 
densities  via  radiative/collisional  steady  state,  and  03(vk2)  number  densities  were  used  to 
estimate  three-body  recombination  rates  and  hence  O  number  densities.  The  inferred  O  and  03 
number  density  profiles  are  compared  to  those  obtained  from  MSIS90  in  Figure  10.  Here  the 
"MSIS"  03  number  densities  were  computed  from  the  MSIS  O,  H,  T,  and  density  profiles 
assuming  steady-state  formation  of  03  by  three-body  recombination  and  loss  of  03  by  reaction 
with  H  and  O.  The  atomic  oxygen  profile  derived  from  03(v)  is  intermediate  between  that  of 
MSIS  and  that  determined  by  resonance  fluorescence  measurements  on  board  the  sensor  module 
[Nelson  et  al,  1992].  The  resonance  fluorescence  measurements  indicate  a  peak  O  number 
density  near  10"  cm'3  between  100  and  105  km,  and  =2.5  x  10l°  cm'3  at  85  km,  with  factor-of- 
two  estimated  uncertainty.  The  derived  03  profile  is  significantly  different  from  the  MSIS 
expectation,  and  suggests  an  elevated  density  profile  due  to  a  dynamic  perturbation.  In 
particular,  the  apparent  03  number  density  near  103  km  is  roughly  twice  that  of  the  MSIS 
estimate. 

Table  1  and  Figure  10 

Estimates  of  the  electron  number  densities  at  apogee  and  maximum  dose  and  a  computed 
secondary  electron  energy  distribution  for  apogee  were  provided  by  F.  Bien  and  J.  Duff  (private 
communication,  1993),  based  on  detailed  modeling  of  the  degradation  of  the  primary  electron 
energy  and  flux  through  inelastic  collisions  in  the  atmosphere.  The  modeled  secondary  electron 
energy  distribution  for  apogee  (115  km)  was  non-Boltzmann,  with  a  peak  between  0.05  and 
0.5  eV  and  a  tail  extending  to  10'7  of  the  peak  at  100  eV.  The  low-energy  portion  of  this 
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distribution  is  similar  to  that  for  an  electron  "temperature"  of  =0.5  eV.  We  expect  the  secondary 
electron  "temperatures"  at  lower  altitudes  to  be  "cooler"  due  to  increased  collision  rates. 

Using  these  results  and  the  approximation  of  steady  state  between  ion  pair  production  and 
electron-ion  recombination,  we  estimated  electron  densities  for  the  remaining  beam-on  cases. 

The  estimated  values  ranged  from  4  x  108  cm'3  at  apogee  (115  km)  to  1.5  x  109  cm'3  at  maximum 
dose  (103  km  ascent).  However,  we  note  that,  for  these  values,  the  e-folding  time  for  electron- 
ion  recombination  is  about  10  to  20  ms,  which  is  a  substantial  fraction  of  the  beam  residence 
times  for  the  apogee  and  descent  cases  (see  Figure  3).  Hence  the  electron-ion  steady  state 
approximation  is  less  accurate  for  these  cases,  and  provides  an  overestimate  of  the  mean  electron 
number  density  during  the  beam  pulse.  However,  for  present  purposes  of  rough  kinetic 
estimates,  this  approximation  is  adequate,  introducing  negligible  error  at  maximum  dosing 
(103  km  ascent)  and  less  than  30%  error  for  104  km  descent. 

To  model  the  electron  impact  excitation  process,  we  require  the  first  30  reactions  listed 
in  Table  1,  which  describe  the  excitation,  deactivation,  and  dissociation  of  03  by  secondary 
electrons,  and  vibrational  deactivation  of  03(v)  by  radiative  and  collisional  processes.  To 
describe  the  electron  impact  excitation  (Reactions  1  to  5),  the  total  rate  coefficient  was  varied  to 
match  the  observed  total  vibrational  excitation,  and  the  partitioning  by  vibrational  level  was 
taken  from  the  laboratory  measurements  [Holtzclaw  et  al.,  1991, 1992].  The  rate  coefficients  for 
superelastic  deactivation  (Reactions  6  to  9)  were  then  obtained  from  microscopic  reversibility 
using  0.5  eV  as  an  approximation  for  the  secondary  electron  "temperature."  Rate  coefficients  for 
the  dissociation  of  03(v)  (Reactions  10  to  15)  were  estimated  by  convolving  the  modeled 
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secondary  electron  energy  distribution  with  the  cross  section  of  Keto  (1981)  for  dissociation  of 
O3(v=0),  where  the  values  for  03(v^  1)  were  determined  by  adjusting  the  cross  section's  threshold 
by  an  amount  corresponding  to  the  vibrational  energy  of  the  molecule. 

The  radiative  rates  and  collisional  deactivation  rate  coefficients  were  taken  from 
laboratory  determinations  as  cited.  The  collisional  deactivation  processes  were  included  for 
completeness,  however  their  rates  are  too  slow  to  contribute  significantly  to  03(v)  losses  on  the 
time  scale  of  the  measurements.  The  total  first-order  loss  rates  in  the  beam  for  03(v=l)  and 
O3(v=0)  between  95  and  105  km  are  estimated  to  be  =20  and  3-8  s'1  respectively.  These  values 
are  too  slow  compared  to  the  beam  processing  rates  to  permit  the  steady-state  approximation  on 
these  species  except  near  maximum  dosing,  so  it  was  necessary  to  numerically  integrate  the 
coupled  rate  equations. 

The  results  of  the  model  calculations  were  compared  to  the  observations  of  beam-excited 
03(v)  at  maximum  dosing  (103  km  ascent)  and  descent  (104  km),  as  well  as  to  the  absence  of 
observable  03(v)  excitation  on  descent  below  100  km.  The  total  03(v)  number  density  observed 
at  maximum  dosing  is  consistent  with  a  total  rate  coefficient  of  =5  x  10*8  cm3/s  for  electron 
impact  excitation  of  03(v);  the  excitation  is  accompanied  by  dissociation  of  about  2/3  of  the 
ambient  03  in  the  beam  volume.  However,  this  model  predicts  essentially  no  difference  between 
the  beam-enhanced  03(v)  number  densities  on  ascent  and  descent,  with  a  predicted  enhancement 
at  100  km  of  a  factor  of  4  in  excess  of  the  detection  limit.  In  order  to  accommodate  the 
excitation  observed  at  104  km  (descent)  and  the  absence  of  observable  excitation  on  descent 
below  that  altitude,  it  is  necessary  to  use  a  total  electron  impact  excitation  rate  coefficient  an 
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order  of  magnitude  smaller,  =5  x  10'9  cm3/s.  This  then  results  in  a  substantial  underprediction  of 
the  excitation  observed  at  maximum  dosing,  suggesting  that  another  reaction  must  be  generating 
03  in  the  beam  volume.  Assuming  that  we  have  accounted  for  all  the  loss  processes,  the 
additional  03  production  rate  required  to  account  for  the  excitation  at  103  km  (ascent)  is 
~3  x  108  cm'V. 

Metastable  Oxygen  Chemistry 

The  Oj  production  rate  deduced  from  the  maximum  dose  data  is  consistent  with  a 
moderately  efficient  (*10%  collision  efficiency)  bimolecular  reaction  between  a  major  ambient 
atmospheric  species  and  a  beam-created  trace  species,  most  likely  metastable.  Bimolecular 
reactions  which  could  potentially  contribute  to  03  formation  and/or  03(v)  excitation  in  electron- 
bombarded  oxygen  or  air  have  been  reviewed  in  previous  papers  [Rawlins  et  al,  1985;  Rawlins, 
1985;  Rawlins  et  al,  1987].  Of  these  reactions,  the  one  most  likely  to  be  responsible  for  the 
EXCEDE  HI  observations  at  high  altitude  is  the  process  of  the  form 

02(A)  +  02  -  Oj  +  O  .  (Rl) 

Here  02(A)  is  a  generic  representation  of  any  or  all  of  the  three  nested  metastable  states  A3X, 
A,3II,  and  c‘l,  which  lie  at  energies  above  4  eV.  Due  to  the  high  energy  of  02(A),  this  process  is 
exoergic,  and  has  been  previously  suggested  for  both  the  auroral  [Rawlins  et  al,  1985;  Rawlins, 
1985]  and  quiescent  [Allen,  1986]  thermosphere.  Since  the  reverse  reaction  provides  adiabatic, 
statistically  favored  pathways  for  the  well-known  reaction  of  O  with  03,  it  is  reasonable  to  expect 
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that  the  forward  reaction  can  occur.  However,  very  little  is  known  about  the  kinetics  of  02(A)  or 
its  production  in  electron-disturbed  air.  In  addition,  previous  measurements  [Kenner  and 
Ogryzlo,  1980]  of  the  rate  coefficients  for  quenching  of  02(A3S)  by  02  indicated  that  the  reaction 
was  too  slow  to  account  for  significant  03  production  in  the  upper  atmosphere,  given  the 
likelihood  for  relatively  low  02(A)  number  densities  even  in  auroral  conditions.  However, 
recently  reported  laboratory  measurements  now  show  that  the  reactions  of  02(A3£,v)  with  02  are 
indeed  highly  efficient,  and  are  likely  to  form  03  [Copeland,  1994;  Copeland  et  al,  1994;  Knutsen 
et  al,  1994;  Yang  et  al,  1993]. 

To  test  the  potential  of  Reaction  (Rl)  as  the  missing  03  production  step,  we  examine  its 
kinetics  in  steady  state  at  maximum  dose  (103  km,  0.5  s).  Based  on  Knutsen  et  al.  (1994),  a 
suitably  conservative  estimate  for  the  rate  coefficient  for  Reaction  (Rl)  is  ~3  x  10'11  cm3/s. 

Given  the  02  number  density  at  103  km  (3.4  x  1012  cm'3)  and  the  apparent  bimolecular  03 
production  rate  of  3  x  108  cm'V1,  the  Oz(A)  number  density  would  have  to  be  about  3  x  106  cm'3 
in  the  beam-irradiated  volume.  At  this  altitude,  the  total  02(A3E)  loss  rate  due  to  radiation  and 
quenching  by  N2, 02,  and  O  (see  Table  1)  is  -400  s'1,  compared  to  the  characteristic  flow  loss 
of  ~2  s'1  (cf.  Fig.  3).  Thus  the  required  02(A)  production  rate  in  the  beam  would  have  to  be 
-1.2  x  109  cm'V1  or  5  x  10'3  molecules  per  ion  pair. 

This  value  is  consistent  with  the  expected  production  rate  of  02(A3E)  in  electron- 
irradiated  air.  Based  on  detailed  modeling  of  energetic  electron  deposition  in  air  and  the 
excitation  cross  section  given  by  Keto  (1981)  has  determined  an  02(A32)  excitation  rate  of  0.02 
molecules  per  ion  pair.  (We  note  that  Keto  (1981)  quoted  this  cross  section  for  excitation  of  the 
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c'Z  state,  but  the  same  cross  section  was  reported  earlier  for  excitation  of  the  A3E  state  by 
Watson  et  al.  (1967).)  Vertical  excitation  from  O2(X3Z,v=0)  populates  the  A-state  on  its 
repulsive  wall  above  the  dissociation  limit,  resulting  in  a  high  dissociation  probability.  Based  on 
Franck-Condon  considerations,  only  about  5  to  10%  of  the  direct  excitation  events  will  result  in 
bound  02(A3I),  with  initial  population  in  high  vibrational  levels.  This  suggests  a  molecular 
excitation  rate  of  =0.002  molecules  per  ion  pair. 

The  excitation  of  02(A3E)  can  be  significantly  augmented  by  the  presence  of  vibrationally 
excited  02(v*l),  which  has  high  Franck-Condon  overlap  with  bound  vibrational  levels  of 
02(A3L,v).  Using  the  Keto  (1981)  cross  section  and  the  computed  secondary  electron  energy 
distribution  described  above,  we  estimate  an  effective  rate  coefficient  of  5  x  10'12  cm3/s  for 
electron-impact  excitation  of  02(X3X,v=l)  to  02(A3Z,v=6-8).  If  we  assume  an  excitation  rate  for 
02(v)  °f  3  quanta  per  ion  pair.  The  main  collisional  loss  for  02(v)  is  deactivation  by  atomic 
oxygen.  Based  on  the  measured  rate  coefficients  for  the  high-pressure  limit  of  O  +  02 
recombination  [Hippier  et  al,  1990]  and  for  O  +  02  isotope  exchange  [Anderson  et  al,  1985],  this 
rate  coefficient  is  probably  no  larger  than  =6  x  10'12  cm3/s.  This  corresponds  to  high-altitude  loss 
rates  of  less  than  2  s'1,  i.e.  less  than  the  flow  loss  even  for  maximum  dose  conditions.  Thus  02(v) 
does  not  reach  steady  state  in  the  irradiated  volume.  For  the  maximum  dose,  a  time-averaged 
median  number  density  for  02(v=l)  is  =  1  x  10n  cm'3.  This  gives  an  02(A,v)  production  rate  of 
7.5  x  10®  cm'V1  or  0.003  molecules  per  ion  pair.  The  total  estimated  02(A)  production  rate  from 
direct  electron  impact  on  O2(v=0,l)  is  then  0.005  molecules  per  ion  pair,  in  excellent  agreement 
with  the  value  required  to  account  for  the  maximum  dose  observations. 
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More  detailed  kinetic  calculations,  performed  by  numerical  integration  of  the  coupled 
differential  equations  represented  by  the  reaction  mechanism  of  Table  1,  confirm  the  above 
analysis.  Example  time-dependent  calculations  for  the  maximum  dose  and  descent  conditions 
are  illustrated  in  Figure  11,  showing  the  effect  of  including  and  omitting  the  02(A)  chemistry.  In 
this  case  we  have  assumed  that  02(A)  reacts  with  02  to  form  O3(v=0),  and  that  all  the  vibrational 
excitation  of  03(v)  comes  from  electron  impact  according  to  the  laboratory  [Holtzclaw  et  al., 
1991, 1992]  vibrational  distribution.  The  maximum-dose  case  is  shown  in  Figure  11a.  In  the 
absence  of  the  02(A)  contribution,  total  ozone  in  the  beam-irradiated  gas  is  depleted  by  a  factor 
of  two,  and  03(v)  levels  are  less  than  half  the  measured  values.  When  production  of  O3(v=0)  by 
02(A)  is  included,  total  ozone  decreases  at  early  times  and  then  increases  substantially  as  the 
02(A)  number  density  increases  beyond  a  certain  level.  The  effect  of  the  02(A)  chemistry 
appears  only  after  about  100  ms  of  reaction  time,  owing  to  the  relatively  slow  production  rate  of 
02(A).  Thus  the  effect  should  not  be  observable  in  spectral  scans  with  shorter  beam  residence 
times,  as  on  descent.  As  shown  in  Figure  lib,  the  104  km  (descent)  scan  has  only  a  factor  of  two 
smaller  ion  pair  production  rate  than  the  maximum  dose  case,  but  the  beam  residence  time  of  less 
than  20  ms  prevents  02(A)  chemistry  from  contributing  to  the  observed  03(v)  excitation.  In  this 
case,  the  03(v)  enhancement  reflects  only  the  initial  excitation  rate  of  03(v)  by  secondary  electron 
impact. 

Figure  11 

The  vibrational  distributions  predicted  for  maximum  dose  are  compared  to  the  observed 
ones  in  Figure  12.  The  calculated  distributions,  which  are  essentially  those  observed  in  the 
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laboratory  for  secondary  electron  impact,  agree  with  the  observations  below  v=3,  but  greatly 
underpredict  the  excitation  for  v=4,5.  This  suggests  the  possibility  of  direct  formation  of 
03(v^3)  via  02(A),  or  possibly  the  energy  transfer  process 

02(b‘S)  +  O3(v=0)  -  02(X'Z)  +  03(v*5)  (R2) 

as  postulated  to  account  for  similar  vibrational  distributions  observed  in  laboratory  microwave 
discharges  [Rawlins  et  al,  1987]. 

Figure  12 

We  have  also  considered  other  possible  beam-initiated  reactions  for  formation  of  03.  The 
reactions  of  02(A)  with  vibrationally  excited  02(v)  and/or  the  metastable  02(a‘A), 


02(A)  +  02(v)  -*  O  +  03(v) 

(R3) 

02(A)  +  02(a'A)-0  +  03(v) 

(R4) 

have  the  necessary  energetics,  and  both  02(v=l)  and  02(a‘A)  are  formed  efficiently  and 
quenched  slowly  in  electron-bombarded  air.  However,  even  at  maximum  dose  these  species  can 
only  reach  number  densities  approximately  1/10  and  1/100  of  the  ambient  02  number  density,  so 
Reactions  (R3)  and  (R4)  would  require  rate  coefficients  near  or  in  excess  of  the  gas  kinetic  limit 
in  order  to  be  responsible  for  all  the  observed  03(v)  excitation.  While  Reaction  (R3)  could 
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possibly  contribute.  Reaction  (R4)  is  unlikely  to  play  a  role.  The  associative  detachment 
reaction. 


O'  +  02(v)  -  03  +  e‘  (R5) 

requires  02(v^3)  to  be  exoergic.  O'  is  produced  at  0.015  molecules  per  ion  pair,  and  is  rapidly 
destroyed  by  associative  detachment  with  O.  Thus  O'  and  02(v;>3)  are  not  likely  to  be  formed  in 
sufficient  abundance  to  produce  03  with  a  reasonable  rate  coefficient.  The  reaction 

O'  +  02(a'A)  ~  Oj  +  e'  (R6) 

has  recently  been  shown  not  to  occur  with  a  significant  rate  coefficient  [Upschulte  et  al.  1994b]. 

6.  Discussion  and  Conclusions 

In  summary,  the  EXCEDE  in  flight  experiment  obtained  high-altitude  data  under  well- 
controlled  dosing  and  measurement  conditions;  the  data  clearly  indicate  a  mechanism  for  the 
excitation  of  03(v)  in  electron-bombarded  air.  Most  of  the  excitation  is  caused  by  inelastic 
collisions  of  secondary  electrons  with  ozone  molecules.  Furthermore,  it  is  likely  that  additional 
03  is  formed  in  the  dosed  volume  by  reactions  of  electron-excited,  metastable  molecular  oxygen 
with  ambient  02.  This  interpretation  of  the  EXCEDE  in  flight  results  is  quantitatively  consistent 
with  laboratory  kinetic  measurements  of  the  excitation  processes. 
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To  further  test  the  comparison  between  the  flight  and  laboratory  results,  we  have 
examined  the  relative  03(v3)  and  C02(v3)  excitation  observed  by  EXCEDE  HI.  Comparisons  of 
the  03  radiances  near  10  fim  and  the  C02  radiances  at  4.3  fim,  accounting  for  the  differences  in 
the  ambient  03  and  C02  number  densities,  clearly  indicate  that  03  excitation  by  electron  impact 
is  roughly  an  order  of  magnitude  more  efficient  than  C02  excitation.  From  inspection  of  the 
beam-excited  C02(v3)  radiances  near  100  km,  we  estimate  the  ratio  of  the  03  and  C02  excitation 
rate  coefficients  to  lie  in  the  range 


7  ^  ^e+O3^e+C02  ^  20. 

This  is  quite  consistent  with  the  laboratory  value  for  this  ratio  of  12  ±  6,  determined  in  the 
LABCEDE  electron  beam  facility  at  =1  mtorr  and  100  K  [Holtzclaw  et  al,  1992].  Since  the 
C02(v3)  excitation  cross  section  peaks  at  =10'16  cm2  near  1  eV  [Lowke  et  al,  1973],  this  ratio 
indicates  a  peak  03(v3)  excitation  cross  section  near  =  10'15  cm2.  Thus  the  electron-impact 
excitation  of  ozone  is  a  highly  efficient  collisional  process,  perhaps  owing  to  the  propensity  of 
ozone  to  attach  low-energy  electrons. 

The  high  excitation  efficiency  of  03(v3)  applies  directly  to  the  auroral  excitation  scenario. 
For  an  IBC  ET  aurora  of  the  type  observed  in  the  HIRIS  rocket  flight  [Rawlins  et  al,  1985; 
Blumberg  et  al.,  1990],  i.e.  =600  kR  at  391.4  nm  near  1 10±10  km,  the  local  ion  pair  production 
rate  exceeds  106  cm'V1,  and  the  electron  number  density  exceeds  106  cm"3.  Under  these 
conditions,  the  excitation  of  03(v)  by  electron  impact  will  proceed  at  a  rate  near  105  cm'V,  or 
perhaps  faster  if  the  auroral  secondary  electron  distribution  is  more  energetic  than  that  produced 
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by  the  EXCEDE  HI  electron  beam.  This  excitation  rate  is  comparable  in  magnitude  to  that 
observed  for  the  unusually  intense  auroral  arc  in  the  HIRIS  flight  above  100  km  [Rawlins  et  al, 
1985].  This  excitation  rate  is  also  comparable  in  magnitude  to  the  production  rate  of  03(v=l)  by 
radiative  excitation  near  100  km,  and  is  well  in  excess  of  the  production  rate  of  03(v;>  1)  by  three- 
body  recombination  above  100  km  [Rawlins,  1985;  Rawlins  et  al,  1993].  Although  more 
detailed  modeling  of  the  auroral  excitation  process  is  warranted,  it  appears  that  direct  excitation 
of  Oj(v)  by  secondary  electron  impact  should  give  rise  to  observable  auroral  enhancements, 
especially  for  strong  auroras  at  altitudes  between  100  and  110  km. 

The  two-body  reaction  (Rl)  can  contribute  to  03  formation  in  both  the  auroral  and  the 
quiescent  upper  atmosphere.  In  auroras,  02(A,A',c)  is  produced  by  electron  impact  excitation  of 
ambient  02  and  amorally  generated  02(v)  as  outlined  above  for  the  electron  beam.  (However, 
due  to  moderately  efficient  quenching  of  02(v)  by  O,  the  02(v)  pathway  may  be  less  significant 
at  auroral  dosing  levels  and  electron  number  densities.)  02(A,A',c)  is  also  formed  by  three-body 
recombination  of  atomic  oxygen,  a  process  which  occurs  in  both  quiescent  and  auroral  condi¬ 
tions.  Ambient  02(A,A',c)  number  densities  maintained  by  O-atom  recombination  are  typically 
102  to  103  cm 3  near  100  km  [Allen,  1986;  Murtagh  et  al,  1986].  Using  the  rate  coefficient  for 
Reaction  (Rl)  discussed  above,  these  number  densities  will  give  two-body  03  production  rates  of 
roughly  104  to  105  cm'V  at  100-105  km.  These  values  are  comparable  to  the  three-body  (O  +  02 
+  M)  recombination  rate  at  this  altitude  [Rawlins,  1985;  Rawlins  et  al,  1993].  Thus  it  appears 
likely  that  metastable  02  via  Reaction  (Rl)  contributes  significantly  to  the  production  of  ozone 
near  100  km,  and  may  be  partly  responsible  for  the  unexpectedly  large  ozone  number  densities 
observed  at  these  altitudes  as  suggested  by  Allen  (1986).  Clearly,  more  detailed  modeling  of  this 
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complex  process  is  warranted,  in  order  to  re-examine  the  prevailing  views  of  odd-oxygen 
chemistry  in  the  lower  thermosphere. 
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Table  1.  List  of  Reactions  and  Rate  Coefficients 


Rate 


Reaction 

Coefficient 

Source 

03(v)  excitation:  electron  impact 

See  text;  LAB  CEDE  v-distribution 

1 

e+O3(0)  =  O3(l)+e- 

4.0  x  10'9  cm3/s 

2 

e  +O3(0)  =  O3(2)+e- 

9.0  x  10,0cm3/s 

3 

e'  +  O3(0)  =  03(3)+  e‘ 

1.7xlO,0cm3/s 

4 

e+O3(0)  =  O3(4)+e- 

3.5  x  10"  cm3/s 

5 

e'  +  O3(0)  =  03(5)+  e‘ 

7.5  x  1012  cm3/s 

03(v)  deactivation:  superelastic 

Reverse  of  excitation;  Te  =  0.5  eV 

6 

e+O3(l)  =  O3(0)+e- 

5.2  x  10‘9  cm3/s 

7 

e'  +  03(2)  =  O3(0)+  e' 

1.5  x  10'9  cm3/s 

8 

e'  +  03(3)  =  O3(0)+  e‘ 

3.5  x  10'10  cm3/s 

9 

e  +  03(4)  =  O3(0)+  e‘ 

9.4  x  10'“  cm3/s 

03(v)  dissociation:  electron  impact 

Cross  section:  Keto  (1981) 

10 

e'  +  O3(0)  =  e'  +  O  +  02 

2.8  x  10'°  cm3/s 

11 

e'  +  03(1)  =  e '  +  O  +  02 

4.1  x  10,ocm3/s 

12 

e  +  03(2)  =  e  +  O  +  02 

6.0  x  10,0cm3/s 

13 

e  +  03(3)  =  e"  +  0  +  02 

9.0  x  10'°  cm3/s 

14 

e'  +  03(4)  =  e  +  O  +  02 

1.4  x  10'9  cm3/s 

15 

e'  +  03(5)  =  e'  +  0  +  02 

2.1  x  10'9  cm3/s 

03(v)  deactivation:  radiative  cascade 

Rawlins  et  al.  (1987) 

16 

03(5)  =  03(4)+  hv 

41.4  s1 

17 

03(4)  =  03(3)+  hv 

35.8  s’1 

18 

03(3)  =  03(2)+  hv 

28.9  s’1 

19 

03(2)  =  03(1)+  hv 

20.7  s1 

20 

03(1)  =  O3(0)+  hv 

11.2  s"1 

03(v)  deactivation:  collisions  with  O 

Rawlins  et  al.  (1987) 

21 

03(5)  +  O  =  02  +  02 

3.6  x  10"  cm3/s 

22 

03(4)  +  O  =  02  +  02 

2.7  x  10  "  cm3/s 

23 

03(3)  +  O  =  02  +  02 

2.0  x  10  "  cm3/s 

24 

03(2)  +  O  =  02  +  02 

1.4  x  10  "  cm3/s 

25 

03(1)  +  0  =  03(0)+0 

8.0  x  1012  cm3/s 

03(v)  Deactivation:  N2, 

02  collisions 

Upschulte  et  al.  (1994a) 

26 

03(5)  +  M  =  03(4)  +  M 

1.0  x  1013  cm3/s 

27 

03(4)  +  M  =  03(3)  +  M 

8.0  x  10'14  cm3/s 

28 

03(3)  +  M  =  03(2)  +  M 

6.0  x  1014  cm3/s 

29 

03(2)  +  M  =  03(1)  +  M 

4.0  x  1014  cm3/s 

30 

03(1)  +  M  =  O3(0)  +  M 

2.0  x  10'14  cm3/s 
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Table  1.  List  of  Reactions  and  Rate  Coefficients  (Continued) 


Reaction 

Rate 

Coefficient 

Source 

31 

02(v),  02(A)  processes 
e'  +  02(0)  =  02(v)  +  e' 

3  (ion  pair)'1 

Gilmore  (1987);  v=l 

32 

e'  +  O2(0)  =  02(A)  +  e' 

0.002  (ion  pair)'1 

Gilmore  (1987);  likely  very  high  v 

33 

e'  +  02(v)  =  02(A)  +  e' 

5.0  x  10'12  cm3/s 

see  text;  likely  v=6-8 

34 

0  +  02(v)  =  O  +  02 

6.0  x  10'12cm3/s 

see  text 

35 

02(A)  +  N2  =  02  +  N2 

2.4  x  10 11  cm3/s 

Knutsen  et  al.  (1994) 

36 

02(A)  +  0  =  02  +  o 

9.0  x  1012  cm3/s 

Kenner  and  Ogryzlo  (1980) 

37 

02(A)  +  02  =  O3(0)+  0 

2.8  x  1011  cm3/s 

Knutsen  et  al.  (1994) 

38 

02(A)  =  02  +  hv 

6.67  s'1 

Saxon  and  Slanger  (1986) 
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Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


Figure  8. 


Figure  9. 


Figure  10. 


Electron  beam  viewing  geometry  (85  to  1 15  km). 

Volumetric  ion  pair  production  rates. 

Gas  residence  times  in  viewed  electron  beam  volume. 

Observed  03(v)  fluorescence  spectra,  (a)  electron  beam  off;  (b)  electron  beam  on. 

Example  of  spectral  fitting  analysis.  Observed  spectrum:  ascent,  102.5  km, 
electron  beam  on.  Computed  spectrum:  (00v3  +  10v3),  (01v3  +  llv3),  18-band 
model. 

03(v)  column  densities,  electron  beam  off. 

03(v)  column  densities,  electron  beam  on,  v=l-3.  Filled  symbols:  ascent;  open 
symbols:  descent.  Curves  are  curvefits  to  beam-off  data  (Fig.  6). 

03(v)  column  enhancement  by  electron  beam. 

03(v)  number  densities  excited  by  electron  beam  in  viewed  beam  volume. 

Estimated  ambient  03  and  O  number  densities,  compared  to  MSIS90  predictions. 
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Figure  11. 


Figure  12. 


Predictions  of  kinetic  model  for  time-dependent  beam  excitation,  (a)  103  km 
ascent;  (b)  104  km  descent.  The  solid  lines  indicate  predictions  including 
formation  of  03  by  02(A)  +  02;  the  dashed  lines  show  the  effects  of  omitting  that 
reaction  from  the  kinetic  scheme.  The  observed,  time-integrated  03(v^l)  number 
densities  are  also  indicated. 

Vibrational  distributions  of  03(v)  excited  in  electron  beam  at  maximum  dose. 

The  solid  and  dashed  lines  indicate  the  model  predictions  for  electron  impact 
excitation  with  and  without  formation  of  additional  03  by  the  02(A)  +  02  reaction. 
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Abstract 

A  high  flux,  wide  area  pulsed  hypersonic  molecular 
beam  facility  has  been  developed  for  the  study  and 
simulation  of  rarefied  flow  phenomenology.  The  beam 
source  operates  via  the  use  of  laser-induced  breakdown  in 
a  confined  volume  of  cold  gas.  The  device  is  small 
scale,  housed  in  a  40  cm  stainless  steel  vacuum  chamber, 
however  the  beam  area  can  be  as  large  as  1000  cm2. 

The  facility  can  be  used  to  study  gas-surface  and  gas-gas 
interactions  under  free  molecular  or  transition  flow 
regime  conditions. 

The  device  was  designed  initially  for  the  production 
of  an  oxygen  atom  beam,  tunable  over  the  velocity  range 
of  5  to  12  lcm/s.  This  beam  has  been  successfully  used 
to  study  material  erosion,  momentum  and  energy 
accommodation  on  surfaces,  contaminant  cleanup,  surface 
glows  and  collisional  excitation  of  gases.  The  source  has 
recently  been  extended  to  produce  beams  of  species  such 
as  N2,  C02  and  H20  in  the  velocity  range  of  2  to 
3.6  km/s.  This  extended  capability  significantly  broadens 
the  range  of  applications  for  the  facility,  for  example,  it 
can  be  applied  to  exhaust  plume  impingement  evaluations. 

The  operating  principles  and  characteristics  of  the 
facility  will  be  presented.  Key  facility  diagnostics,  such 
as  a  time  of  flight  mass  spectrometer,  quartz  crystal 
microbalances,  and  silver  film  scatterometers  will  be 
described,  and  representative  data  applicable  to  various 
rarefied  flow  phenomena  will  be  presented. 

Introduction 

We  developed  a  pulsed  fast  oxygen  atom  beam  in  the 
late  1980’s1  to  provide  a  ground  test  facility  to  qualify 
materials  to  be  used  in  low  earth  orbit  (LEO).  This  was 
in  response  to  a  number  of  early  space  shuttle 
observations2'6  of  material  erosion  and  property  change 
resulting  from  interactions  with  the  LEO  ambient 
atmosphere,  which  is  dominantly  atomic  oxygen.  These 
interactions  occur  at  orbital  velocities,  -  8  km/s,  and  for 
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a  while  there  was  a  frenzied  activity  to  develop  atomic 
oxygen  beams  which  operate  at  this  velocity.  The 
manifold  of  such  sources  has  been  reviewed  previously6 
and  to  our  knowledge  our  source  is  the  only  system  that 
provides  a  neutral  oxygen  atom  beam  with  sustained 
high  flux  and  wide  area  operating  capability  at  the 
desired  velocity  of  8  km/s.  Indeed,  in  several  tests 
groups  of  materials  having  dimensions  of  2.5  x  2.5  cm 
have  been  simultaneously  irradiated  to  8  km/s  oxygen 
atom  fluences  of  1022  cm'2. 

Over  the  years  this  fast  beam  technology  has  found 
many  applications  beyond  material  testing.  For 
example,  the  O-atom  beam  has  been  used  to  study  and 
simulate  the  visible  shuttle  glow,7'9  demonstrating  this 
glow  results  from  the  interaction  of  ambient  oxygen 
atoms  with  surface-adsorbed  NO.  Erosion  induced 
infrared  glows  have  also  been  measured. 10  The  O-atom 
beam  has  also  been  used  to  study  LEO  ambient 
atmosphere  scouring  of  surface  contaminants11  and  to 
measure  fast  oxygen  atom  surface  accommodation 
coefficients12  and  scattering  distributions.  Lastly,  the 
O-atom  beam  has  been  used  in  a  crossed  beam 
configuration  to  study  gas-gas  interactions  which  lead  to 
ultraviolet,  visible,  and  infrared  emissions. 13-16 

The  basic  fast  beam  technology  has  been  extended  to 
other  atomic  systems:  rare  gases  for  contaminant 
cleaning  applications;  nitrogen  atoms  for  kinetic  studies; 
and  chlorine  and  fluorine  atoms  for  semiconductor 
etching  applications.  Most  recently  the  technique  has 
been  extended  to  form  wide  area  beams  of  molecular 
species  for  plume  impingement  simulation. 

Several  of  these  activities  have  been  reviewed 
recently17  and  others  are  in  the  development  stage.  In 
this  paper  we  will  provide  an  overview  of  the  fast  beam 
technology,  including  a  review  of  diagnostic  capabilities, 
then  discuss  new  work  on  thin  silver  film  oxidation  and 
O-atom  scattering  distributions.  Lastly,  we  will  present 
preliminary  results  on  the  formation  of  lower  velocity 
molecular  beams. 
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The  Fast  Atom  Source 


Table  1.  Fast  O-Beam  Properties 


The  hyperthermal  oxygen  atom  source  and 
experimental  configuration  have  been  described  in  detail 
previously1'17  and  will  be  discussed  only  briefly  here.  A 
schematic  of  the  experimental  system  is  shown  in  Fig.  1. 
It  is  comprised  of  two  stainless  steel  six-way  crosses, 
including  an  8  in.  cross  source  chamber  housing  a  pulsed 
oxygen  valve/nozzle  assembly,  connected  to  a  16  in. 
cross  expansion  chamber.  The  system  is  pumped  by  a 
cryopump  (CVI,  TM250)  attached  to  the  large  cross 
which  maintains  a  base  pressure  of  3  x  KF7  Torr. 


The  oxygen  atoms  are  generated  in  a  pulsed  laser 
discharge  of  pure  02  as  described  by  Caledonia  et  al.1  A 
10  J/pulse  C02  laser  is  focused  with  a  100  cm  focal 
length  BaF2  lens  into  the  throat  of  a  20-deg  full  angle, 
12.5  cm  long  conical  expansion  nozzle  which  has  been 
partially  filled  with  by  a  pulsed  beam  valve  (General 
Valve).  The  plasma,  ignited  at  the  throat,  expands  out 
the  nozzle  dissociating  the  molecular  oxygen  in  front  of 
it.  The  nozzle  was  designed  to  allow  ion  and  electron 
recombination  while  the  slower  kinetics  of  atom/atom 
recombination  maintains  a  highly  dissociated  beam.  The 
result,  at  8  km/s,  is  an  approximately  50  fis  pulse  of 
highly  dissociated  oxygen  (>  80%  atoms)  with  less  than 
1  %  ion  content  (the  beam  is  charge  neutral).  The  general 
beam  properties  are  provided  in  Table  1  while  the 
diagnostics  used  in  evaluating  these  properties  are 
summarized  in  Table  2. 

The  beam  velocity  (5  to  12  km/s)  is  selected  by 
varying  the  time  delay  between  the  pulsing  of  the 
valve  and  the  triggering  of  the  CO^  laser.  The  delay 
determines  the  mass  of  Oj  processed  by  the  pulsed 
discharge  which  is  inversely  related  to  the  beam  velocity. 
The  velocity  is  determined  by  measuring  the  time  of 


Velocity 

8  km/s  ±  15  %  (5  to  12  km/s 
range) 

Fluence 

— 1018  O-atoms/pulse,  3  Hz 

Composition 

>80%  oxygen  atoms 

Size 

Expandable  to  >  1000  cm2  area 

Charge 

content 

<  1  %  ions  (controllable  by 
pseudo-Helmholtz  coils) 

Metastable 

content 

0(1D)  concentration  <0.4% 

Temperature 

T  *  300  K 

flight  of  the  pulse  between  two  radiometers  attached  to 
the  source  chamber,  mounted  along  the  beam  axis  and 
separated  by  7.6  cm,  and  filtered  to  monitor  the  5S-5P 
oxygen  atom  line  (777  nm)  (see  Fig.  1). 


Table  2.  O-Beam  Diagnostics 


Velocity,  temporal 
pulse  shape 

Radiometers 

Mass  spectrometer  (standard 
quadrupole  and  TOF) 

Pulse  fluence 

Mass  flow  meter 

Thin  silver  film  actinometer 

Beam  composition 

Mass  spectrometer 

Excited  state 
concentrations 

spectrometer 

Bandpass  filters 

Charge  content 
(beam  is  charge 
neutral) 

Langmuir  probe 

Beam  shape 

Erosive  witness  samples 
Positionable  Langmuir  probe 

2D  visible  imaging 

Beam  translational 
temperature 

Slitted  thin  silver  film  target 
chamber 

We  have  estimated  the  flux  of  oxygen  atoms 
generated  per  pulse  by  measuring  the  mass  of  O2 
processed  per  pulse.  This  was  determined  by  using  the 
known  time  delay  between  the  pulsing  of  the  oxygen 
valve  and  the  firing  of  the  laser,  monitoring  the  flow  of 
oxygen  with  a  mass  flow  meter  and  measuring  the 
temporal  profile  of  the  molecular  oxygen  pulse  (laser 
off)  with  a  pressure  transducer  aligned  along  the  beam 
axis.  With  this  measure  of  the  O2  processed  per  pulse, 
the  flux  at  any  point  down  stream  can  be  estimated  if  the 


2 

166 


expansion  characteristics  of  the  beam  are  known.  The 
expansion  properties  have  been  determined 
previously13,17  by  measuring  the  fast  O-atom  erosive 
mass  loss  of  polyethylene  samples  mounted  75  cm 
downstream  from  the  nozzle  throat  as  a  function  of  radial 
position.  The  expansion  half  angle  (to  1/e  times  the  on- 
axis  flux)  was  found  to  be  14  deg.  These  calibrations 
have  been  cross  checked  by  measuring  the  absolute 
oxidation  rates  of  materials  and  comparing  them  with 
equivalent  space  flight  measurements.  More  recently,  as 
discussed  in  the  next  section,  the  oxidation  rate  of  thin 
silver  films  has  been  used  as  a  calibration  standard. 


Although  the  beam  is  charge  neutral,  it  does  have  a 
small  but  finite  ion  content.  The  charge  concentration 
can  be  decreased  by  up  to  a  factor  of  eight  using  a  pair  of 
pseudo-Helmholz  coils  as  shown  in  Fig.  2.  Thus  any 
effects  due  to  charged  rather  than  neutral  species  can  be 
delineated.  The  local  charge  concentration  can  be 
measured  through  use  of  a  retractable  Langmuir  probe, 
Fig.  2.  The  probe  is  positionable  both  horizontally  and 
vertically  and  it  has  been  found  that  the  spatial  and 
temporal  shape  of  the  ion  distribution  matches  that  of  the 
neutral  distribution.17 


Fig.  2.  Schematic  of  beam  ion  control  system. 


We  have  recently  measured  the  0(1D)  concentration 
of  the  beam  by  monitoring  the  emission  from  the 
forbidden  0(JD  3P)  transition  using  a  specially 
constructed  bank  of  narrow-band  interference  filters.  The 

O^D)  concentration  was  found  to  be  £  0.4%  for  a 
beam  velocity  of  8  km/s  and  to  increase  monotonically 
with  beam  velocity  over  the  velocity  range  of  6  to 
12  km/s.  This  concentration  is  sufficiently  low  as  to 
render  metastable  effects  negligible  in  the  various  studies 
mentioned  above. 


We  have  also  evaluated  the  translational  temperature 
of  the  beam.  This  was  done  by  passing  the  O-beam 
through  a  0.15  cm  diameter  circular  aperture  and 
measuring  the  oxidation  pattern  on  a  thin  silver  film 
witness  plate  placed  5  cm  behind  the  aperture.  Any 
growth  in  the  beam  beyond  the  slit  width  is  due  to  the 
Maxwellian  motion  of  the  O  atoms  and  can  be  related  to 
an  effective  temperature  for  an  8  km/s  velocity  beam. 
This  temperature  was  deduced  to  be  —  300  K  at  a 
distance  of  85  cm  from  the  nozzle  throat.  The  beam 
temperature  will  be  higher  closer  to  the  throat  and  is 
calculated  to  be  —  1000  K  at  the  nozzle  exit  plane. 


Thin  Silver  Film  O-Atom  Diag 


There  has  been  a  need  for  some  time  for  a  general 
purpose  monitor  for  the  flux  of  atomic  oxygen.  Thin 
silver  films  were  a  prime  candidate  because  they  have 
the  virtue  of  being  oxidized  efficiently  by  oxygen  atoms 
but  not  by  oxygen  molecules.  A  number  of  the  film 
properties  change  upon  oxidation.  Specifically,  the  film 
goes  from  being  electrically  conductive  to  insulative,  the 
film  transmission  characteristics  for  visible  light  go  from 
opaque  to  transparent  and,  of  course,  the  film  gains 
mass.  Each  of  these  properties  may  be  exploited  to 
evaluate  the  flux  of  oxygen  atoms  impinging  the  film 
surface.  Early  uses  of  this  diagnostic  as  an  in-situ 
ambient  O-atom  diagnostic18"22  relied  upon  a  temporal 
monitor  of  the  resistivity  (conductivity)  of  the  films 
although  recently  the  use  of  silver-coated  quartz  crystal 
monitors  (QCM)  have  been  discussed.23 


Thin  silver  films  have  also  been  used  in  recent  years 
for  LEO  applications.  Gregory  and  Peters24  have  flown 
scatterometers  designed  to  monitor  the  angular 
distribution  of  ambient  oxygen  atoms  scattered  off  of 
targets  in  LEO.  These  scatterometers  use  a  thin  silver 
film  witness  plate  with  the  oxide  thickness  diagnosed  by 
visible  light  transmission  measurements.  Cross  and 
Blais25  have  developed  silver  actinometers  to  monitor 
erosion  rates  and  to  calibrate  their  laboratory  fast  oxygen 
atom  beam  (2  to  3  eV).  As  in  the  earlier  work,  these 
actinometers  use  the  change  in  resistance  in  a  thin  silver 
film  to  deduce  impinging  oxygen  atom  fluxes. 


Prior  to  our  work,  there  does  not  appear  to  have 
been  any  systematic  study  of  the  phenomenology  of  the 
interaction  between  energetic  oxygen  atoms  and  thin 
silver  films.  In  particular,  for  high  velocity  interaction, 
we  could  locate  no  information  on  the  oxidation 
mechanism,  its  efficiency  and  velocity  and  temperature 
dependence,  and  on  the  effects  of  surface  contamination 


3 

167 


and  film  thickness.  We,  therefore,  undertook  such 
investigations  using  our  fast  oxygen  atom  source.  The 
results  of  this  study  are  described  in  detail  elsewhere26 
and  summarized  briefly  below. 

We  have  employed  our  fast  O-atom  beam  to  study 
oxidation  rates  of  thin  silver  films  vacuum  deposited  on 
crystals  (various  thicknesses  of  up  to  1  fim)  and  on 
kapton  sheets  (600 A  A g  films).  The  resistance  method 
was  used  to  monitor  the  oxidation  rate  of  the  latter  films 
but  the  data  scatter  was  so  large  that  our  analysis  was 
confined  to  the  QCM  results. 

Oxidation  of  the  silver  QCM  ciystal  results  in  a  mass 
gain  and  a  corresponding  decrease  in  the  oscillator 
frequency.  The  QCM  used  in  this  work  (from  Sycon 
Instruments)  incorporates  a  6  MHz  quartz  crystal  with  a 
sensitivity  of  1.25  x  10r8  g/cm2/Hz.  A  frequency  shift  of 
1  Hz  therefore  corresponds  to  the  oxidation  of  0.8 A  of 
silver  to  AgO.  The  mass  gain  during  oxidation  was 
monitored  by  sending  the  output  of  the  QCM  controller 
(voltage  proportional  to  mass  loading)  to  a  Data 
Translations  D/A  board  for  collection  and  monitoring  in 
real  time  with  a  personal  computer. 

The  QCM  was  mounted  directly  onto  a  temperature- 
controlled  aluminum  block  and  thus  the  quartz  crystal, 
which  is  in  intimate  thermal  contact  with  the  QCM  base, 
was  maintained  to  ±0.5°C.  Temperature  regulation  of 
the  QCM  was  necessary  since  the  oscillator  frequency  is 
temperature  dependent.  Therefore,  a  change  in 
temperature  could  be  erroneously  interpreted  as  a  mass 
gain  or  loss.  Near  20°C,  the  temperature  dependence 
results  in  a  frequency  shift  of  +0.2  Hz/deg  for  the  AT- 
cut  of  this  QCM  crystal  (35  deg,  11  min),  but  at  60°C, 
the  sensitivity  is  +6  Hz/deg. 

Typical  oxidation  histories  for  irradiation  by  8  km/s 
oxygen  atoms  at  room  temperature  are  shown  in  Fig.  3. 
Shown  is  oxidation  depth  versus  the  number  of  impinging 
O-atom  pulses  as  determined  by  the  measured  mass 
increase  for  two  film  thicknesses.  Two  important 
features  can  be  seen.  First,  no  oxidation  occurs  during 
the  first  two  hundred  pulses.  This  is  typically  observed 
and  we  relate  it  to  the  required  removal  of  a 
contamination  layer  on  the  silver  film.  Indeed,  if  we  first 
irradiate  the  film  with  several  hundred  pulses  of  a  fast 
argon  beam  and  then  shortly  thereafter  impinge  it  with 
fast  oxygen  atoms,  we  find  that  oxidation  initiation  is 
immediate.  Second,  after  oxidation  initiation,  the 
oxidation  depth  increases  linearly  with  O-atom  fluence  for 
approximately  the  first  250A  of  the  silver  film  and  then 
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Fig.  3.  Oxidation  behavior  of  thin  silver  films 
coated  on  QCM  crystals  (O-atom  velocity  =  8  km/s, 
room  temperature). 

the  oxidation  rate  slows  down.  We  believe  the  oxidation 
rate  becomes  diffusion  controlled  at  these  thicknesses 
and  thus  only  very  thin  films  are  useful  as  a  linear 
calibration  source.  Note  that  the  observed  oxidation  rate 
continues  to  decrease  very  dramatically  in  the  case  of  the 
1  fim  thick  film. 

Note  that  in  the  interpretation  of  the  QCM  data  it  is 
important  to  understand  the  oxidation  process.  During 
the  course  of  our  work  we  realized  that  the  poly¬ 
crystalline  silver  oxidized  to  the  peroxide,  AgO,  rather 
than  Ag20.  This  was  determined  by  measuring  the  mass 
uptake  of  a  totally  oxidized  film  with  known  silver  load¬ 
ing,  and  verified  both  by  heating  the  film  to  200°C  to 
form  Ag20  and  by  X-ray  diffraction  analysis  of  the  film. 

We  also  examined  the  temperature  dependence  of  the 
silver  oxidation  rate  over  the  range  of  0°  to  85  °C. 

These  results  are  shown  in  Fig.  4,  contrasted  with 
simultaneous  measurements  performed  on  a  poly- 
crystalline  carbon  film  also  coated  on  a  QCM  crystal 
(for  convenience,  all  results  are  normalized  to  a  silver 
reactivity  of  unity  at  273  K).  Although  there  is  some 
scatter  in  the  data,  the  silver  oxidation  rate  is  inde¬ 
pendent  of  temperature,  unlike  that  for  carbon.  Thus, 
we  believe  that  these  silver  films  will  provide  a  good 
standard  for  relative  investigation  of  high  velocity 
oxidation  behavior  of  materials.  Indeed,  we  have  found 
that  the  silver  film  oxidation  efficiency  is  also  inde¬ 
pendent  of  O-atom  velocity  over  the  velocity  range  of 
8  to  12  km/s. 

We  now  have  developed  a  clear  understanding  of  the 
behavior  of  thin  silver  films  subjected  to  high  velocity 
oxygen  atom  irradiation  and  believe  that  they  will  prove 
a  valuable  diagnostic  tool  within  their  demonstrated 
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Temperature  (K) 
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Fig.  4.  Measured  silver  and  carbon  reactivity 
vs  temperature  (O-Atom  velocity  8  kan/s). 


limitations.  One  such  diagnostic  application,  the 
measurement  of  the  scattering  distribution  for  high 
velocity  oxygen  atoms  impinging  on  selected  surfaces,  is 
described  in  the  next  section. 

O-Atom/Surface  Scattering  Experiments 

Knowledge  of  the  scattering  behavior  of  ambient 
oxygen  atoms  striking  spacecraft  surfaces  in  LEO  is 
important  for  two  reasons:  1)  the  degree  of  elasticity  of 
the  collisions  affects  spacecraft  lift  and  drag;  and  2)  the 
scattered  atoms  can  cause  deleterious  effects  to  materials 
not  impinged  by  the  ram  flow. 

Silver  film  scatterometers  designed  to  measure  the 
atomic  oxygen  scattering  distribution  off  various 
materials  have  been  flown  on-orbit  in  shuttle  experiments 
by  John  Gregory24  and  other  investigators.  These  on- 
orbit  measurements  show  that  the  atomic  oxygen 
scattering  is  non-specular.  Our  laboratory  measurements 
represent  the  first  attempt  to  confirm  these  observations 
in  the  laboratory  in  an  8  km/s  atomic  oxygen  flow.  The 
representative  materials  chosen  for  this  study  are  clear 
and  black  anodized  aluminums  provided  by  NASA/ 
Johnson  Space  Center,  and  microscope  slide  glass.  These 
materials  are  relevant  as  most  spacecraft  are  constructed 
from  aluminum  to  minimize  mass,  and  the  largest 
exposed  area  on  a  spacecraft  are  the  glass  covered  solar 
arrays.  We  observed  that  scattering  off  all  three 
materials  is  non-specular.  The  details  of  the 
scatterometer  design,  silver  film  cleaning  procedures,  and 
atomic  oxygen  scattering  measurements  are  described 
briefly  below. 


surfaces.  In  this  device  the  beam  passes  through  a  small 
slit  into  a  chamber  containing  the  target  material.  The 
O-atoms  scattered  from  the  target  ultimately  interact  with 
a  silver  film  "witness  plate"  positioned  at  the  front  of  the 
chamber  (around  the  slit).  The  distance  between  target 
and  slit  must  be  large  enough  to  allow  sufficient  spread 
of  the  scattered  beam  to  provide  good  measurement 
resolution  on  the  witness  plate  for  the  scattered  angular 
distribution.  This  measurement  is  complicated  however, 
in  that  the  beam  does  grow  in  extent  between  the  slit  and 
the  target  because  of  its  random  or  Maxwellian  velocity 
component.  The  ratio  between  the  thermal  and  directed 
velocities  defines  a  beam  spread  angle  which  must  be 
known  in  order  to  interpret  the  witness  plate  data.  All 
this  information  is  available  for  our  device  as  described 
previously. 


We  have  designed  and  built  silver  film 
scatterometers  specific  to  our  experiment.  Two 
scatterometers  were  constructed  from  6061-T6  aluminum 
to  hold  up  to  three  samples  in  individual  compartments 
at  either  a  90  or  45  deg  angle  of  incidence  to  the  8  km/s 
atomic  oxygen  flow.  The  90  deg  scatterometer  is  shown 
in  Fig.  5.  The  silver  film  to  be  oxidized  by  the 
scattered  atomic  oxygen  was  attached  to  the  2.222  cm 
radius  surface  by  adhesive  which  securely  held  the  film 
flat,  and  yet  permitted  easy  removal.  A  0.25  cm  wide 
by  1.25  cm  long  slit  was  cut  in  the  film  prior  to 
mounting  to  restrict  the  atomic  oxygen  flow  into  the 
scatterometer.  The  glass  and  aluminum  samples  were 
cut  to  fit  the  individual  sample  chambers  and  attached  to 
the  silver  coated  kapton  covered  rear  surface  with 
double  faced  tape.  A  key  feature  of  our  scatterometers 
not  found  in  other  designs  is  the  addition  of  a  hinged 
cover  to  permit  cleanup  via  argon  atom  bombardment  of 
both  the  sample  and  silver  film  surfaces  to  remove  the 
surface  contamination  layers  in  the  chamber  immediately 
prior  to  the  oxygen  scattering  measurement.  We  have 


We  used  a  modified  version  of  Gregory’s 
scatterometer  in  our  laboratory  experiments  to  measure 
the  oxygen  atom  angular  scattering  distributions  from 


Fig.  5.  Normal  incidence  spectrometer  used  to  measure 
distribution  of  fast  oxygen  atoms  scattered  off  of 
selected  materials. 
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found  that  this  was  a  key  feature  for  successful 
accomplishment  of  the  experiments. 

Our  witness  samples,  500A  silver  film  on  2  mil 
mylar,  were  prepared  by  Sheldahl.  In  this  application, 
the  films  were  analyzed  for  oxidation  by  measuring 
optical  transmission  with  a  computer  page  scanner 
equipped  with  a  transmission  attachment.  As  mentioned 
earlier,  silver  oxide  is  transparent  although  silver  is 
opaque.  The  fast  O-atom  beam  was  used  to  calibrate  the 
film  transparency,  e.g.,  optical  thickness  versus  atomic 
oxygen  fluence,  and  the  scattering  measurement 

experiment  duration  was  limited  to  ensure  that  <  250A 
depth  of  the  film  was  oxidized,  i.e.,  to  the  range  of  linear 
oxidation  behavior. 

The  resulting  measured  scattered  oxygen  atom 
distributions  versus  angle  from  the  slit  are  shown  in 
Figs.  6  and  7  for  normal  and  45  deg  incidence 
respectively.  The  results  are  averaged  over  the  width  of 
the  witness  film  and  are  presented  as  a  percentage  of  the 
incident  flux  per  degree.  The  angular  resolution  of  the 
data  is  about  2.5  deg.  The  absolute  calibration  for  this 
data  is  uncertain  to  about  a  factor  of  two  but  this 
uncertainty  can  be  greatly  reduced  in  future  work  by 
simultaneously  monitoring  the  O-atom  flux  during  the 
scattering  measurement. 


Angle  from  Slit  (deg) 

C-5282 


Fig.  6.  Oxygen  scattering  distribution  at  0  deg  angle 
of  incidence  for  three  materials.  O-atom  velocity  is 
8  km/s.  Low  measurements  at  positive  angles  on  black 
anodized  aluminum  are  due  to  improper  witness  film 
cleaning. 

The  scattering  distribution  is  quite  flat,  suggesting 
isotropic  scattering  for  all  three  materials.  No  strong 
specular  regions  are  apparent.  The  central  spike  that  is 
seen  in  the  figures  should  be  ignored.  It  is  due  to  the 
presence  of  a  slit  in  the  silver  film  and  the  fluence 
indicated  is  an  artifact  of  the  analysis  technique. 


Angle  from  Slit  (deg) 
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Fig.  7.  Oxygen  scattering  distribution  at  45  deg  angle 


of  incidence  for  three  materials.  V  =  8  km/s. 


The  importance  of  working  with  clean  films  was 
accidentally  demonstrated  in  the  data  for  black  anodized 
aluminum  shown  in  Fig.  6.  It  turns  out  that  the  door  of 
the  scatterometer  opened  too  far  during  the  exposure  of 
the  black  anodized  sample  and  one  side  of  the  silver  was 
not  scoured  by  the  atom  beam.  Thus,  the  anomolously 
low  apparent  scattered  fluence. 


There  is  one  interesting  observation  on  the  data 
presented  in  Figs.  6  and  7  which  can  only  be  treated  as 
a  conjecture  at  present.  That  is  that  the  apparent 
fractional  scattering  of  oxygen  atoms  is  largest  for  black 
anodized  aluminum,  smaller  for  clear  anodized 
aluminum  and  smallest  for  glass.  The  difference  is 
small  in  the  case  of  normal  incidence  but  the  scattered 
fluences  differ  by  20%  in  the  45  deg  incidence 
measurements.  These  measurements  were  taken  to 
evaluate  the  scattering  distribution,  not  the  magnitude  of 
the  scattered  fluence,  and  thus  the  incident  fluence  was 
not  highly  controlled  (  — 10%).  Nevertheless, 
considering  the  order  the  data  was  taken  in,  it  would  be 
highly  coincidental  that  the  apparent  scattering  efficiency 
of  the  three  materials  would  be  artificially  preserved  at 
two  different  impingement  angles.  If  this  is  a  true 
variation,  it  could  imply  a  higher  efficiency  of  O-atom 
recombination  on  glass  than  on  aluminum. 


We  note  that  an  isotropic  scattering  distribution  need 
not  imply  high  surface  accommodation.  The  roughness 
of  the  surfaces  tested  could  provide  for  specular 
scattering  over  a  wide  range  of  effective  angles  of 
incidence.  The  nature  of  the  scattering  is  best  deduced 
by  measuring  the  velocity  of  the  scattered  atoms  using  a 
time-of-flight  mass  spectrometer  (TOFMS). 


Our  approach  was  to  mount  a  Comstock  system  Eli 
mass  spectrometer  onto  a  rotating  ISO  400  flange  on  the 
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top  of  our  FAST-2  test  chamber.  The  experiment  was 
configured  to  provide  sufficient  ion  signal  based  upon 
expected  scattering  densities  compared  to  the  operating 
background  pressure  of  3  x  10*7  Torr. 

Experiments  were  performed  to  observe  scattered 
atomic  oxygen  after  impact  with  a  sample  of  clear- 
anodized  aluminum  as  provided  by  NASA  Johnson.  The 
experimental  setup  is  shown  in  Fig.  8.  The  TOFMS  was 
mounted  on  the  top-most  flange  of  a  40.6  cm  diameter 
six-way  cross  which  comprises  the  target  chamber  in  the 
FAST-2  atomic  oxygen  test  facility.  The  TOFMS 
ionization  region,  flight  tube,  and  detection  region  were 
differentially  pumped  with  a  turbo-molecular  pump  to 
pressures  of  *  3  x  1(T7  Torr  (ultimately  determined  by 
the  background  pressure  in  the  FAST-2  facility).  The 
TOFMS  was  mounted  13.3  cm  off  centerline  on  a  flange 
modified  to  allow  360-deg  rotation  around  the  vertical 
centerline  of  the  six-way  cross.  Two  apertures  in  the 
shroud  surrounding  the  ionization  region  of  the  TOFMS 
were  aligned  to  define  an  axis  through  the  center  of  the 
chamber.  The  use  of  two  apertures,  one  on  either  side  of 
the  ionization  region,  allows  passage  of  incoming  species 
through  the  mass  spectrometer.  This  reduces  the 
likelihood  of  a  density  gradient  occurring  in  the  region 
around  the  ionizer  as  would  almost  certainly  occur  in  the 
absence  of  a  second  aperture. 


Turbo 
Molecular 
Pump  — 
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Translating 
Rotating 
Shaft 


Time-of-Flight 
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Fig.  8.  Experimental  arrangement  for  measurement 
of  O-atom  scattering  with  a  time-of-flight 
mass  spectrometer. 


The  target  was  mounted  75  cm  from  the  O-atom 
source  in  the  center  of  the  chamber  on  a  shaft  which 
could  be  translated  along  the  vertical  axis  and  rotated 
with  respect  to  the  O-atom  beam.  Consequently,  the 
target  could  be  removed  from  the  O-atom  beam  for 


measurement  of  background  mass  spectra.  The  angle  of 
impact  for  the  atomic  oxygen  could  be  varied  from  0  to 
90  deg  by  rotation  of  the  target.  Since  the  TOFMS  was 
always  aligned  with  a  path  through  the  center  of  the 
chamber,  the  center  of  the  target  was  always  observed 
by  the  TOFMS  regardless  of  target  rotation  or  placement 
of  the  TOFMS  with  respect  to  the  O-atom  beam 
centerline.  The  target  was  oversized  at  5  x  7.5  cm  to 
provide  ease  in  aligning  the  spectrometer.  The  center  of 
the  TOFMS  ionizer  was  positioned  12.5  cm  from  the 
target. 

To  minimize  contamination  from  test  chamber- 
scattered  oxygen  atoms,  an  aluminum  plate  with  two 
apertures  3.2  cm  in  diameter  was  placed  42  cm 
downstream  from  the  throat  of  the  O-atom  source. 

These  apertures  served  to  limit  oxygen  flow  into  the 
target  chamber  and  thus  limit  background  signal.  One 
aperture  was  placed  directly  on  centerline  to  allow 
oxygen  to  impact  the  target,  the  second  was  placed  to 
allow  laser  access  to  the  nozzle  throat.  The  O-atom 
source  chamber  was  equipped  with  a  second  cryopump 
to  allow  for  efficient  removal  of  the  bulk  of  the  oxygen 
between  laser  shots. 

Since  the  ionization  efficiency  in  any  mass 
spectrometer  is  small,  the  minor  ion  content  of  the 
O-atom  beam  can  completely  obscure  any  signal 
generated  by  ionization  of  atomic  oxygen  neutrals.  Ions 
were  prevented  from  entering  the  mass  spectrometer  by 
placing  a  90%  transmitting  tungsten  mesh  screen 
(~0.3  mm  mesh  pitch)  in  front  of  each  aperture.  A 
bias  of  +20V  was  found  to  be  sufficient  to  eliminate  the 
ion  signal  from  the  beam. 

Two  sets  of  experiments  have  been  performed.  In 
the  first,  the  target  was  oriented  perpendicular  to  the 
O-atom  source  and  the  TOFMS  was  positioned  40  deg 
off-axis  from  the  beam.  The  second  experiment  was 
performed  with  the  target  tilted  45  deg  and  the  TOFMS 
oriented  90  deg  with  respect  to  the  O-atom  beam. 
Consequently,  to  the  extent  that  surface  roughness  does 
not  blur  the  angular  scattering  of  the  O-atoms,  the 
second  experiment  was  arranged  to  maximize  the 
probability  of  detecting  "specularly"  scattered  atoms. 

The  results  from  the  first  experiment  described 
above  are  shown  in  Fig.  9  where  both  the  mass  16  and 
mass  32  (O  and  O2)  signals  are  displayed.  For  mass  16, 
the  indicated  maximum  number  density  occurs  at  a 
delays  of  *  300  /xs  after  the  laser  fires  or  approximately 
210  ^s  after  the  oxygen  impacts  the  target.  This 
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Delay  After  Atom  Impact  on  Surface  (ps)  C  528s 


Fig.  9.  Temporally  resolved  signals  from  atomic  oxygen 
and  molecular  oxygen.  The  target  was  positioned 
normal  to  the  O-atom  pulse  and  the  TOFMS  was  placed 
40  deg  off-normal  with  respect  to  the  target.  The 
O-atom  nominal  velocity  was  8  km/s  corresponding  to  a 
time-of-arrival  of  95  /is  after  the  laser  shot. 

corresponds  to  a  most  probable  O-atom  velocity  of  7  x 
104  cm/s  or  a  kinetic  energy  of  «  0.04  eV.  It  should  be 
noted  that  the  mass  16  signal  was  quite  prominent  at  its 
maximum  being  *  5  times  larger  that  the  signal  present 
with  the  target  retracted.  The  corresponding  mass  32 
data  show  a  most  probable  time-of-arrival  of  ^  300  /is 
although  this  is  less  well  defined. 

These  data  are  consistent  with  the  hypothesis  that  the 
incoming  O-atom  kinetic  energy  is  largely  accommodated 
by  the  aluminum  surface  and  that  the  atoms  leave  the 
surface  at  near  thermal  velocities.  The  large  mass  32 
signal,  which  indicates  a  similar  maximum  number 
density  of  Oj  when  the  larger  O2  ionization  cross-section 
is  accounted  for  (a  factor  of  —  1.6),  further  suggests  a 
significant  probability  for  surface-assisted  O-atom 
recombination  on  the  aluminum  since  the  O2  content  in 

the  incident  beam  is  £  20%. 


Similar  results  are  observed  for  signals  generated 
with  the  target  tilted  45  deg  and  the  mass  spectrometer 
collection  axis  normal  with  respect  to  the  O-atom  beam. 
In  neither  experiment  was  a  prominent  signal  for  mass 
16  observed  at  the  short  times  consistent  with  elastic 
scattering  from,  the  aluminum.  These  experiments  do 
suggest  that  a  small  fraction  of  the  atoms  scattered  from 
the  surface  leave  with  velocities  consistent  with  near- 
elastic  scattering.  For  example,  a  signal  for  16, 
significantly  above  that  observed  with  the  target 
recessed,  was  observed  only  45  /is  after  the  nominal 
time  of  O-atom  impact.  This  corresponds  to  a  leading 
velocity  of  3  x  105  cm/s. 

In  summary,  these  experiments  suggest  that:  1)  the 
great  majority  of  atomic  oxygen  impacting  the  aluminum 
was  accommodated  by  this  surface;  2)  there  is  a 
substantial  range  in  the  extent  to  which  the  O-atom 
kinetic  energy  is  accommodated  (based  on  the  wide 
range  of  indicated  leaving  velocities);  and  3)  there  is 
substantial  O-atom  recombination  occurring  on  the 
aluminum  surface  as  indicated  by  the  prominent  mass  32 
signals. 

While  these  data  suggest  the  preceding  conclusions, 
they  are  by  no  means  definitive  and  they  must  be  viewed 
as  preliminary.  There  are  several  uncertainties  in  the 
experiment  which  could  greatly  influence  their 
interpretation.  First,  the  extent  to  which  species  leaving 
the  surface  are  detained  in  the  mass  spectrometer  (i.e., 
the  extent  to  which  they  "rattle  around"  before  being 
ionized)  is  not  presently  known.  The  layout  of  the 
experiment  allows  this  possibility  since  the  target  was 
large  enough  (*  5  x  8  cm)  to  produce  scattered  species 
with  an  angle  of  incidence  sufficiently  oblique  onto  the 
TOFMS  to  preclude  clean  passage  on  one  pass  through 
the  instrument.  This  "piling  up"  in  the  TOFMS  would 
have  the  effect  of  lengthening  the  duration  of  the  O-atom 
pulse.  A  second  effect  would  be  an  apparent  increase  in 
the  02  content  since  collisions  with  the  interior  TOFMS 
surfaces  would  promote  O-atom  recombination  into  O2. 
These  uncertainties  will  be  evaluated  in  future  work. 

In  any  event,  the  scatterometer  and  TOFMS 
measurements  are  consistent  with  each  other  both  in 
terms  of  the  clear  inelasticity /isotropicity  of  the 
scattering  and  also  in  terms  of  the  evidence  for 
significant  surface  induced  atomic  recombination  in  these 
high  velocity  gas  surface  interactions.  It  is  clear  that  the 
combination  of  these  two  diagnostics  provides  a 
powerful  tool  for  the  investigation  of  such  scattering 
phenomena. 
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Lower  Velocity  Molecular  Beams 


20  cm 


The  last  topic  in  this  paper  is  a  preliminary  study  of 
the  extension  of  our  laser  breakdown  technique  to  form 
beams  of  molecular  species  at  lower  velocities,  e.g.,  2  to 
4  km/s.  This  effort  was  originally  driven  by  the  desire  to 
simulate  exhaust  plume  impingement  effects.  The  effect 
of  exhaust  plume  impingement  on  spacecraft  surfaces  has 
been  an  issue  for  many  years.  The  need  to  evaluate  these 
effects  at  relatively  large  plume  expansion  distances  and 
low  pressures  has  required  the  use  of  costly  test  facilities 
with  high  pumping  capacity.  In  this  work,  we  have 
investigated  the  potential  of  using  the  laser  breakdown 
technique  to  develop  a  small-scale  facility  which  provides 
beams  of  exhaust  species  that  can  be  directed  onto 
surfaces.  This  beam  can  be  operated  over  a  range  of 
density/velocity  to  study  plume  impingement  effects. 

Our  effort  was  directed  towards  the  design, 
construction,  implementation,  and  testing  of  a  system  that 
would  generate  a  pulsed  beam/jet  that  simulates  the 
3.5  km/s  exhaust  of  a  hypergolic  bipropellant  thruster. 

Our  ultimate  goal  would  be  to  obtain  a  velocity  controlled 
beam  with  the  mole  fractions  of  species  listed  in 
Table  3. 5  Single  species  simulation  is  critical  for  plume 
impingement  model  evaluation  and  development  and  this 
was  investigated  first. 

Table  3.  Example  Shuttle  Thruster  Exhaust  Species 


Species 

Mole  Fraction 

h2o 

0.33 

n2 

0.31 

h2 

0.17 

CO 

0.13 

C02 

0.036 

H 

0.015 

°2 

0.0004 

A  design  modification  of  our  usual  breakdown  nozzle 
was  required  for  a  multi  component,  slow  velocity  beam. 
The  existing  nozzle  allows  injection  of  only  one  gas 
stream.  This  severely  restricts  synthesis  of  the  proper 
gas  mixture  by  requiring  sources  of  premixed  gases  for 
the  injector.  We  therefore  modified  the  usual  nozzle  to 
accept  a  second  pulsed  valve  and  gas  injector  as  shown  in 
Fig.  10.  This  allowed  us  to  controllably  mix  gases  in  the 
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Fig.  10.  Schematic  of  modified  nozzle. 


nozzle.  The  relative  fraction  of  the  two  gas  streams 
could  be  controlled  on  a  pulse-to-pulse  basis  for  testing 
of  exhaust  beam  relative  mole  fractions.  This  design  is 
not  optimum  for  good  gas  mixing  as  discussed  later.  A 
more  appropriate  design  would  have  both  valves  inletting 
into  the  nozzle  throat  area. 

Our  usual  nozzle  length,  and  laser  power  in 
conjunction  with  the  maximum  operating  pressure  of  the 
pulsed  valve  restricts  the  minimum  operating  velocity  of 
the  breakdown  source  to  about  5  km/s.  In  order  to  slow 
the  exhaust  beam  even  further  it  is  necessary  to  add 
more  gas  or  reduce  the  laser  power.  Reducing  the  laser 
power  translates  into  a  reduced  flux  and  density  of  the 
exhaust  beam.  We  chose  to  keep  the  laser  power  high 
to  generate  the  maximum  exhaust  density  for  plume 
impingement  studies,  and  thus  lengthened  the  existing 
nozzle  so  as  to  add  more  gas  to  the  breakdown.  A 
cylindrical  nozzle  extension,  20  cm  long  was  added  to 
the  12.7  cm  long  10  deg  half  angle  nozzle  to  confine  the 
longer  gas  pulses.  The  modified  nozzle  is  shown 
schematically  in  Fig.  10. 

The  standard  velocity  diagnostic  for  the  fast  O  beam 
source  is  a  pair  of  radiometers  viewing  the  exhaust  at 
known  distances  downstream  of  the  nozzle  exit  plane. 
Visible  emission  arising  from  the  ion-electron 
recombination  of  residual  charge  is  observed  as  the 
exhaust  beam  passes  the  radiometer’s  narrow  field  of 
view.  A  simple  timing  analysis  of  the  emission  arising 
in  the  radiometers  in  combination  with  a  known 
separation  between  the  radiometers  allows  the  beam 
velocity  to  be  measured.  This  technique  lacks  the 
sensitivity  required  to  detect  the  beam  passage  in  the 
current  configuration.  First,  slower  beams  produce  less 
ionization  in  the  beam  exhaust,  which  translates  into 
weaker  visible  emissions,  and  second,  the  nozzle 
extension  quenches  some  of  the  ions  and  electrons, 
which  further  reduces  the  intensity  of  the  visible 
emissions. 

For  the  velocities  of  interest  in  this  work,  we  found 
that  piezo-electric  dynamic  pressure  transducers  are  the 
preferred  diagnostic.  The  velocity  evaluated  in  this 
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Relative  Signal  (V) 


manner  was  cross-checked  through  measurements  of  the 
time  of  arrival  of  the  pulse  at  the  mass  spectrometer 
positioned  at  the  end  of  the  test  chamber. 

We  concentrated  on  generating  3  to  4  km  s’1  exhaust 
plumes  containing  the  relevant  species  in  approximately 
the  ratio  predicted  by  equilibrium  calculations  of  the 
combustion  products  in  a  hypergolic  bi-propellant 
thruster.  Beams  containing  individual  constituents  were 
also  generated.  A  quadrupole  mass  spectrometer 
(Balzers)  was  used  to  evaluate  the  species  components  of 
the  beam.  Experiments  were  performed  with  single  and 
double  gas  injection  into  the  breakdown  nozzle,  and  with 
liquid  water  injection  at  the  nozzle  throat.  Pure  gas 
experiments  using  N2  and  C02  exhibited  more  than  88  % 
molecular  species  in  the  plume  at  velocities  around  3.0  to 
3.6  km  s'1.  A  typical  mass  spectrometric  scan  for  the 
C02  beam  is  shown  in  Fig.  11.  The  beam  width  is  seen 
to  be  several  hundred  fis.  This  result  is  typical  for  the 
pure  gas  experiments.  Injection  of  liquid  water  at  the 
nozzle  throat  also  proved  to  be  an  excellent  means  of 
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Fig.  11.  Mass  spectrometer  signal  for  C02  and 
CO  observed  for  a  3.6  km/s  beam  formed  by  laser 
breakdown  in  pure  C02.  The  CO  signal  has  been 
corrected  for  the  C02  cracking  fraction. 


i.e.,  with  greater  than  75%  H20  in  the  beam.  Several 
experiments  using  gas  mixtures  and  mixtures  of  liquid 
water  with  gases  were  also  performed  in  an  attempt  to 
generate  the  exact  mole  fractions  of  exhaust  species. 

The  run  matrix  is  summarized  in  Table  4.  The  results 
indicate  that  the  generation  of  an  exact  simulation  of  all 
the  species  at  the  correct  velocity  may  be  possible. 
However,  the  matrix  of  gases,  liquids,  and  setup 
parameters  are  so  broad  that  the  experiments  to  perfect 
such  a  simulation  were  beyond  the  scope  of  this 
preliminary  effort.  Interestingly,  we  found  out  that  the 
mixture  species  remained  largely  undissociated  in  the 
fast  beam.  This  observation  is  likely  due  to  the  large 
fraction  of  the  gas  pulse  that  is  not  processed  by  the 
laser  breakdown,  but  by  the  blast  wave.  This 
observation  further  reinforces  the  belief  that  the  whole 
exhaust  plume  can  ultimately  be  simulated.  The  beam 
velocity  in  these  multi-injection  experiments  was  easily 
tuned  but  we  observed  some  temporal  separation  of 
species  in  the  exhaust  pulse.  We  attribute  this  to  poor 
mixing  in  the  nozzle  and  believe  that  an  improved  design 
can  eliminate  this  problem. 

Simulations  of  exhaust  plumes  using  fast  molecular 
beams  of  a  single  or  pure  gas  are  more  valuable  to  our 
understanding  of  plume  impingement  effects  than  a 
simulation  of  the  exact  exhaust  composition.  The  exact 
beam  composition  will  determine  the  final  surface 
heating  and  momentum  accommodation  information  for 
the  experimental  conditions,  but  such  a  measurement  is 
too  complex  to  provide  insights  to  the  mechanisms  or 
key  components  involved.  On  the  other  hand,  the  single 
component  beams  provide  the  laboratory  tool  to  evaluate 
plume  impingement  heating,  momentum  transfer,  and 
even  surface  chemical  effects  under  controlled  conditions 
with  straightforward  interpretation.  Most  importantly, 
the  individual  species  beams  can  be  used  to  evaluate  the 
angular  dependent  momentum  accommodation 
coefficients. 

In  this  preliminary  work,  we  have  demonstrated  the 
ability  to  create  2  to  4  km/s  velocity  pulsed  beams  of  a 
variety  of  molecule  species.  These  beams  have  not  yet 
been  fully  calibrated  in  terms  of  both  flux  and  spatial 
dimensions,  however,  we  have  performed  a  "plume 
impingement  heating"  experiment  to  demonstrate  their 
utility.  Specifically,  we  have  experimentally  measured 
the  heating  due  to  the  impingement  of  3.6  km/s  carbon 
dioxide  and  nitrogen  plumes  on  a  thin  thermally  isolated 
nickel  foil  using  a  previously  developed  technique12  that 
is  illustrated  in  Fig.  12.  The  plume  gases  carry  kinetic 
energy  (and  possibly  chemical  and/or  thermal  enthalpy 
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Table  4.  Gas  Mixture  Experiments 


Gas  -  1% 

Gas  -2% 

Gas  -  3% 

Gas  -  4% 

Nozzle  Throat  Valve 

Side  Valve 

Nr50 

O2-50 

S> 

1 

8 

H2-62 

00 

Nr50 

O2-50 

Monomethy  lamine- 1 00 

H20  (Uq)-100 

H2-50 

N2-50 

H20  (Liq)-100 

H2-80 

N2-20 

Fast 

Beam 


Scattered 

Laser 

Radiation 


.—  Thermally 
Isolated 
Thin  Target 


Radiation 

_  ~r  _ 

Radiation 

o- 

—  1  — 

*  T0 

0<e<  1 
Sample 
Specific 

e  =  1 

Blackened 

Surface 

Surface 

Heating 

By  Energy 
Accom. 

— ►To  Optical 
Pyrometer 

Heating 

- 

Fig.  12.  Schematic  of  the  plume  impingement 
heating  experiment. 


but  these  are  negligible  in  comparison  to  the  kinetic 
energy  under  the  test  conditions).  Upon  striking  a 
surface  these  gases  will  deposit  some  fraction  of  this 
energy  into  the  surface  as  sensible  heat.  This  fractional 
deposition  is  defined  as  the  energy  accommodation 
coefficient.  For  normal  incidence,  without  chemical 
heating,  this  quantity  can  be  directly  related  to  the 
momentum  accommodation  coefficient.  The  purpose  of 
this  study  was  to  demonstrate  that  the  pulsed  beams  were 
sufficiently  intense  to  allow  employment  of  this 
diagnostic. 


A  thermally  isolated  foil  in  vacuum  can  only  lose 
heat  by  radiation.  It  must  ultimately  reach  a  steady  state 
temperature  that  reflects  a  balance  between  plume  energy 
accommodation  and  radiative  cooling.  Our  experimental 
procedure  involves  the  repetitive  exposure  of  the  front 
surface  of  a  nickel  foil  to  the  simulated  plume,  and 
recording  the  black  painted  back  surface  foil  temperature 
rise  due  to  the  plume  impingement  heating.  The 
temperature  was  meausred  with  an  Omega  Model  OS-602 


optical  pyrometer.  The  details  of  our  approach  and  the 
relationship  between  temperature  increase  and 
accommodation  coefficient  are  provided  in  Ref.  12. 

We  performed  heating  studies  on  a  nickel  foil  with 
both  N2  and  C02  beams.  A  0.6°  to  0.7°C  temperature 
rise  was  observed  in  the  foil  temperature  when  it  was 
exposed  to  a  3  Hz  3.6  km/s  beam  of  either  C02  or  N2 
as  shown  in  Fig.  13.  Scattered  10.6  micron  C02  laser 
radiation  from  the  nozzle  will  also  contribute  to  the  foil 
heating,  however  this  effect  is  small  and  was  measured 
to  be  0.2  C  by  placing  a  very  thin  polyethylene  film  in 
front  of  the  foil  which  allows  the  scattered  laser 
radiation  to  heat  the  foil  but  prevents  beam 
impingement.  The  signal  to  noise  of  the  data  is  clearly 
large  enough  to  provide  useful  data  on  plume 
impingement  heating  once  the  beams  are  calibrated. 


Fig.  13.  Typical  temperature  variation  on  nickel  foil 
observed  during  beam  impingement  experiment. 


These  individual  constituent  beams  should  provide  a 
powerful  tool  in  evaluating  plume  impingement  heating 
and  momentum  transfer  parameters  for  various 
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materials.  We  have  demonstrated  that  the  beams  are 
sufficiently  intense  to  produce  measurable  temperature 
changes  on  isolated  targets.  This  has  been  done  without 
an  optimized  system.  We  fully  anticipate  that  similar 
measurements  of  momentum  transfer  can  be  achieved  in 
our  system. 

Conclusions 

The  fast  O-atom  beam  source  has  proven  to  be  a  very 
versatile  facility  for  the  study  and  simulation  of  a  wide 
variety  of  phenomena  of  relevance  to  the  interaction  of 
the  ambient  atmosphere  with  spacecraft  (and  their  local 
environment)  operating  in  low  earth  orbit.  The 
development  and  phenomenological  understanding  of  thin 
film  silver  probes  and  the  implementation  of  a  time-of- 
flight  mass  spectrometer  within  the  test  facility  have 
significantly  enhanced  the  diagnostic  capabilities  and 
measurement  domain  of  the  system. 

The  recent  extension  of  the  laser  breakdown 
technique  to  form  lower  velocity,  2  to  4  km/s,  beams  of 
selected  molecular  species  and  mixtures  will  provide  the 
ability  to  examine  other  rarefied  flow  phenomena,  such  as 
exhaust  plume  impingement  effects,  in  a  small-scale 
laboratory  facility. 
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ABSTRACT 

We  have  used  Physical  Sciences  Inc/s  (PSD  pulsed 
8  km/s  oxygen  atom  source  to  perform  a  number  of 
measurements  of  material  interactions  of  interest  to 
low  earth  orbit  applications.  In  particular,  we  have 
developed  a  new  system  to  measure  the  scattering 
distribution  resulting  from  8  km/s  oxygen  atoms 
impinging  on  selected  spacecraft  materials.  We  will 
present  data  for  scattering  from  glass  and  black  and 
clear  anodized  aluminum.  We  will  also  present  results 
on:  the  fast  atom  oxidation  mechanism  for  thin  silver 
films;  the  O-atom  velocity  dependent  erosion  rates  for 
selected  materials;  and  the  oxidation  dependence  on 
material  temperature. 

Key  words:  oxygen  atoms,  LEO,  Ag,  material  oxidation 

1 .  INTRODUCTION 

In  the  mid-1  980's  we  developed  a  pulsed  fast  oxygen 
atom  beam  to  provide  a  ground  test  facility  to  qualify 
materials  to  be  used  in  low  earth  orbit  {LEO)  (Ref.  1 ). 

To  our  knowledge,  ours  is  the  only  system  that 
provides  a  high  flux  neutral  oxygen  atom  beam  with 
wide  area  operating  capability  at  the  desired  velocity  of 
8  km/s.  Sets  of  materials  (2.5  x  2.5  cm  samples)  have 
been  simultaneously  irradiated  to  8  km/s  oxygen  atom 
fluences  of  6  x  lO^1  cm'2  in  our  facility. 

Flight  tests  of  materials  are  preferable  but  considerably 
more  expensive  than  ground  evaluation.  Furthermore, 
it  is  difficult  to  acquire  material  aging  data  from  flight 
tests  which,  LDEF  notwithstanding,  are  generally 
limited  to  less  than  2  weeks  experiment  duration. 
Lastly,  the  flight  data  represent  an  average  over  the 
orbital  environmental  conditions  of  the  flight.  Many 
materials  exhibit  an  increased  rate  of  oxidative 
behavior  with  increased  temperature,  e.g.  Ref.  2. 
Materials  can  vary  in  temperatures  between  -40 °C  to 
+  65°C  as  the  surface  exposure  varies  from  deep 
space  to  solar  illumination  (pathological  materials 
excluded).  Thus  a  measurement  over  a  short  shuttle 
experiment  may  not  be  reliably  extended  to  a  longer 
mission  with  different  orbits,  even  if  the  aging  effects 
can  be  linearly  extrapolated. 

The  use  of  ground  test  data  has  been  problematical 
both  because  facility  definitions  have  been  vague  and 
because  standardized  measurement  techniques  have 
not  been  applied.  The  need  for  a  measurement 
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protocol  is  clear.  Banks  et  al.  have  recently  discussed 
this  for  asher  systems  (Ref.  3)  and  Minton  (Ref.  4)  has 
recently  presented  version  1  of  a  detailed  testing 
protocol. 

PSI  has  performed  commercial  evaluation  of  LEO 
materials  since  1 987  and  has  developed  a 
measurement  protocol  very  similar  to  that  outlined  by 
Minton.  The  first  step  in  such  a  protocol  is  the  clear 
definition  of  the  capabilities  and  operating  properties  of 
the  fast  oxygen  atom  facility.  Since  the  PSI  0-atom 
beam  source  has  been  used  for  material  evaluation  for 
the  past  7  years,  a  significant  effort  has  been 
dedicated  to  characterizing  its  operating  properties  in 
detail.  The  operating  principles  and  characteristics  of 
our  device  are  described  below.  This  discussion  is 
followed  by  some  representative  results  of  material 
behavior  under  simulated  LEO  conditions. 

1 .1  The  fast  atom  source 

The  hyperthermal  oxygen  atom  source  and 
experimental  configuration  have  been  described 
previously  (Refs.  1,5,6)  and  will  be  discussed  only 
briefly  here.  A  schematic  of  the  experimental  system 
is  shown  in  Figure  1 .  It  is  comprised  of  two  stainless 
steel  six-way  crosses,  including  an  8  in.  cross  source 
chamber  housing  a  pulsed  oxygen  valve/nozzle 
assembly,  connected  to  a  16  in.  cross  expansion 
chamber.  The  system  is  pumped  by  a  cryopump 
*  attached  to  the  large  cross  which  maintains  a  base 
pressure  of  3  x  10~  Torr. 


Figure  1 .  Schematic  of  the  fast- 2  system. 
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Geophysics  Directorate  of  the  USAF  Phillips  Laboratory  and  by  the  NASA  Lyndon  B.  Johnson  Space  Center, 
monitored  by  Mr.  Steve  Fitzgerald. 
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The  oxygen  atoms  are  generated  in  a  pulsed  laser 
discharge  of  pure  02  as  described  by  Caledonia  et  al. 
(Ref.  1).  A  12-J/pulse  C02  laser  is  focused  with  a 
1 00  cm  focal  length  BaF2  lens  into  the  throat  of  a 
20  deg  full  angle,  12.5  cm  long  conical  expansion 
nozzle  which  has  been  partially  filled  with  02  by  a 
pulsed  beam  valve.  The  plasma,  ignited  at  the  throat, 
expands  out  the  nozzle  dissociating  the  molecular 
oxygen  in  front  of  it.  The  nozzle  was  designed  to 
allow  ion  and  electron  recombination  while  the  slower 
kinetics  of  atom/atom  recombination  maintains  a  highly 
dissociated  beam.  The  result,  at  8  km/s,  is  an 
approximately  50  ps  pulse  of  highly  dissociated  oxygen 
(>  80%  atoms)  with  less  than  1%  ion  content  (the 
beam  is  charge  neutral).  The  general  beam  properties 
are  provided  in  Table  1  while  the  diagnostics  used  in 
evaluating  these  properties  are  summarized  in  Table  2. 


Table  1.  Fast  0-beam  properties. 


Velocity 

8  km/s  ±  15%  (5  to  12  km/s 
range) 

Fluence 

~1018  O-atoms/pulse,  3  Hz 

Composition 

>80%  oxygen  atoms 

Size 

Expandable  to  >  1 000  cm2  area 

Charge 

content 

<  1  %  ions  (controllable  by 
pseudo-Helmholtz  coils) 

Metastable 

content 

O^D)  concentration  <0.4% 

Temperature 

T  =  300  K 

VUV/UV 

content 

One  photon  per  1 04  0  atoms 
(similar  to  LEO) 

The  beam  velocity  (5  to  12  km/s)  is  selected  by 
varying  the  time  delay  between  the  pulsing  of  the  02 
valve  and  the  triggering  of  the  C02  laser.  The  delay 
determines  the  mass  of  02  processed  by  the  pulsed 
discharge  which  is  inversely  related  to  the  beam 
velocity.  The  velocity  is  determined  by  measuring  the 
time  of  flight  of  the  pulse  between  two  radiometers 
attached  to  the  source  chamber,  mounted  along  the 
beam  axis  and  separated  by  7.6  cm,  and  filtered  to 
monitor  the  5S-5P  oxygen  atom  line  (777  nm)  (see 
Figure  1). 

We  have  estimated  the  flux  of  oxygen  atoms  generated 
per  pulse  by  measuring  the  mass  of  02  processed  per 
pulse  and  the  expansion  characteristics  of  the  beam. 
The  expansion  properties  have  been  determined 
previously  (Refs.  5,6)  by  measuring  the  fast  O-atom 
erosive  mass  loss  of  polyethylene  samples  mounted 
75  cm  downstream  from  the  nozzle  throat  as  a 
function  of  radial  position.  The  expansion  half  angle 
(to  l/e  times  the  on-axis  flux)  was  found  to  be  1 4  deg. 
These  calibrations  have  been  cross  checked  by 
measuring  the  absolute  oxidation  rates  of  materials  and 
comparing  them  with  equivalent  space  flight 
measurements.  More  recently,  as  discussed  in  the 
next  section,  the  oxidation  rate  of  thin  silver  films  has 
been  used  as  a  calibration  standard. 

Although  the  beam  is  charge  neutral,  it  does  have  a 
small  but  finite  ion  content.  The  charge  concentration 
can  be  decreased  by  up  to  a  factor  of  eight  using  a  pair 
of  pseudo-Helmholtz  coils.  Thus  any  effects  due  to 
charged  rather  than  neutral  species  can  be  delineated. 
Control  of  the  ionic  concentration  is  critical  in  that 
some  materials  are  efficiently  oxidized  by  O  +  and  not 
by  O.  The  local  charge  concentration  can  be  measured 
through  use  of  a  retractable  Langmuir  probe.  The 
probe  is  positionable  both  horizontally  and  vertically 
and  it  has  been  found  that  the  spatial  and  temporal 
shape  of  the  ion  distribution  matches  that  of  the 
neutral  distribution  (Ref.  6). 


Table  2.  O-Beam  diagnostics 


Velocity,  temporal  pulse  shape 

Radiometers,  mass  spectrometer  (standard  quadrupole  and 

TOF) 

Pulse  fluence 

Mass  flow  meter,  thin  silver  film  actinometer 

Beam  composition 

Mass  spectrometer 

Excited  state  concentrations 

Spectrometer,  bandpass  filters 

Charge  content  (beam  is  charge  neutral) 

Langmuir  probe 

Beam  shape 

Erosive  witness  samples,  positionable  Langmuir  probe,  2D 
visible  imaging 

Beam  translational  temperature 

Slitted  thin  silver  film  target  chamber 

We  have  recently  measured  the  O^D)  concentration  of 
the  beam  bv  monitoring  the  emission  from  the 
forbidden  0(  D  3P)  transition  using  a  specially 
constructed  bank  of  narrow-band  interference  filters. 
The  O^D)  concentration  was  found  to  be  <  0.4%  for 
a  beam  velocity  of  8  km/s  and  to  increase 
monotonicaily  with  beam  velocity  over  the  velocity 
range  of  6  to  12  km/s. 


We  have  also  evaluated  the  translational  temperature 
of  the  beam.  This  was  done  by  passing  the  O-beam 
through  a  0.15  cm  diameter  circular  aperture  and 
measuring  the  oxidation  pattern  on  a  thin  silver  film 
witness  plate  placed  5  cm  behind  the  aperture.  Any 
growth  in  the  beam  beyond  the  slit  width  is  due  to  the 
Maxwellian  motion  of  the  0  atoms  and  can  be  related 
to  an  effective  temperature  for  the  8  km/s  velocity 
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beam.  This  temperature  was  deduced  to  be  -  300  K 
at  a  distance  of  85  cm  from  the  nozzle  throat.  The 
beam  temperature  will  be  higher  closer  to  the  throat 
and  is  calculated  to  be  -  1 000  K  at  the  nozzle  exit 
plane. 

Lastly,  we  have  measured  the  VUV/UV  content  of  the 
beam  and  find  that  it  is  somewhat  less  then  the 
equivalent  LEO  value  per  oxygen  atom.  Since  the 
beam  is  pulsed,  the  VUV  radiation  can  be  chopped 
allowing  independent  examination  of  VUV  effects  on 
materials. 

Our  O-beam  facility  has  been  employed  on  several 
occasions  to  pre-test  materials  which  were 
subsequently  flown  in  space  experiments.  In  general, 
the  comparison  between  flight  data  and  ground  test 
data  has  been  quite  good.  The  most  recent  large 
ground  test  was  for  the  SEE  program  on  materials  for 
the  EOIM-3  flight  (Ref.  7}.  We  irradiated  84  materials 
to  a  fluence  of  2  x  1020  cm'2  8  km/s  oxygen  atoms 
and  followed  a  careful  protocol  of  pre-  and  post-test 
weighing  taking  account  of  water  uptake.  The  Jet 
Propulsion  Laboratory  performed  subsequent  surface 
analysis.  With  a  few  understood  exceptions,  the 
comparison  between  our  ground  test  results  for  these 
materials  and  the  flight  observations  was  quite  good 
(Ref.  7),  demonstrating  our  ability  to  simulate  LEO 
conditions. 

2.  EXAMPLE  MATERIAL  STUDIES 

Material  oxidation  rates  can  be  significantly  affected 
both  by  the  test  environment  and  the  test  system 
parameters.  In  this  section  we  will  briefly  examine  the 
effect  of  sample  temperature  and  O-atom  velocity  on 
two  materials,  a  thin  silver  film  and  //-crystalline 
carbon.  In  the  case  of  the  thin  silver  film  the  effect  of 
two  additional  parameters,  test  duration  and  surface 
contamination  will  be  addressed. 

We  have  performed  investigations  of  the  temperature 
dependent  oxidative  behavior  of  several  materials  over 
the  temperature  range  of  0°  to  85 °C.  An  example  of 
such  studies  is  our  observations  of  the  oxidative 
behavior  of  thin  films  of  silver  (600A)  (Ref.  8)  and 
//-crystalline  carbon  (1  //m)  which  were  vacuum 
deposited  onto  the  crystals  of  quartz  crystal  monitors 
(QCM's).  QCM's  are  the  diagnostic  of  choice  for 
monitoring  mass  change  as  the  measurements  are 
performed  continuously  in-situ.  Unfortunately  in  many 
instances  the  materials  to  be  tested  cannot  be  coated 
on  a  crystal.  Then  weight  change  is  monitored  by  pre- 
and  post-testing  weighing  using  our  Mettler  analytical 
balance  and  Cahn  micro-balance.  (Here  a  careful 
protocol  must  be  followed  to  account  for  moisture 
uptake). 

The  silver  and  carbon  films  are  good  examples  for 
consideration  in  that  the  silver  gains  mass  upon 
oxidation  while  the  carbon  loses  mass.  In  these 
experiments,  the  QCM  was  mounted  directly  onto  a 
temperature-controlled  aluminum  block  and  thus  the 
quartz  crystal,  which  is  in  intimate  thermal  contact 
with  the  QCM  base,  was  maintained  to  ±0.5°C 
(material  coupons  can  be  monitored  in  the  same 
manner).  Temperature  regulation  of  the  QCM  was 
necessary  since  the  oscillator  frequency  is  temperature 


dependent.  Therefore,  a  change  in  temperature  could 
be  erroneously  interpreted  as  a  mass  gain  or  loss. 

The  silver  films  are  particularly  interesting  since 
because  of  their  high  oxidation  efficiency  we  feel  that 
they  can  be  used  as  a  standard  for  laboratory  flux 
calibrations.  Furthermore,  they  have  been  successfully 
used  in  both  flight  (Ref.  9)  and  laboratory  (Ref.  6) 
studies  to  examine  the  distributions  or  hyperthermal 
oxygen  atoms  scattered  off  of  selected  surfaces.  The 
oxidation  behavior  of  these  films  exhibits  several 
intriguing  effects.  Our  studies  in  this  are  described  in 
detail  in  Ref.  8  and  presented  briefly  below. 

Typical  thin  film  silver  oxidation  histories  for  irradiation 
by  8  km/s  oxygen  atoms  at  room  temperature  are 
shown  in  Figure  2.  Shown  is  oxidation  depth  versus 
the  number  of  impinging  0-atom  pulses  as  determined 
by  the  measured  mass  increase  for  two  film 
thicknesses.  Two  important  features  can  be  seen. 
First,  no  oxidation  occurs  during  the  first  two  hundred 
pulses.  This  is  typically  observed  and  we  relate  it  to 
the  required  removal  of  a  contamination  layer  on  the 
silver  film.  Indeed,  if  we  first  irradiate  the  film  with 
several  hundred  pulses  of  a  fast  argon  beam  and  then 
shortly  thereafter  impinge  it  with  fast  oxygen  atoms, 
we  find  that  oxidation  initiation  is  immediate.  Second, 
after  oxidation  initiation,  the  oxidation  depth  increases 
linearly  with  O-atom  fluence  for  approximately  the  first 
250A  of  the  silver  film  and  then  the  oxidation  rate 
slows  down.  We  believe  the  oxidation  rate  becomes 
diffusion  controlled  at  these  thicknesses  and  thus  only 
very  thin  films  are  useful  as  a  linear  calibration  source. 
Note  that  the  observed  oxidation  rate  continues  to 
decrease  very  dramatically  in  the  case  of  the  1  //m 
thick  film. 


Figure  2.  Oxidation  behavior  of  thin  silver  films 
coated  on  QCM  crystals  ( O-atom  velocity  = 

8  km/s ,  room  temperature}. 


Note  that  in  the  interpretation  of  the  QCM  data  it  is 
important  to  understand  the  oxidation  process.  During 
the  course  of  our  work  we  realized  that  the  poly¬ 
crystalline  silver  oxidized  to  the  peroxide,  AgO,  rather 
than  Ag20.  This  was  determined  by  measuring  the 
mass  uptake  of  a  totally  oxidized  film  with  known 
silver  loading,  and  verified  both  by  heating  the  film  to 
200°C  to  form  Ag20  and  by  X-ray  diffraction  analysis 
of  the  film. 
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These  limited  tests  on  silver  films  demonstrate  a 
number  of  interesting  features.  From  the  perspective 


of  a  ground  test  facility,  the  two  that  are  most 
important  are  the  demonstrated  role  of  surface 
contamination  and  the  observation  of  the  time 
(thickness)  variation  in  the  oxidative  behavior. 

In  further  studies  we  examined  the  temperature 
dependence  of  the  silver  and  carbon  oxidation  rate 
over  the  range  of  0°  to  85 °C.  These  results  are 
shown  in  Figure  3.  (For  convenience,  all  results  are 
normalized  to  a  silver  reactivity  of  unity  at  273  K). 
Although  there  is  some  scatter  in  the  data,  the  silver 
oxidation  rate  is  apparently  independent  of 
temperature,  unlike  that  for  carbon.  Thus,  we  believe 
that  these  silver  films  will  provide  a  good  standard  for 
relative  investigation  of  high  velocity  oxidation 
behavior  of  materials.  On  the  other  hand,  the 
oxidation  rate  for  this  particular  form  of  carbon  varies 
significantly  with  temperature,  a  critical  effect  for 
predicting  lifetime  if  this  material  were  to  be  used  in  an 
engineering  design. 


Figure  3.  Measured  silver  and  carbon  reactivity 
versus  temperature  (O-atom  velocity  8  km/s ). 


We  have  also  investigated  the  effect  of  O-atom 
velocity  on  the  oxidation  rate.  This  study  is 
complicated  by  the  fact  that  the  number  of  oxygen 
atoms  per  pulse  in  our  beam  decreases  as  the  beam 
velocity  increases.  Rather  than  go  into  the  details  of 
our  velocity  dependent  calibration  we  will  only  discuss 
the  relative  oxidation  behavior  of  the  thin  silver  film 
and  carbon  samples  as  the  beam  velocity  is  varied. 

The  carbon  and  silver  film  samples  were  irradiated 
simultaneously  at  velocities  of  8,  10  and  1 2  km/s  at  a 
pulse  frequency  of  2  Hz.  The  observed  oxidation  depth 
versus  time  of  the  thin  silver  film  is  shown  in  Figure  4. 


Figure  4.  Observed  silver  oxidation  depth 
versus  time  for  three  O-atom  velocities. 
Pulse  repetition  frequency  was  2  Hz. 


The  decreasing  oxidation  rate  with  increasing  velocity 
correlates  with  the  decreased  number  of  oxygen  atoms 
produced  per  pulse  as  the  velocity  is  increased.  Indeed 
within  the  uncertainty  of  our  calibration  the  silver  film 
oxidation  efficiency  remained  constant  over  the 
velocity  range.  (This  is  not  surprising  since  the 
oxidation  efficiency  is  near  unity.) 

This  result  is  contrasted  with  the  observed  etched 
depth  histories  in  the  carbon  sample  as  shown  in 
Figure  5.  Here  it  can  be  seen  that  the  etch  rate 
increased  with  velocity  even  though  the  O-atom 
fluence  decreased.  Thus  the  carbon  sample  studied 
shows  a  strongly  velocity  dependent  oxidation  rate. 
We  have  observed  a  similar  trend  over  the  velocity 
range  of  6  to  11  km/s  in  other  materials,  such  as 
kapton  (Ref.  10). 


Figure  5.  Observed  etch  depth  of  carbon  versus 
time  for  three  O-atom  velocities.  Pulse 
repetition  frequency  was  2  Hz. 


3.  SCATTERING  STUDIES 

We  employed  the  oxidative  behavior  of  thin  silver  films 
to  measure  the  scattering  distribution  of  8  km/s 
oxygen  atoms  impinging  selected  surfaces.  These 
scattering  distributions  are  of  interest  to  the  materials 
community  since  non-ram  surfaces  can  be  oxidized  by 
scattered  ambient  oxygen  atoms. 

Silver  film  scatterometers  designed  to  measure  the 
atomic  oxygen  scattering  distribution  off  various 
materials  have  been  flown  on-orbit  in  shuttle 
experiments  by  John  Gregory  (Ref.  9)  and  other 
investigators.  These  on-orbit  measurements  show  that 
the  atomic  oxygen  scattering  is  largely  non-specular. 
Our  measurements  attempt  to  confirm  these 
observations  in  the  laboratory  in  an  8  km/s  atomic 
oxygen  flow.  The  representative  materials  chosen  for 
this  study  are  clear  and  black  anodized  aluminums  and 
glass.  As  described  previously  (Ref.  6),  we  observed 
that  scattering  off  all  three  materials  is  largely  non- 
specular.  T*-e  details  of  the  scatterometer  design 
atomic  oxygen  scattering  measurements  are  described 
briefly  below. 

We  used  a  modified  version  of  Gregory's  scatterometer 
in  our  laboratory  experiments  to  measure  the  oxygen 
atom  angular  scattering  distributions  from  surfaces.  In 
this  device  the  beam  passes  through  a  small  slit  into  a 
chamber  containing  the  target  material.  The  O-atoms 
scattered  from  the  target  ultimately  interact  with  a 
silver  film  "witness  plate"  positioned  at  the  front  of  the 
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chamber  (around  the  slit).  The  distance  between 
target  and  slit  must  be  large  enough  to  allow  sufficient 
spread  of  the  scattered  beam  to  provide  good 
measurement  resolution  on  the  witness  plate  for  the 
scattered  angular  distribution. 

We  designed  and  built  silver  film  scatterometers 
specific  to  our  experiment.  Two  scatterometers  were 
constructed  to  hold  up  to  three  samples  in  individual 
compartments  at  either  a  90  or  45  deg  angle  of 
incidence  to  the  8  km/s  atomic  oxygen  flow.  The 
90  deg  scatterometer  is  shown  in  Figure  6.  The  silver 
film  to  be  oxidized  by  the  scattered  atomic  oxygen  was 
attached  to  the  2.22  cm  radius  surface.  A  0.25  cm 
wide  by  1 .25  cm  long  slit  was  cut  in  the  film  prior  to 
mounting  to  restrict  the  atomic  oxygen  flow  into  the 
scatterometer.  The  glass  and  aluminum  samples  were 
cut  to  fit  the  individual  sample  chambers  and  attached 
to  the  silver  coated  kapton  covered  rear  surface.  A 
key  feature  of  our  scatterometers  not  found  in  other 
designs  is  the  addition  of  a  hinged  cover  to  permit 
cleanup  via  argon  atom  bombardment  of  both  the 
sample  and  silver  film  surfaces  to  remove  the  surface 
contamination  layers  in  the  chamber  immediately  prior 
to  the  oxygen  scattering  measurement.  We  have 
found  that  this  was  a  key  feature  for  successful 
accomplishment  of  the  experiments. 


a  factor  of  two  but  this  uncertainty  can  be  greatly 
reduced  in  future  work  by  simultaneously  monitoring 
the  O-atom  flux  during  the  scattering  measurement. 


Figure  7.  Oxygen  scattering  distribution  at  norma l 
incidence  for  three  materials .  O-atom  velocity  is 
8  km/s.  Low  measurements  at  positive  angles  on 
black  anodized  aluminum  are  due  to  improper 
witness  film  cleaning. 

The  scattering  distribution  is  quite  flat,  suggesting 
largely  isotropic  scattering  for  all  three  materials. 
Similar  results  were  observed  at  45  deg  incidence  with 
no  strong  specular  regions  apparent.  The  central  spike 
that  is  seen  in  the  figure  should  be  ignored.  It  is  due 
to  the  presence  of  a  slit  in  the  silver  film  and  the 
fluence  indicated  is  an  artifact  of  the  analysis 
technique. 

The  importance  of  working  with  clean  films  was 
accidentally  demonstrated  in  the  data  for  black 
anodized  aluminum  shown  in  Figure  7.  It  turns  out 
that  the  door  of  the  scatterometer  opened  too  far 
during  the  exposure  of  the  black  anodized  sample  and 
one  side  of  the  silver  was  not  scoured  by  the  atom 
beam.  Thus,  the  anomolously  low  apparent  scattered 
fluence. 


Figure  6.  Normal  incidence  spectrometer  used  to 
measure  distribution  of  fast  oxygen  atoms 
scattered  off  of  selected  materials. 

Our  witness  samples,  500A  silver  film  on  2  mil  mylar, 
were  prepared  by  Sheldahl.  In  this  application,  the 
films  were  analyzed  for  oxidation  by  measuring 
changes  in  film  transmission  (silver  oxide  is  transparent 
although  silver  is  opaque).  The  fast  O-atom  beam  was 
used  to  calibrate  the  film  transparency,  e.g.,  optical 
thickness  versus  atomic  oxygen  fluence,  and  the 
scattering  measurement  experiment  duration  was 

limited  to  ensure  that  <  250A  depth  of  the  film  was 
oxidized,  i.e.,  to  the  range  of  linear  oxidation  behavior. 

The  resulting  measured  scattered  oxygen  atom 
distributions  versus  angle  from  the  slit  is  shown  in 
Figure  7  for  normal  incidence.  All  measurements  were 
taken  in  the  central  chamber.  The  results  are  averaged 
over  the  width  of  the  witness  film  and  are  presented  as 
a  percentage  of  the  incident  flux  per  degree.  The 
angular  resolution  of  the  data  is  about  2.5  deg.  The 
absolute  calibration  for  this  data  is  uncertain  to  about 


Subsequent  measurements  made  at  a  few  scattering 
angles  with  a  time-of-flight  mass  spectrometer  showed 
that,  although  the  velocity  spread  of  the  scattered 
atoms  was  quite  broad,  a  significant  degree  of 
thermalization  occurred.  The  details  of  these  measure¬ 
ments  are  presented  in  Ref.  6. 

4.  SUMMARY 

Ground  testing  provides  a  critical  adjunct  to  the 
evaluation  and  validation  of  material  applicability  for 
low  earth  orbit  operations.  Care  must  be  taken, 
however,  to  ensure  that  the  ground  test  facility  can 
provide  valid  test  results.  In  the  case  of  hypersonic 
oxygen  atom  testing,  proper  material  evaluation  and 
data  interpretation  require  both  a  clear  understanding 
of  the  characteristics  of  the  test  facility  as  well  as  a 
standardized  measurement  protocol. 

Details  of  the  RSI  fast  oxygen  atom  source  have  been 
provided.  This  system  has  been  used  to  evaluate 
materials  since  1  987,  over  2000  material  samples  have 
been  tested.  Good  comparison  with  space  test  results 
has  been  demonstrated.  All  major  beam  parameters, 
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e.g.,  VUV,  metastable  content,  beam  velocity  spread, 
have  been  evaluated  and  a  standard  measurement 
protocol  has  been  applied. 

Example  test  data  has  been  presented  for  8  km/s 
oxygen  atom  beam  irradiation  of  a  thin  silver  film  and 
a  sample  of  micro-crystalline  carbon.  Typical 
measurements  of  the  scattering  distributions  of  8  km/s 
oxygen  atoms  impinging  representative  spacecraft 
surfaces  were  also  shown.  Although  these  data  are  of 
interest  in  their  own  right,  they  also  demonstrate  the 
importance  of  controlling  test  parameters  such  as 
contamination  level,  sample  temperature,  and  beam 
velocity. 
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ABSTRACT 

We  have  measured  the  Einstein  coefficient  for  the  forbidden  atomic  transition,  N(2P-4S), 
at  346.6  nm,  by  determining  absolute  photon  emission  rates  at  346.6  nm  when  known  number 
densities  of  N(2P)  are  in  the  detector  field  of  view.  The  N(2P)  is  generated  in  a  discharge-flow 
reactor  and  its  number  density  is  determined  by  vacuum  ultraviolet  resonance  absorption  at 
174  nm.  The  absolute  response  of  the  monochromator  is  determined  from  measurements  of  the 
O/NO  air  afterglow.  This  first  experimental  determination,  (5.4  ±1 .4)  x  10'3  s‘\  agrees  well  with 
published  theoretical  calculations. 

Key  words:  Einstein  coefficient,  metastable  nitrogen  atoms,  aurora 
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1.  INTRODUCTION 


Metastable  N(2P)  nitrogen  atoms  are  a  source  of  bright  auroral  radiation  at  both  346.6  nm 
in  the  ultraviolet,  from  the  N(2P-4S)  transition,  and  at  1040  nm  in  the  near  infrared,  from  the 
N(2P-2D)  transition.1'3  Observations  of  N(2P)  emission  intensities  in  comparison  with  intensities 
of  N2+(B-X)  or  N2(C-B)  transitions  have  been  used  to  estimate  the  depth  of  penetration  of  auroral 
electrons2  and  the  altitude  profiles  of  atomic  oxygen  number  densities.3  Both  of  these  measure¬ 
ments  rely  on  the  fact  that  long-lived  N(2P)  is  quenched  rather  efficiently  by  atomic  oxygen  under 
auroral  conditions,  whereas  the  molecular  emissions  are  not  because  of  their  short  radiative 
lifetimes. 

Proper  analysis  of  N(2P)  emission  measurements  in  aurorae  require  accurate  knowledge 
of  the  N(2P  -  2D,  4S)  transition  probabilities  and  the  rate  coefficient  for  N(2P)  quenching  by 
atomic  oxygen.  We  determined  the  latter  experimentally  several  years  ago 4  The  former, 
however,  is  known  only  through  two  theoretical  calculations  published  some  30  years  apart.5-6 
Although  the  two  calculations  agree  within  12%,  our  observations  of  N(2P)  excitation  in  the 
energy  transfer  reaction  between  N2(A)  and  N(4S)7  suggested  the  calculation  of  the  N(2P-4S) 
transition  probability  could  possibly  be  erroneous,  and  an  experimental  determination  of  that 
transition  probability  was  warranted. 


2.  EXPERIMENTAL 

The  intensity  of  an  electronic  transition  is  given  by  the  product  of  the  number  density  of 
species  in  the  emitting  state  and  the  Einstein  coefficient  for  spontaneous  radiation: 

I  =  NA,  (1) 

where  I  is  the  absolute  photon  emission  rate  in  units  of  photons  cm'3  s'1,  N  is  the  number  density 
in  species  cm'3,  and  A  is  the  Einstein  coefficient  in  units  of  s'1.  The  basis  of  our  measurement, 
therefore,  is  to  determine  absolute  photon  emission  rates  and  number  densities  of  emitting 
species  simultaneously.  The  ratio  of  these  two  quantities,  then,  is  the  Einstein  coefficient  we 
seek. 


Our  experiments  used  the  5  cm  diameter  discharge-flow  reactor  shown  schematically  in 
Figure  1 .  N(2P)  was  generated  in  the  upstream  section  of  the  reactor  in  a  microwave  discharge 
sustained  in  a  flow  of  N2  dilute  in  either  helium  or  argon  at  total  pressures  between  0.3  and 
1 .5  Torr.  The  N(2P)  was  then  detected  in  emission  at  346.6  nm  using  a  0.5  m  monochromator 
equipped  for  photoelectric  detection  and  also  by  vacuum  ultraviolet  resonance  absorption  at 
174  nm.  The  absorption  measurements  are  made  across  the  central  diameter  of  the  flow  reactor. 
The  source  of  the  resonance  radiation  was  a  microwave  excited  discharge  in  helium  with  a  trace 
of  added  air.  On  the  opposite  side  of  the  flow  reactor,  a  0.2  m  vacuum  monochromator  and 
solar-blind  photomultiplier  detected  the  radiation.  Normal  to  the  diameter  defined  by  the 
resonance  absorption  system,  a  0.6  cm  diameter  fiber-optic  bundle  collected  346.6  nm  radiation 
emanating  from  the  center  of  the  flow  reactor  and  transported  it  to  the  0.5  m  monochromator. 
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The  number  density  of  the  N(2P)  atoms  varies  radially  across  the  flow  reactor,  being 
highest  in  the  center  and  zero  at  the  walls.  The  absorption  and  emission  diagnostics,  therefore, 
both  sample  average  number  densities  in  the  reactor.  Because  the  system  and  flows  are 
cylindrically  symmetric,  however,  and  because  the  two  diagnostics  have  similar  fields  of  view 
across  the  flow  reactor  (narrow  and  collimated),  the  average  number  densities  they  monitor  are 
the  same.  Thus,  the  non-uniform  N(2P)  number  densities  have  no  effect  on  the  experimental 
results. 


The  0.5  m  monochromator  was  calibrated  for  absolute  photon  emission  rate  measure¬ 
ments  by  observing  the  intensity  of  the  O/NO  air  afterglow  when  known  amounts  of  O  and  NO 
were  added  to  the  reactor.  Titrating  N  atoms  with  NO,  using  procedures  we  have  detailed 
previously,8,9  produced  the  O.  Our  air  afterglow  calibrations,  made  at  580  nm,  used  a  chemi¬ 
luminescence  rate  coefficient  of  1.17  x  10'19  cm3  molecule'1  s'1  nm'1.  This  value  is  an  average 
based  upon  our  critical  analysis10  of  the  published  literature1118  on  the  air  afterglow.  We  think 
this  value  is  accurate  within  ±10%. 

We  established  the  relative  spectral  response  of  the  system  in  situ  to  extrapolate  the 
calibration  made  at  580  nm  to  346.6  nm.  We  used  a  branching  ratio  technique  based  upon 
measurements  of  NO(B  2II  —  X  2II)  band  intensities  and  the  NO(B-X)  transition  probabilities  we 
reported19  (which  were  confirmed  theoretically  by  Langhoff,  et  al.20  and  experimentally  by  Luque 
and  Crosley21)  to  determine  the  response  variation  between  250  and  500  nm,  and  O/NO  air 
afterglow  intensity  measurements  to  establish  the  relative  spectral  response  between  400  and 
850.  The  two  techniques  agreed  to  within  ±5%  in  the  region  of  spectral  overlap. 

We  have  discussed  our  use  of  VUV  resonance  absorption  techniques  to  determine  atomic 
number  densities  in  a  discharge-flow  reactor  previously.22,23  The  174.3  nm  transition  that  is 
absorbed  by  N(2P)  is  actually  a  multiplet  consisting  of  four  separate  transitions,  two  doublets 
separated  by  about  0.25  nm  with  a  separation  between  the  lines  in  each  doublet  of  about 
0.0008  nm  or  0.26  cm'1.  The  Doppler  width  of  each  line  is  about  0. 17  cm'1,  small  enough  that 
the  unresolved  multiplet  can  be  treated  as  a  composite  of  four  independent  transitions. 

Our  measurements  were  made  with  a  resolution  of  1  nm  to  ensure  that  all  lines  of  the 
multiplet  were  transmitted  equally  through  the  instrument.  The  equation  describing  the  fractional 
absorption,  A=(I0-I)/Io,  then  becomes24 


A  = 


(2) 


where  Q  is  the  relative  intensity  of  each  multiplet  component,  koi  the  absorption  coefficient  at 
line  center,  oo„  the  reduced  frequency  of  the  line,  a  is  the  ratio  of  the  width  of  the  line  from  the 
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resonance  lamp  to  that  of  the  absorption  line,  and  f  is  the  absorption  path  length.  The  reduced 
frequency,  co,  is  given  in  terms  of  the  separation  from  the  center-line  frequency,  v0,  by 


O) 


2(v-vo) 


(3) 


where  AvD  is  the  full  width  at  half  maximum  of  the  line  as  determined  by  Doppler  broadening: 


AVS'‘)  =  -f  ■ 


2RT  ln2 
M 


(4) 


X0  being  the  center  wavelength  of  the  transition,  R  the  gas-law  constant,  T  the  translational 
temperature,  and  M  the  mass  of  the  atom. 

Doppler  broadening  is  the  only  significant  mechanism  of  line  broadening  in  a  low  power, 
low  pressure  microwave  discharge  resonance  lamp  and  of  the  absorption  line  as  well  (collision 
broadening  contributes  less  than  1%  to  the  total  linewidth  in  3  Torr  of  helium).25  As  a  result  the 
ratio  of  the  widths  of  the  lamp  line  to  the  absorption  line  is  just  the  square  root  of  the  ratio  of 
translational  temperatures  in  the  lamp  and  the  flow  reactor.  The  translational  temperature  in  the 
flow  reactor  is  equal  to  the  laboratory  temperature,  300  K.  We  determined  the  translational 
temperature  in  the  resonance  lamp  to  be  500  K  from  measurements  of  the  band  profiles  of 
N2(C  3IIu  -  B  3Hg)  transitions.  Figure  2  shows  a  comparison  between  a  fit  of  the  observed 
N2(C-B)  bands  and  band  contours  calculated  assuming  rotational  temperatures  of  400,  500,  and 
600  K.  The  rotational  profile  calculated  from  a  500  K  temperature  best  reproduces  the  data. 

The  uncertainty  in  this  value,  ±50  K,  contributes  to  an  uncertainty  in  the  overall  absorption 
measurement  of  only  ±3%. 

The  absorption  coefficient  at  line  center  is 


k  =  — 

°  Av„N 


ln2  Vdu 


n  8  k  d 


AulN 


(5) 


where  du  and  J,  are  the  degeneracies  of  the  upper  and  lower  states  in  each  of  the  four  transitions 
in  the  multiplet,  A  the  Einstein  coefficient,  and  N  the  number  density  of  absorbing  atoms.  We 
used  transition  probabilities  based  upon  the  radiative  lifetime  for  the  upper  state  of  the  multiplet 
of  1.8  ns  as  measured  by  Hutchison.26  Table  1  summarizes  the  relevant  spectroscopic  data. 


3.  RESULTS 


Figure  3  shows  typical  results  of  measurements  of  the  346.6  nm  intensity  as  a  function  of 
N(2P)  number  density.  The  N(2P)  number  densities  were  determined  by  measuring  the  fractional 
absorption  at  174  nm  and  using  the  results  of  a  numerical  integration  of  Eq.  (2)  to  relate  the 
fractional  absorption  to  N(2P)  number  density.  Since  the  absolute  photon  emission  rate  is  the 
product  of  the  number  density  of  emitters  and  the  Einstein  coefficient  for  the  transition,  the 
Einstein  coefficient  for  the  346.6  nm  transition  is  just  the  slope  of  the  line  through  the  data  in 
Figure  3. 

The  average  result  from  several  measurements  in  both  Ar  and  He  carrier  gas  is  5.4  x 
10'3  s'1.  We  estimate  the  uncertainty  in  this  value  to  be  ±26%  based  upon  the  square  root  of  the 
sum  of  the  squares  of  the  various  uncertainties  in  the  experimental  parameters:  5%  in  the  raw 
data  (ratio  of  intensity  to  number  density),  12%  in  the  number  density  diagnostic  (lamp 
temperature  and  transition  probabilities),  10%  in  the  air  afterglow  rate  coefficient,  and  20%  in 
the  photometric  calibration  of  the  monochromator. 

4.  DISCUSSION 

This  Einstein  coefficient  was  first  calculated  by  Garstang  to  be  5.4  x  10'3  s'1  over  40  years 
ago.5  More  recently,  Godefroid  and  Froese-Fischer6  calculated  a  value  of  4.8  x  10'3  s'1.  Both  of 
these  values  agree  well  with  our  experimental  determination.  We  had  speculated  these  calculated 
values  might  be  erroneous  because  we  determined  a  rate  coefficient  for  N(2P)  formation  in  the 
energy-transfer  reaction  between  N2(A)  and  N(4S)  to  be  five  times  larger  than  the  rate  coefficient 
for  N2(A)  quenching  by  N(4S).7  The  most  likely  explanation  for  this  discrepancy  seemed  to  be 
that  the  Einstein  coefficient  used  to  determine  either  N2(A)  or  N(2P)  number  density  was  in  error. 
We  have  considered  the  issue  of  the  N2(A)  Einstein  coefficient,27  and  determined  that  previously 
accepted  values  are  not  seriously  in  error.28  The  results  of  this  investigation  indicate  that  the 
same  can  be  said  for  the  N(2P-4S)  transition.  Thus  the  discrepancy  in  our  previous  measurements 
must  result  from  an  additional  metastable  in  our  flow  reactor  that  we  have  so  far  been  unable  to 
detect.  Evidence  for  this  additional  metastable  has  been  found  in  other  studies  as  well.29, 30 
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Figure  Captions 


Figure  1. 
Figure  2. 

Figure  3. 


Schematic  view  of  discharge-flow  reactor  used  in  the  experiments. 

Comparison  of  the  spectrum  of  N2(C  3IIU  -  B  3IIg)  emission  from  the  lamp  used  in 
the  absorption  measurements  with  synthetic  spectra  generated  assuming  different 
rotational  temperatures  for  the  N2(C)  emission.  The  best  fit  rotational  temperature, 
500  K,  corresponds  to  the  kinetic  temperature  in  the  lamp. 

Typical  experimental  results  showing  the  relationship  between  observed  N(2P  -  4S) 
emission  intensities  and  measured  N(2P)  number  densities. 
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Table  1.  Spectroscopic  Data  for  the  N(2p3  2P°  -  2p2  3s  2P)  Multiplet  at  174  nm 


Wavelength 

(nm) 

■ 

■ 

AU| 

(s') 

■ 

^o/^total 

(cm2) 

1 

174.272 

2 

a 

0.247  x  108 

1 

0.984  x  1013 

0.328  x  1013 

2 

174.273 

4 

D 

1.233  x  108 

5 

2.455  x  10-13 

1.637  x  1013 

3 

174.525 

2 

2 

0.987  x  108 

2 

1.974  x  1013 

0.658  x  10  '3 

4 

174.526 

n 

2 

0.493  x  108 

1 

0.493  x  10  ' 3 

0.329  x  10'3 
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Observations  on  the  O/NO  Air  Afterglow 

Lawrence  G.  Piper’ 
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Andover,  MA,  01810-1077 

ABSTRACT 

We  have  reviewed  the  literature  on  the  quantitative  aspects  of  the  O/NO  air  afterglow  to 
reconcile  conflicts  and  confusions  in  the  literature.  Most  previous  reports  of  the  relative 
distribution  of  intensity  as  a  function  of  wavelength  agree  reasonably  well  for  wavelengths 
shorter  than  700  nra.  At  longer  wavelengths,  however,  consensus  on  the  distribution  is  lacking. 
Consequently,  we  report  new  measurements  of  this  distribution  between  400  and  1750  nm.  The 
absolute  values  of  the  O/NO  emission  rate  coefficient  have  been  reported  in  a  variety  formats, 
which  renders  comparison  of  the  different  measurements  difficult.  We  have  reanalyzed  the 
previous  reports  and  find  a  consensus  value  to  be  (1.17  ±  0.1 1)  x  10'19  cm3  molecule'1  s"1  nm'1  at  a 
wavelength  of  580  nm  and  room  temperature  (300  K).  The  total  rate  coefficient  integrated  over 
all  wavelengths  out  to  3500  nm  is  8.9  x  10'17  cm3  molecule'1  s'1. 

Key  Words:  air  afterglow,  chemiluminescence,  kinetics,  photometry  standards 
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1.  INTRODUCTION 


The  air  afterglow  has  been  observed  in  discharge  systems  for  more  than  a  century,  and 
has  been  observed  in  the  upper  atmospheric  night  glow  for  at  least  half  that  long.  Kaufman1 
reviewed  most  of  the  pertinent  facts  related  to  the  air  afterglow  about  25  years  ago,  and  more 
recently,  Adler-Golden  proposed  a  model  for  it.2,3  Basically,  mixtures  of  atomic  oxygen  and 
nitric  oxide  combine  in  an  apparent  two-body  reaction  to  emit  a  continuum  radiation  which 
extends  from  about  390  nm  out  to  beyond  3000  nm  in  the  IR.  The  intensity  of  this  emission  is 
directly  proportional  to  the  product  of  the  number  densities  of  atomic  oxygen  and  nitric  oxide, 
and  is  independent  of  pressure  of  bath  gas  for  system  pressures  greater  than  about  0.2  Torr: 

^O/NO  =  ^5/No[0][NO].  (1) 

At  lower  pressures  the  intensity  for  constant  number  densities  of  O  and  NO  drops  off  and 
plateaus  to  a  value  about  one  tenth  the  high  pressure  limit  for  pressures  below  about  a  mTorr. 
Although  in  the  Torr  range,  the  intensity  is  independent  of  pressure,  it  does  vary  to  some  extent 
with  the  nature  of  the  bath  gas.  For  example,  intensities  in  a  helium  bath  gas  are  somewhat 
higher  than  those  in  nitrogen;  oxygen  or  argon,  all  of  which  are  roughly  comparable,  and, 
conversely,  intensities  are  somewhat  lower  in  a  C02  or  SF6  atmosphere. 

Since  the  pioneering  observations  of  Fontijn  et  al.,4  the  air  afterglow  has  been  used  as  a 
quantitative  photometric  standard  in  discharge-flow  systems,  and  the  validity  of  a  large  number 
of  kinetic  measurements  depend  upon  the  basic  accuracy  of  the  air  afterglow  rate  coefficient.  In 
the  20  years  subsequent  to  Fontijn  et  al.'s  publication,  at  least  seven  other  papers  presented  values 
for  the  air  afterglow  rate  coefficient.511  Direct  comparison  of  the  results  from  the  various  groups 
is  difficult  because  the  authors  have  chosen  a  variety  of  formats  in  which  to  present  their  results. 

2 

204 


Although  some  have  presented  an  absolute  value  tied  to  a  particular  wavelength,  most  have 
reported  the  integral  of  the  intensity  over  some  wavelength  interval,  one  which  generally  differs 
from  one  paper  to  another.  When  we  reduced  all  reported  measurements  to  a  value  at  a  common 
wavelength  and  temperature,  we  found  that  reported  values  varied  by  a  factor  of  two  and  a  half. 
Five  of  the  eight  reports,  however,  are  within  ±  15%  of  each  other. 

Reports  of  the  relative  intensity  of  the  air  afterglow  emission  as  a  function  of  wavelength 
are  reasonably  consistent  at  wavelengths  shorter  than  700  nm,4,5,7'12  but  diverge  considerably  at 
longer  wavelengths.4^'8, 10,  u' 13  Our  observations  indicate  that  some  of  the  disagreement  may 
result  from  systematic  errors  in  the  response  calibrations  of  the  spectral  systems  used  for  the 
measurements. 

The  purpose  of  the  present  work  was  to  examine  the  literature  to  place  all  reported 
measurements  on  a  common  basis  so  as  to  allow  us  to  recommend  a  consensus  value  for  the  air 
afterglow  rate  coefficient.  In  addition,  we  have  measured  the  variation  in  the  relative  intensity  of 
the  air  afterglow  emission  as  a  function  of  wavelength  over  the  range  between  400  and  1750  nm. 
Our  measurements  allow  us  to  determine  which  of  the  previous  measurements  are  reliable.  Since 
the  air  afterglow  is  rather  easily  generated,  its  distribution  as  a  function  of  wavelength  has  at 
times  been  used  to  provide  relative  response  calibrations  for  photometric  systems.  Naturally 
using  an  incorrect  distribution  will  result  in  invalid  results.  In  fact  it  was  our  initial  inability  to 
obtain  reasonable  results  for  our  measurements  of  branching  ratios  from  common  upper  vibronic 
levels  to  various  lower  ones  in  the  nitrogen  first  positive  system,14  when  we  used  an  air  afterglow 
derived  spectrometer  response,  which  prompted  this  investigation. 
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2.  EXPERIMENTAL 


Our  measurements  of  the  air  afterglow  spectral  distribution  used  a  2-in.  diameter 
discharge  flow  reactor.  The  apparatus  was  configured  so  that  the  monochromator  looked  through 
a  CaF2  window  down  the  axis  of  the  flow  reactor  without  any  intervening  optics  between  the 
monochromator  entrance  slit  and  the  flow  reactor  window.  This  procedure  eliminates  any 
distortions  to  the  intensity  distribution  that  might  accrue  as  a  result  of  chromatic  aberrations  in 
refractive  optical  elements.  The  reactants  were  prepared  and  mixed  in  a  side  arm  of  the  reactor 
prior  to  their  flowing  into  the  observation  region. 

We  used  a  variety  of  approaches  to  generate  the  air  afterglow.  All  techniques  involved 
generating  at  least  one  of  the  free  radicals  in  a  microwave  discharge.  In  some  instances,  both  free 
radicals  were  generated  in  the  same  discharge.  The  sources  included  discharging  02  dilute  in  Ar 
or  He  and  adding  NO  to  the  discharge  effluent;  discharging  N2  dilute  in  Ar  or  He  and  adding 
excess  NO,  so  as  to  convert  all  the  N  atoms  generated  in  the  discharge  to  O  atoms  with  some 
residual  NO  remaining;  discharging  mixtures  of  02  dilute  in  N2  or  N2  dilute  in  02;  and 
discharging  N20  dilute  in  helium.15  All  sources  gave  the  same  relative  spectral  distribution. 

The  monochromator  used  was  a  0.5  m  focal  length  Minuteman  instrument  using  different 
grating  and  detector  combinations  to  cover  the  wavelength  regions  400  to  880  nm  and  650  to 
1750  nm.  At  the  shorter  wavelengths,  we  used  a  1200  g  mm'1  grating  blazed  at  250  nm  and  a 
thermoelectrically  cooled  HTV  R-943-02  photomultiplier.  Signals  were  processed  using  either 
a  photon-counting  rate  meter  or  picoammeter.  For  the  longer  wavelength  region,  we  used  a 
300  g  mm'1  grating  blazed  at  1000  nm  and  a  liquid  nitrogen  cooled  intrinsic  germanium  detector. 
For  these  measurements,  light  from  the  reactor  was  chopped  prior  to  its  entering  the  mono¬ 
chromator  and  a  phase  sensitive  amplifier  processed  the  signals. 
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The  relative  response  of  the  monochromator  was  calibrated  by  scanning  the  output  of 
several  different  radiation  standards  whose  emission  was  reflected  into  the  monochromator  off  a 
diffusely  reflecting  screen.  For  our  early  measurements  we  used  an  aluminum  plate  coated  with 
several  millimeters  of  BaS04  as  the  diffuse  reflector,  but  later  on  used  a  Spectralon®  target.  We 
used  three  different  quartz  halogen  lamps  with  spectral  calibrations  traceable  to  the  National 
Institute  of  Standards  and  Technology,  and,  in  addition,  a  1000°C  black  body. 

3.  RESULTS 

Relative  Spectral  Distribution 

Figure  1  shows  how  the  intensity  of  the  air  afterglow  varies  as  a  function  of  wavelength 
between  400  and  1750  nm.  This  distribution  results  from  the  average  of  a  number  of  spectra 
taken  over  a  5-yr  time  period  using  a  variety  of  excitation  methods  for  generating  the  air 
afterglow.  Distributions  between  400  and  880  nm  and  650  and  1750  nm  were  normalized  to  give 
common  values  between  700  and  850  nm.  The  relative  agreement  over  this  wavelength  region 
was  within  2%  for  the  two  detector  systems. 

In  general  our  results  for  the  relative  intensity  distribution  in  the  visible  are  reasonably 
consistent  with  previous  measurements.  Ratios  of  the  other  distributions  to  ours  are  constant,  for 
the  most  part,  to  better  than  ±  10%  over  the  wavelength  range  450  to  750  nm,  although  there 
appear  to  be  some  systematic  drifts  in  a  few  of  the  ratios.  We  have  plotted  these  ratios  in 
Figure  2.  A  completely  flat  line  would  correspond  to  perfect  agreement  with  our  results.  Prior  to 
making  the  ratios,  the  distributions  were  all  normalized  to  have  a  common  integral  value  between 
450  and  750  nm.  The  various  data  sets  plotted  in  Figure  2  have  been  offset  from  each  other  for 
clarity  and  error  bars  illustrating  variations  of  ±  5%  have  been  attached  to  some.  Most  notable  in 
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the  figure  is  that  our  results  agree  with  those  of  Paulsen  et  al.12  (curve  2  in  Figure  2)  and  Golomb 
and  Brown7  (curve  1)  to  better  than  ±  2%  over  the  whole  wavelength  region,  and  with  the  results 
of  Woolsey  et  al.8  (curve  3)  and  Bradbum  and  Lilenfeld11  (curve  4)  to  about  ±  5%. 

Our  measurements  indicate  the  peak  in  the  intensity  distribution  is  shifted  slightly  to  the 
red  compared  to  some  determinations.  We  find  the  peak  in  the  distribution  to  be  around  670  nm 
as  compared  to  the  600  nm  reported  by  Fontijn  et  al.4  (curve  8),Vanpee  et  al.5  (curve  6),  Pravilov 
and  Smirnova9  (curve  5),  and  Sutoh  et  al.10  (curve  7).  Our  findings,  however,  are  supported  by 
the  other  groups  reporting  this  distribution.7,  “■ 12 

Figure  3  shows  the  ratio  of  intensity  at  each  wavelength  in  the  IR  region  of  the  spectrum 
as  reported  by  the  various  groups  to  that  which  we  determined  plotted  as  a  function  of  wave¬ 
length.  The  distributions  were  again  normalized  to  have  a  common  integrated  intensity  over 
an  appropriate  wavelength  interval,  and  have  again  been  offset  from  each  other  for  clarity.  We 
see  that  the  ratio  of  the  measurements  of  Vanpee  et  al.5  (curve  4  in  Figure  3)  and  of  Bradbum  and 
Lilenfeld"  (curve  3)  to  ours  are  nearly  constant  over  the  complete  wavelength  range,  showing 
very  good  agreement  between  the  distributions  determined  by  the  three  groups.  The  general  drift 
downward  in  the  curves  related  to  the  data  reported  by  Fontijn  et  al.4  (curve  7),  Sutoh  et  al.10 
(curve  6),  and  Woolsey  et  al.8  (curve  5)  indicate  these  groups  all  under  determined  the  amount  of 
radiation  in  this  region. 

To  the  red  of  1250  nm,  we  plot  ratios  of  data  from  Golde  et  al.6  (curve  2  in  Figure  3)  and 
Stair  and  Kennealy13  (curve  1).  The  fact  that  the  two  curves  are  essentially  parallel  indicates 
excellent  agreement  between  the  results  from  these  two  groups.  Their  results,  however,  indicate 
a  systematic  variation  with  respect  to  ours.  This  variation  might  result  from  some  problem  in 
calibration  units.  We  find  that  if  we  adjust  both  sets  of  data  by  multiplying  by  wavelength,  the 
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intensity  ratios  become  constant  to  about  ±  2%  over  the  common  wavelength  region.  Such  a 
correction  would  obtain  if  the  actual  plotted  distributions  had  been  in  units  of  Watts  rather  than 
photon  per  second. 

Interestingly,  the  data  of  Sutoh  et  al.  also  show  this  trend  to  the  red  of  900  nm.  Their 
distribution  between  450  and  900  nm  is  reasonably  consistent  with  our  measurements,  but  then 
begins  to  diverge  considerably  from  ours  in  the  IR.  If  we  weight  their  published  distribution  to 
the  red  of  900  nm  by  wavelength,  however,  their  distribution  and  ours  become  congruent.  Such  a 
discrepancy  would  be  consistent  with  their  having  changed  calibration  units  for  their  IR  detector. 

While  this  suggestion  of  a  systematic  error  is  appealing,  Golde  et  al.  are  quite  explicit  that 
their  results  are  in  photon  units.  The  distribution  of  Stair  and  Kennealy  was  originally  reported 
to  be  units  of  Watts  per  wavenumber,  and  we  were  careful  to  convert  to  a  consistent  set  of  units. 
The  fact  that  our  results  agree  excellently  with  those  of  Vanpee  et  al.  and  with  Bradbum  and 
Lilenfeld  as  well  as  the  fact  that  our  earlier  work  on  the  transition  moment  variation  of  the 
nitrogen  first  positive  band  system  gave  reasonable  results  only  when  the  relative  monochromator 
response  function  used  to  derive  the  present  data,  was  used,  indicates  to  us  that  our  results  are 
essentially  correct. 

Table  1  lists  our  recommended  relative  distribution  as  a  function  of  wavelength  over  the 
range  between  400  and  1750  nm. 

Absolute  Rate  Coefficient 

The  absolute  rate  coefficient  for  the  air  afterglow  has  been  reported  by  eight  different 
groups.  Unfortunately,  inconsistencies  in  how  the  various  groups  reported  their  results,  as  well 
as  differences  in  the  conditions  employed,  make  comparison  between  the  various  results  difficult. 
We  have  analyzed  each  of  the  results  to  determine  a  consensus  value  for  the  rate  coefficient  at  a 
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single  wavelength.  The  rate  coefficient  at  other  wavelengths,  then,  can  be  determined  by 
multiplying  the  value  at  the  reference  wavelength  by  the  ratio  of  the  intensity  from  the  relative 
distribution  at  the  desired  wavelength  to  that  at  the  reference  wavelength.  We  have  chosen  our 
reference  wavelength  to  be  580  nm.  This  wavelength  corresponds  to  a  peak  in  the  nitrogen-first 
positive  radiation  when  it  is  generated  by  N-atom  recombination.  It  is  one  of  the  wavelengths 
most  convenient  to  use  for  calibrating  optical  systems  when  using  the  air  afterglow  technique. 
We  have  described  our  procedures  for  doing  this  in  some  detail  previously.15'17 

Fontijn  et  al.4  published  the  first  measurements  on  the  absolute  photon  yield  for  the  air 
afterglow  reaction.  They  used  an  actinometric  technique  whereby  light  from  the  air  afterglow 
was  sent  through  a  cell  containing  a  solution  of  K3Fe(C204)3.  The  light  from  the  reaction 
converted  the  ferric  ion  in  solution  to  ferrous  ion,  and  the  final  concentration  of  ferrous  ion  was 
subsequently  determined  spectrophotometrically.  The  quantum  yield  for  the  conversion  process 
as  a  function  of  wavelength  had  been  determined  previously.  These  measurements  along  with 
their  determination  of  the  distribution  of  intensity  of  the  air  afterglow  as  a  function  of  wavelength 
led  them  to  report  the  total  air  afterglow  rate  coefficient,  integrated  over  the  extent  of  their 
spectral  distribution,  measured  by  them  to  be  385  to  1400  nm,  as  6.4  x  10'17  cm3  molecule'1  s'1. 
Since  subsequent  observations  indicated  that  they  seriously  underestimated  the  intensity  of  the  air 
afterglow  emission  at  longer  wavelengths,  this  number  is  generally  considered  to  be  low. 

Their  actinometric  measurements  were  sensitive  to  the  air  afterglow  only  from  its  onset 
at  390  nm  out  to  about  510  nm,  where  the  actinometric  quantum  yield  drops  off.  Thus  one  can 
expect  reasonable  accuracy  for  their  value  at  shorter  wavelengths.  By  multiplying  their  total 
rate  coefficient  by  the  ratio  of  their  relative  intensity  in  a  one  nanometer  interval  around 
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580  nm  to  the  integral  over  the  intensity  distribution  they  published,  we  find  that  their 
measurements  indicate  a  value  for  the  air  afterglow  rate  coefficient  at  580  nm  of  1.36  x  10'19  cm3 
molecule'1  s'1  nm'1. 

Vanpee  et  al.5  investigated  the  air  afterglow  in  a  free  jet  expansion  of  O  and  NO.  They 
characterized  the  jet’s  physical  characteristics  and  used  a  photometric  system  they  had  calibrated 
against  several  tungsten  ribbon  lamps.  They  also  demonstrated  that  the  O/NO  intensity 
distribution  extended  out  to  at  least  2000  nm.  Their  published  result  was  6.8  x  10'17  cm3 
molecule'*  s'1  over  the  wavelength  range  between  385  and  2000  nm  and  at  a  temperature  of 
367  K.  In  order  to  compare  their  results  with  others,  one  must  not  only  correct  for  a  different 
wavelength  range,  but  also  for  a  different  temperature.  Golomb  and  Brown7  report  a  temperature 
dependence  which  incorporates  not  only  their  own  measurements  but  also  those  of  several  other 
groups.  Using  this  value  to  correct  Vanpee  et  al.’s  measurements  to  300  K,  and  further 
correcting  the  Vanpee  et  al.  value  to  our  reference  wavelength,  results  in  a  rate  coefficient  at 
580  nm  of  1.59  x  10'19  cm3  molecule'1  s'1  nm"1. 

Golde  et  al.6  investigated  the  air  afterglow  in  the  IR  between  1280  and  3300  nm.  They 
calibrated  their  optical  system  absolutely  against  a  black  body  source  and  determined  that  the 
rate  coefficient  for  the  air  afterglow  emission  at  1510  nm  was  2.4  x  10"2°  cm3  molecule'1  s'1  nm"1. 
If  we  multiply  this  value  by  the  ratio  of  intensity  at  580  nm  to  that  at  1510  nm,  as  determined  in 
this  work,  we  find  a  value  for  the  air  afterglow  rate  coefficient  at  580  nm  of  1.16  x  10'19  cm3 
molecule'1  s'1  nm'1. 

Golomb  and  Brown7  conducted  an  extensive  series  of  measurements  on  the  air  afterglow 
emission  rate  coefficient  as  a  function  of  temperature  and  pressure.  They  used  standard  radiance 
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sources  to  calibrate  their  optical  system.  They  reported  a  rate  coefficient  of  3.7  x  10'17  cm3 
molecule*1  s*1  over  the  wavelength  range  between  400  and  800  nm  at  300  K.  This  translates  to  a 
value  of  1.14  x  10'19  cm3  molecule'1  s'1  nm'1  at  580  nm. 

Woolsey  et  al.8  developed  a  novel  radiation  calibration  source  which  mimicked  the 
geometry  one  would  observe  when  looking  along  the  axis  of  a  discharge  flow  reactor.18  They 
used  this  to  calibrate  their  optical  system  and  subsequently  measured  rate  coefficients  for  the  air 
afterglow  emission  rate  at  a  number  of  wavelengths  between  530  and  1400  nm.  Their  results 
give  a  value  of  0.72  x  10'19  cm3  molecule'1  s'1  nm'1  at  580  nm. 

Pravilov  and  Smirnova9  measured  the  absolute  air  afterglow  rate  coefficient  between 
430  and  630  nm,  but  reported  a  value  of  8.0  x  10'17  cm3  molecule'1  s'1  over  the  wavelength 
range  between  395  and  2000  nm  at  300  K.  They  reported  a  relative  intensity  distribution 
between  400  and  about  950  nm,  and  extended  the  distribution  into  the  1R  using  the  distribution 
Vanpee  et  al.  had  published.  Using  their  composite  curve,  we  infer  a  value  from  their  data  of 
1.27  x  10'19  cm3  molecule'1  s'1  nm'1  at  580  nm.  This  measurement,  however,  was  made  in  a 
helium  buffer  gas  which  has  a  slightly  higher  third  body  efficiency  than  does  oxygen,  nitrogen,  or 
argon,  all  of  which  produce  similar  air  afterglow  intensities  given  comparable  concentrations  of 
O  and  NO  in  one’s  reactor.  Kaufman’s  review1  indicates  an  enhancement  over  argon  as  a  third 
body  of  1 .21  ±  0.04.  Making  this  correction  leads  to  a  value  for  the  air  afterglow  rate  coefficient 
of  1.05  x  10"19cm3  molecule'1  s'1  nm'1. 

Sutoh  et  al.10  measured  the  absolute  value  of  the  air  afterglow  rate  coefficient  at  562  nm 
using  a  photometer  for  a  detector  (essentially  an  interference  filter/photomultiplier  combination 
with  a  carefully  controlled  field  of  view).  They  calibrated  their  photometer  using  an  absolute 


10 

2i2 


brightness  standard.  Their  result  was  1.07  x  10'19  cm3  molecule'1  s'1  nm'1  at  562  nm  which 
translates  to  a  value  of  1.14  x  10'19  cm3  molecule'1  s'1  nm'1  at  580  nm. 

Finally  Bradbum  and  Lilenfeld11  compared  the  intensity  of  the  air  afterglow  emission  at 
1270  nm  produced  by  known  concentrations  of  O  and  NO  with  the  intensity  of  emission  from  02 
(a  *Ag)  when  known  concentrations  of  singlet  oxygen  were  in  their  reactor.  They  measured  the 
singlet  oxygen  via  EPR  for  which  an  absolute  calibration  could  be  determined  by  comparing 
singlet  oxygen  signals  with  those  produced  by  known  concentrations  of  normal,  ground  state, 
triplet  oxygen.  Using  the  widely  accepted  value  for  the  radiative  lifetime  for  02  (a  lAg)  of 
64  min,19  they  could  place  their  fluorescence  measurements  on  an  absolute  basis.  Their  result, 
0.62  x  10"19  cm3  molecule'1  s'1  nm'1  at  1270  nm,  translates  to  a  value  of  1.86  x  10'19  cm3 
molecule'1  s'1  nm'1  at  580  nm.  The  accuracy  of  this  value,  naturally,  depends  upon  the  accuracy 
of  the  lifetime  for  02  (a  ’Ag).  We  discuss  this  issue  below. 

Table  2  summarizes  the  values  for  the  air  afterglow  rate  coefficient  at  580  nm  and  room 
temperature  we  have  determined  from  analyzing  the  results  of  each  group.  The  data  divide  into 
roughly  three  ensembles.  Five  of  the  groups  obtained  results  which  are  remarkably  consistent, 
two  reports  are  considerably  higher  than  the  bulk  of  the  measurements,  while  Woolsey  et  al.’s 
determination  is  much  lower  than  any  of  the  other  reported  values.  An  average  of  the  results  of 
Fontijn,  et  al.,  Golde  et  al.,  Golomb  and  Brown,  Pravilov  and  Smimova,  and  Sutoh  et  al.  yields  a 
value  of  (1.17  ±  0.1 1)  x  10'19  cm3  molecule'1  s'1  nm'1  at  580  nm,  where  the  error  bar  represents 
one  standard  deviation.  All  of  the  measurements  included  in  the  average  are  within  two  standard 
deviations  from  the  mean,  whereas  the  other  three  values  are  all  four  or  more  standard  deviations 
removed  from  this  mean  value. 
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4.  DISCUSSION 


Most  of  the  determinations  of  the  absolute  rate  coefficient  for  the  air  afterglow  depend 
upon  calibrations  determined  by  comparison  with  standard  lamps.  Such  calibrations  are  straight 
forward  in  principle,  but  are  often  difficult  to  implement  in  practice.  The  reason  for  this  is  that 
the  geometry  of  the  optical  system  is  generally  quite  different  when  viewing  an  extended  source 
such  as  the  air  afterglow  compared  to  a  compact  filament  of  a  standard  lamp.  The  generally  good 
agreement  between  most  of  the  groups  using  such  calibration  procedures  indicates  their 
geometrical  corrections  have  been  handled  correctly. 

Of  the  two  measurements  which  do  not  use  direct  comparison  with  a  standard  photo¬ 
metric  source,  one  agrees  well  with  the  bulk  of  the  observations  while  the  other  does  not. 

Fontijn  et  al.’s  determination  was,  in  effect,  subsequently  confirmed  by  Fontijn  and  Lee20  using 
a  second  chemiluminescent  comparison.  That  is  to  say,  they  “calibrated”  an  optical  system  using 
the  air  afterglow  and  using  the  phosphorescence  from  the  liquid  luminol  reaction,  whose 
quantum  yield  had  been  established  previously  by  several  groups,  and  obtained  calibration 
factors  agreeing  within  10%. 

As  we  mentioned  above,  the  determination  of  Bradbum  and  Lilenfeld  depends  on  the 
validity  of  the  Einstein  coefficient  for  spontaneous  emission  from  02  (a  ‘Ag).  For  30  years  the 
value  reported  by  Badger  et  al.,  2.58  x  104  s'1,  was  accepted  without  serious  question,  although 
several  estimates  based  on  aeronomic  observations  tended  to  favor  a  somewhat  lower  value.21 
The  remeasurement  of  this  value  by  Lin  et  al.,22  which  was  subsequently  superceded  by  a  more 
careful  measurement  from  the  same  group  by  Hsu  et  al.,23  appeared  to  confirm  Badger  et  al.’s 
measurement.  Mlynczak  and  Nesbitt,21  however,  showed  that  Lee  et  al.  and  Hsu  et  al.  converted 
their  line  strength  data  into  a  band  Einstein  coefficient  incorrectly  when  they  used  the  wrong 
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value  for  a  degeneracy  ratio  between  ground  and  excited  states  of  the  transition.  When  the 
proper  degeneracy  ratio  is  used,  the  Einstein  coefficient  inferred  from  Hsu  et  al.’s  data  is  actually 
1.47  x  W4  s'1. 

Using  this  corrected  value  for  the  02  (a  *Ag)  Einstein  coefficient  to  analyze  the  measure¬ 
ments  of  Bradbum  and  Lilenfeld  would  result  in  a  value  of  the  air  afterglow  rate  coefficient  at 
580  nm  of  1.04  x  10'19  cm3  molecule'1  s'1  nm'1,  in  agreement  with  the  bulk  of  the  measurements 
listed  in  Table  1.  Alternatively,  if  one  considers  the  value  for  the  air  afterglow  rate  coefficient 
we  propose  here  to  be  correct,  then  one  could  look  upon  Bradbum  and  Lilenfeld’s  work  as 
essentially  an  alternative  measurement  of  the  Einstein  coefficient  for  spontaneous  radiation  from 
02  (a  *Ag).  One  would  then  derive  a  value  of  1.65  x  1CT4  s'1  for  this  transition  probability, 
agreeing  within  12%  with  the  re-analyzed  Hsu  et  al.  value,  and,  thereby,  confirming  Mlynczak 
and  Nesbitt’s  contention  that  the  long  accepted  value  of  Badger  et  al.  is  indeed  incorrect.  This 
issue  is  still  a  contentious  one  in  the  aeronomic  community.24 

The  air  afterglow  rate  coefficient  integrated  over  the  wavelength  range  studied  here, 

400  to  1750  nm,  is  8.15  x  10'17  cm3  molecule'1  s'1.  If  we  normalize  the  relative  distribution 
reported  by  Golde  et  al.  to  outs  at  1650  nm,  we  can  integrate  the  rate  coefficient  out  to  3500  nm. 
This  results  in  a  total  rate  coefficient  for  emission  in  the  recombination  of  O  with  NO  of 
8.9  x  10'17  cm3  molecule'1  s'1.  This  value  agrees  excellently  with  the  9.4  x  10'17  cm3  molecule'1  s'1 
proposed  by  Golde  et  al.  from  a  combination  of  their  results  in  the  IR  with  those  of  Fontijn  et  al. 
in  the  visible. 
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Table  1 .  Distribution  of  Air  Afterglow  Intensities  as  a  Function  of  Wavelength 


Wavelength 

Relative 

Absolute  Rate  Coefficient 

(nm) 

Intensity 

(10'19  cm3  molecule'1  s'1  nm'1) 

400 

1.3 

410 

4.1 

420 

7.4 

0.09 

430 

12.6 

0.15 

440 

— 

■■■■ 

450 

22.5 

0.28 

460 

29.0 

0.36 

470 

36.5 

0.45 

480 

44.2 

0.54 

490 

50.1 

0.62 

500 

56.8 

0.70 

510 

63.1 

0.78 

520 

70.4 

0.87 

530 

74.4 

0.92 

540 

80.4 

0.99 

550 

84.5 

1.04 

560 

88.5 

1.09 

570 

92.0 

1.13 

580 

95.1 

1.17 

590 

97.8 

1.20 

600 

100.0 

1.23 

610 

101.9 

1.25 

620 

103.4 

1.27 

630 

104.6 

1.29 

640 

105.4 

1.30 

650 

105.9 

1.30 

660 

106.2 

1.31 

670 

106.2 

1.31 

680 

105.9 

1.30 

690 

105.4 

1.30 

700 

104.7 

1.29 

710 

103.8 

1.28 

720 

102.7 

1.26 

730 

101.5 

1.25 

16 

218 


Table  1.  Distribution  of  Air  Afterglow  Intensities  as  a  Function  of  Wavelength  (Continued) 


Wavelength 

Relative 

Absolute  Rate  Coefficient 

(nm) 

■ESESm 

(10'19  cm3  molecule’1  s’1  nm'1) 

740 

1.23 

750 

98.6 

1.21 

760 

96.6 

1.19 

770 

94.6 

1.16 

780 

92.5 

1.14 

790 

90.4 

1.11 

88.3 

1.09 

810 

86.1 

1.06 

820 

84.0 

1.03 

830 

81.9 

1.01 

840 

79.7 

0.98 

850 

77.6 

0.95 

860 

75.5 

0.93 

870 

73.4 

0.90 

880 

72.5 

0.89 

890 

70.9 

0.87 

69.2 

0.85 

950 

61.7 

0.76 

1000 

55.2 

0.68 

1050 

49.5 

0.61 

1100 

44.6 

0.55 

1150 

40.3 

0.50 

36.5 

0.45 

1250 

33.2 

0.41 

1300 

30.3 

0.37 

1350 

27.6 

0.34 

1400 

25.1 

0.31 

1450 

22.6 

0.28 

1500 

20.2 

0.25 

1550 

17.7 

0.22 

1600 

15.3 

0.19 

1650 

14.3 

0.18 

1700 

13.1 

0.16 

1750 

12.2 

0.15 
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Table  2.  Determinations  of  the  Absolute  Air  Afterglow  Rate  Coefficient 

at  580  nm  and  300  K 


Group 

Year 

Absolute  Rate  Coefficient 
(1019  cm3  molecule'1  s'1  nm'1) 

Fontijn  et  al.4 

1964 

1.36 

Vanpee  et  al.5 

1971 

1.59 

Golde  et  al.6 

1973 

1.16 

Golomb  and  Brown7 

1975 

1.14 

Woolsey  et  al.8 

1977 

0.72 

Pravilov  and  Smirnova9 

1978 

1.05 

Sutoh  et  al.10 

1980 

1.14 

Bradbum  and  Lilenfeld11 

1988 

1.86 
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Figure  Captions 


1.  Relative  distribution  of  the  intensity  of  the  air  afterglow  as  a  function  of  wavelength. 

2.  Ratio  of  the  air  afterglow  intensity  distributions  between  420  and  900  nm  reported  in  the 
literature  to  that  determined  in  this  work.  The  curves  have  been  offset  from  each  other 
for  clarity.  The  bars  on  some  of  the  distributions  represent  variations  of  ±  5%  in  the 
ratios.  The  numbered  curves  are  as  follows:  1  -  Golomb  and  Brown7;  2  -  Paulsen  et  al.12; 
3  -  Woolsey  et  al.8;  4  -  Bradbum  and  Lilenfeld11;  5  -  Pravilov  and  Smirnova9;  6  -  Vanpee 
et  al.5;  7  -  Sutoh  et  al.10;  and  8  -  Fontijn  et  al.4 

3.  Ratio  of  the  air  afterglow  intensity  distributions  between  800  and  1800  nm  reported  in  the 
literature  to  that  determined  in  this  work.  The  curves  have  been  offset  from  each  other 
for  clarity.  The  bars  on  some  of  the  distributions  represent  variations  of  ±  5%  in  the 
ratios.  The  numbered  curves  are  as  follows:  1  -  Stair  and  Kennealy13;  2  -  Golde  et  al.6; 

3  -  Bradbum  and  Lilenfeld11;  4  -  Vanpee  et  al.5;  5  -  Woolsey  et  al.8;  6  -  Sutoh  et  al.10; 
and  7  -  Fontijn  et  al.4 


19 

221 


APPENDIX  9 


CIRRIS  1A  Ozone  Airelow  Analysis 


225/226 


VG93-194 


OBJECTIVES 


C 

0 

0  Q 


£  CL 
0  0 


0  0  CO 

■J=  0 


li=  O 
O  0 


0  Q 

O  0 
CD  C 


0  S 
o  o 


>,  0  iS 
£  'g  o 
c  0  x 
0  QC  LU 


0  0 
c  CL 

E  E 
c  o 

0  O 

0  i 

Q 


228 


Improve  profiles,  distributions,  and  relaxation  rates  to  permit 
model  refinement 
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Preliminary  comparison  with  SPIRE  and  model 
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CONTINUING  ANALYSIS  EFFORTS 
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Provide  profiles,  distributions,  and  rates  to  permit  predictive 
atmospheric  radiance  model  refinement 


